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Abstract 



The mass mj of the top quark t is now measured to be 174 ±5 GeV. Improvement of the present 
error Amt — 5 GeV is important not only for its own sake. For example with Arrit ~ 0.2 GeV, 
the Higgs mass ttih is predicted from electroweak precision measurements with an error of 
AmH/TTiH ~ 20%, which is now ~ 60%. There is a consensus that rrit is best measured at 
future e'^e~ Linear Colliders (LC) by using the line shape of the total cross section cTtot(v^ ~ 
2mt; e+e" —>■ tt) near the ti threshold. However, it was reported recently that the Next-to- 
Next-to-Leading-Order (NNLO) correction to cxtot is as large as the NLO correction to it. This 
imphes that the convergence of the series is not good, and this theoretical uncertainty affects 
the precision of the determination of rrit at LC. In the thesis, we improve the convergence 
of the perturbative series of (Xtot by choosing appropriate expansion parameters for both nit 
and the strong gauge coupling a^. As for the mass parameter, the pole mass was used in the 
previous studies. However, because of the infrared (IR) structure of QCD, the pole mass of 
a quark is defined only with an ambiguity of C(Aqcd) ~ 1 GeV. Instead, we use another 
mass parameter called potential-subtracted mass, which is less affected by the IR region of 
QCD. As for the gauge coupling, we resum a certain class of higher order corrections by using 
renormalization group. With these two prescriptions, the perturbative convergence of Utot 
is much improved. Correspondingly, the resulting ambiguity for the determination of rrit is 
estimated to be ~ 0.1 GeV, which is not larger than the expected statistical error at future 
LC. We also study the top-quark momentum distribution da/dpt. The NNLO correction for 
the line-shape of da/dpt turns out to be fairly small compared to that for citot- 

Aside from the contribution of Standard Model (SM), we also study that of new physics. 
Since Cabbibo-Kobayashi-Maskawa matrix elements for the top quark are almost diagonal 
{\Vtb\ ^ 1), CP- violating effect in top-quark sector is highly suppressed in SM. In the thesis, 
we study the effect of all types of anomalous CP violating top-quark Electric Dipole Moment 
(EDM) interactions near the threshold for the first time. If they are measured to be non-zero, 
it would immediately imply new physics. We show that one can disentangle tig-, tij- and 
tiZ-EDMs by the dependence of CP-odd observables on the CM-energy and on the degree of 
initial state polarization. The sensitivity to chromo-EDM at LC turns out to be better than 
that at LC in the open top region {^/s ^ ^rrit), but to be worse than that at LHC by a factor 
(9(1/10). The sensitivities to electroweak-EDMs near the threshold are similar to those in the 
open top region at LC. 
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Chapter 1 
Introduction 



1.1 Top quark 

As more experimental data is accumulated, it has become clearer that the Standard Model 
(SM) describes the physics around and below the ElectroWeak (EW) scale very well. One of 
the great triumph of the SM is that the top quark t is discovered at Tevatron |]l| within the 
range predicted from EW precision measurements. Currently 



For the indirect fit, the central value and the first uncertainty are for niH = rnz, and the second 
is the shift from changing rriH to 300 GeV. Here mn is the Higgs mass. 

Despite of its success, few people consider the SM as the ultimate theory because of its 
"un-naturalness" . There are several issues and possible solutions to them. First, the gauge 
symmetry of SM is a direct product of three distinct groups, and gauge anomaly cancellation 
is highly non-trivial. These "un-natural" aspects are explained if physics at higher energies is 
governed by a Grand Unified Theory (GUT). In fact, the three gauge couplings unify around 
10^^ GeV in the minimal supersymmetric SM, if the soft SUSY breaking scale is around 1 TeV. 
Weak-scale Super Symmetry (SUSY) itself is motivated to stabilize the EW scale against 
radiative corrections, which is also an "un-natural" aspect of the SM. 

Another one is flavor repetition, which seems redundant, and its breaking by masses and 
mixings. In this respect, top quark is highly exceptional, since the mass of the next heaviest 
fermion b is m?, = 4.1-4.4 GeV 0, which is 40 times smaller than m^. Thus one can expect 
that top quark may play an important role for the physics around TeV scale, especially for 
flavor and Higgs physics. In fact, among the ingredients of the SM, the mechanism that breaks 
the electroweak symmetry is the only one that is not verified experimentally. The situation is 
similar also for top quark. Except its mass rrit, almost no property of the top quark is known 
experimentally. Thus there are still possibilities of its "anomalous" property. For example, it 
might not be merely a "heavy quark", which means that the generation is not repeating in 
fact; it might be a composite of subquarks, or it might condensate to trigger EW symmetry 
breaking which means that there is a certain strong attractive force between tt. Or there 
might be a substantial mixing with a "fourth generation quark" , which can be vector-like. On 
the other hand, even if top quark itself were merely a "heavy quark", it is expected to be 
sensitive to the mechanism of SI]{2)l x U(l)y breaking, because of its large mass and the large 
mass difference between t and b quarks. 



mt = 174.3 ± 5.1 GeV 
mt = 170 ± 7(+14) GeV 



from direct measurement, 
from indirect SM EW fit. 
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Some properties of top quark will be measured at LHC, which operates before e+e~ linear 
colliders does. With its high CM energy and high luminosity, it is expected that certain clue 
to new physics will be discovered at LHC. However precise measurement at hadron colliders 
may not be possible, since the CM energy of each elementary process is not fixed. Thus it is 
lepton colliders that can measure the properties of top quark etc. precisely, after the discovery 
at Tevatron, or hadron colliders in general. 



1.1.1 Top quark pair production in e+e collisions near threshold 

Studies on top quark at e+e~ colliders can be classified into two groups depending on the CM 
energy ^/s] one is at the threshold region {^/s ~ 2mt), and the other is at the open top region 
(a/s ^ 2mt). One prominent feature of the threshold region, which makes the distinction 
above, is multiple gluon rescattering between t and i, which is effectively not suppressed by 
powers of a^- In fact, top quark is an ideal probe of the short distance behavior of QCD. This 
is because of its large mass rrit and its large decay width Ft, which are related as follows: 

r, = ^^1.5GeV. 

This expression holds at tree level in the limit mw/frit 0, and is for the SM, where almost 
every t decays into hW^ . Here \Vth\ — 1 is used, which is indeed the case when the unitarity 
of VcKM among the first three generations is assumed. When the distance between t and 
t is shorter than the typical non-perturbative scale ~ 1/Aqcd ~ 1/(1 GeV), gluon exchanges 
between them can be approximated by the Coulomb potential| V = —Cpag/r, where Cp = 4/3 
is the quadratic Casimir for the fundamental representation of SU(3). Due to the large mass, 
the Bohr radiusQ of the tt pair is much smaller than the typical non-perturbative scale: 
tb = l/{Cpasmt/2) ^ 1/(20 GeV) ^ 1/Aqcd- This means that the QCD potential between 
t and i is surely well within the perturbative region. And due to the large decay width, 
a ti pair decays before it hadronizes that is, the fluctuation of top quark momentum^ 
y/nitTt — 15 GeV S> Aqcd due to the off-shell effect acts as an IR cut-off. Thus a tt pair 
produced in e~^e~ collisions can be treated perturbatively from the production to the decay, 
even near the threshold, where the QCD effect is enhanced. Thus the threshold region is one 
of the best places for the study of perturbative QCD. 

There are many studies on ti pair production near the threshold in e~^e~ collisions P,P|-P^, 
including ours. Some of them are reviewed briefly in the next two sections, and in the former 
halves of Chapters ^ and ^ at some length. Concisely speaking, we (i) reduce the theoretical 
uncertainty for the determination of the top-quark mass using the process e~^e~ — > tt near the 
threshold, and (ii) study the sensitivity of e+e~ colliders with ^/s ~ 2mt to all the anomalous 
Electric Dipole Moments in t-t-gauge boson interactions for the first time. 



1.2 Overview of tt production cross section 



-CF)as/r, where Ca = 3. 



^ More detailed expression is shown in Section B.l. 

^ This expression is for a colo r-singlet qq. For a co lor-octet state, it is F = 
^ It is defined in Section B.3, See also Figure 2.4, 

^ The peak momentum of the top quark is ~ \ \/ 'mt{^/s — 2nit + 1 GeV -I- iTt)\- See Figure 19. The fluctu- 



ation given here is for near the threshold. 



1.2 Overview of tt production cross section 



1.2.1 LO and NLO calculations of dtot ^^id da /dp 

In e'^e~ collisions with ^/s ~ 2mq where q = c,b, there are sharp resonances in the total 
production cross section crtot(s) due to the qq bound states such as J/\I' or T. The locations of 
those bound state poles are used to probe the QCD potential at r ~ l/{asmg), which is near 
the boundary of the perturbative description of QCD . On the other hand, due to the large 
decay width of top quark Tf, (would-be) ti bound-states merge into a single slight rise of o"tot(s) 
near the IS' resonance. Thus in order to calculate crtot(s) just below the threshold {y/s < 2mt) 
reliablyQ, one needs to sum over a host of resonances. This summation is conveniently done 
using the Green function for a tt pair M : 



G{r,r') = (r 



1 



H-{E + tTt) 



r 



where H is the Non-Relativistic (NR) Hamiltonian for tt, E = ^/s — 2mt is the NR energy, and 
ipn and ipk energy eigenfunctions for En{< 0) and E^ = k'^/mt{> 0), respectively. Since 
top quarks are non-relativistic near the threshold, it is convenient to use a non-relativistic 
approximation at the leading order. There is another advantage for the NR Green-function 
treatment. Near the threshold, it is known that the naive perturbative expansion in as breaks 
down ("threshold singularity"). In fact, n-th power of as is accompanied by l/Pt"^ {m < n), 
where (3t is the velocity of t and t at their CM frame. The situation is similar to certain pertur- 
bation series where powers of are accompanied by powers of log(g/yu); a stable perturbative 
expansion is obtained by summing up (potentially large) leading log's. In other words, one 
should identify J2ni^s{fJ') ln((3'//i))" to be the leading order. Likewise, one should identify 
^„(«s/A)" to be the leading order near the threshold. On the other hand, this summation 
is realized^ automatically in the Green function for Coulomb potential between t and i. Thus 
the Green function method provides an efficient way (i) to sum up infinite number of broad 
resonances, and (ii) to sum up powers of as/Pt- After these summations, the cross section is 
finite[] even at the threshold, which correspond^ to (3t = 0. Order counting becomes easier by 
noting as = /{'i'nhc) = 0(l/c) and Pt = Vt/c = 0(l/c), which means ag/Pt = C'(l). This is 
the Leading Order (LO). Likewise the Next-to-Leading Order (NLO) is 0{l/c), and so on. 

The total cross section at the leading order was calculated in [|ri|] with the QCD potential 
V{r) that is Renormalization Group (RG) improved in the coordinate space. It triggered several 
studies both by experimentalists and by theorists [|T^-T^]. A part of NLO corrections was taken 



into account in their analysis. After a while, complete NLO corrections to tt production and 
their decay cross section were calculated [|TB],0. They included the Final-State Interactions 
(FSIs), by which we mean gluon exchange between t and b, t and b, and b and b. It was also 



^ The production cross section grows rapidly near the threshold, because the leading contribution is S'-wave. 
The mass of top quark is determined by measuring this rise of CTtot- Thus the precise calculation of CTtot near the 
threshold is required for this reason. Note also that the "binding energy" for the IS* state is ^ p^/mt ^ 2 GeV, 



where = 1/''b is Bohr momentum. See Section B.3. 

^ The explicit expression is given in Eq. ( 2.46 ). Note that ^Ejmt = A + O{0f)- 
See Figure IJ. 

® Due to finite decay width Fj, velocity (it of to p qu ark is not determined by the energy. Its typical fluctuation 
is /3t ~ ^jYt/mt ~ 0.1 ~ Us- See just below Eq. ( |2.12 ). This is consistent with our order counting. The symbol 
c is speed of light. 
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shown that ahhough these FSIs modify the momentum distribution of t, their corrections to 
the total cross section cancel to NLO. 

Top-anti-top pair production near the threshold in e~^e~ collisions not only provides the 
place for a precise study of perturbative QCD, but also was known to be the best environment 
to determine the top quark mass m^. This is because Monte Carlo detector simulations [0,|,18 



based on NLO calculations |]I1|-[T^ showed that can be determined with statistical error^ 
(ArritY^'^^ — 0.2 GeV using the energy dependence of i?-ratio and the top quark momentum 
distribution da /dp. It was also shown that QCD gauge coupling ^^(m^) can be measured with 
statistical error 0.005. It was assumed that 11 energy points are sampled with 1 fb~^ for each 
point. A similar MC study is shown in [Q. These MC studies are based on Next-to-Leading 
Order (NLO = 0{l/c)) theoretical calculations and assume no theoretical uncertainties. 

1.2.2 NNLO calculations of dtot 

However recently the NNLO corrections to the total cross section atot has been calculated by 
several groups pBI-p^F^ and it turned out to be unexpectedly large. In fact the magnitude of 



the NNLO correction is similar to that of the NLO correction. This is problematic since it is 
desirable to determine the top-quark mass in high precision in order to pin down contributions 
from new physics (and/or Higgs). For example [Q for wide range of Higgs mass tuh, the 
uncertainty (Am/^/m^^) of Higgs mass extracted from radiative corrections depends on the 
uncertainty of rrit, and is 57% for the present precision, while it reduces to 17% for Linear 
Colliders (LCs) assuming {AmtY^'^^ ~ 0.4 GeV. 

1.2.3 Theoretical progress related to dtot 

After the calculations of the NLO corrections to citot, several pieces of theoretical progress have 
been made with the precise calculation for heavy quark- antiquark system One is the use 



of low energy effective theories called NRQCD and pNRQCD [^,|^, which enabled 

the calculations of the NNLO corrections explained above, and the other is the cancellation of 
the ambiguities in the pole mass of a quark and in the coordinate space QCD potential; those 
ambiguities are due to the infrared structure of QCD, which is called "renormalon ambigu- 
ity" |55-|60[]. The latter enables us to improve the perturbative convergence of atot, which is 



made in p9| , on which a part of the thesis is based. Here we explain these two topics briefly. 

At NNLO and beyond, non-relativistic calculations become complicated. There are several 
reasons. Typical one is that the normalization of the ti current ju depends on the momentum 
transfer; in other words, anomalous dimension of ju is non-zero. This makes the use of low- 
energy effective theories particularly powerful. Schematically, the reduction goes as follows. 
The starting point is QCD, which is fully relativistic. By expanding the fields for non-relativistic 
particles in 1/c, one obtains NRQCD [^,|3^. At NNLO (and lower orders) further reduction is 
possible, since no real gluon is emitted^ at these orders. Thus "soft" gluons can be integrated 
out to result in a ti potential [^,0. It is known that this formalism, called potential NRQCD 
or pNRQCD, reproduces correctly the energy shift of the hydrogen atom and positronium to 



10 



9 On the other hand, LHC may give (Amt]f^~ 2 GeV. 
Their results are reproduced in Section [2.6.2 . Comprehensive summary (including our work) to the end 



of 1999 is available in |g3[. 

This may be understood by the fact that the vertex for transverse gluon gx in Coulomb gauge is suppressed 
by 0{l/c) compared to the vertex for Coulomb gluon gc- 
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ma"^ X For actual calculations, one does not need to know the precise "normalization", 

or the matching coefficients, of the each operator in pNRQCD Lagrangian, a priori. Only 
a combination of them is relevant for, say, a total production cross section citot- It can be 
determined by calculating atot in both perturbative QCD and pNRQCD, and by demanding 
those two results coincide. This procedure called "direct matching" |^ is possible^ since for 
Q^s ^ A ^ 1 both perturbative QCD and pNRQCD are applicable. 

On the other hand, "renormalon ambiguity" |55-^ is related to the convergence of a 
perturbative expansion in QFT, which is in fact an asymptotic expansion. One of the techniques 
to sum up an asymptotic series S{a) with respect to a is Borel summation]^, where the sum 
is expressed in terms of integration of e"*/"- S(t) from t = to t = oo. Here S{t) is the 
function related to S{a). The fact that S{a) diverges is expressed by the poles of S(t) on the 
integration path t = O-oo, which are called "renormalon poles". One can step aside those poles 
by modifying the integration path; there may be two ways for each pole. Ambiguities due to 
those poles may be estimated by the difference between the two choices of the modification. 
Due to the factor e~*/", the severest ambiguity comes from the "renormalon pole" nearest to 
the origin [t = 0). Now, it is known that the pole mass mpoie suffers from the renormalon 
poles which originate from the IR structure of QCD. However it was shown [^, |60|] 

that in the combination 2mpoie + V{r), where V{r) is the coordinate space QCD potential, the 
severest renormalon pole is cancelled; this means V{r) also suffer from the renormalon pole, 
and the location of the pole is the same as for mpoie and the residue is twice and the sign is 
opposite. This suggests that one should not attempt to determine pole mass. Instead another 
mass scheme that do not sensitive to IR physics of QCD, such as MS mass, should be used. 
Several such "short-distance" mass schemes are proposed in literature; potential-subtracted 



mass 



591, IS* mass [ETI], and kinetic mass |62 



1.2.4 Our contributions 

As was explained above, the previous NNLO calculations of atot showed large corrections. How- 
ever, since the convergence of a perturbative series changes depending on the renormalization 
scheme^, there is a possibility for finding a scheme where cxtot converges better. This is the 
subject of the latter half of Chapter |[ There we implement two prescriptions; one is to use 
the mass scheme that do not suffer from IR renormalon ambiguity; and the other is to sum up 
(potentially large) leading logarithms by using renormalization group. Both of these reduce the 
theoretical uncertainty of cxtot- With this improvement, we estimate the theoretical uncertainty 
[AmtY^ of the top quark mass determination is reduced to ~ 0.1 GeV, which is smaller than 
the expected statistical error. 

We also calculate the NNLO correction for da /dp from the rescattering between t and i. 
However it is not considered that the rescattering between t and b etc., which may modify the 
distribution to the similar extent. 



Of course one needs the result for perturbative QCD. For the case we deal, 0{^a^) corrections to a{e'^e~ 
7* — > tt) are calculated in 

Here a may be considered as a cou pling. Asymptotic behavior of S{a) can be inferred by using renormal- 
ization group, for example. See Section 2.7.3 for detail. 

A different renormalization scheme (and renormalization scale) corresponds to a different way of summing 
up a perturbation series. 
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1.3 Overview of top-quark anomalous EDMs 

Of course top quark is (potentially) sensitive not only to the SM dynamics but also to the 
physics beyond the SM. One particularly clear signal is CP violation in top-quark sector. 
This is because the SM contribution is suppressed to many orders below the current or near- 
future experimental reachQ. Thus a non-zero expectation value of a CP-odd observable means 
contribution from new physics immediately^^. In fact there are plenty of new sources for CP 
violation once the particle content of the SM is extended. Many models, including Minimal 
SuperSymmetric SM (MSSM) or multi Higgs-doublets model in general, induce CP violating 
t-t-gauge boson interactions at one-loop level. In the CM frame of ti pair, or a particle-anti- 
particle system in general, a CP-odd observable related to them is odd under St ^ St, and 
vice versa, since CV transformation exchanges their spins s at the frame. Also in this respect 
top quark is excellent since the polarization of t is not disturbed by hadronization, and the 
information of Pt is inherited to its decay products; especially, the charged lepton is emitted 
with the angular distributiorf^ oc 1-|- \Pt\ cos 6', where 6 is the angle between the polarization of 
t and the momentum of Thus the CP violations in the top-quark sector can be measured by 
the difference — Pj between the polarizations of t and t, which can be measured statistically 
from the directions of the charged leptons. 

Since many extensions of the SM contain new sources of CP violation, we parameterize their 
effects by effective couplings in a model independent way. Among CP violating interactions 
of t, Electric Dipole Moment (EDM) interaction is the lowest dimensional operator]^. Thus in 
terms of an effective Lagrangian, CP violating effects in the top quark sector are parameterized 
by several anomalous EDMs of top quark at the first approximation. There are several of them: 
ttg-, tt'y-, ttZ-, and tbW-EDM. The first three affect the production (and rescattering) process 
of tt, while the last one the decay process. 



1.3.1 Previous studies 

Many studies^ have been done for the anomalous EDMs of top quarks. One major topic 
is how well those EDMs are measured in future colliders. Some study for e'^e" colliders with 
\/i 2mt, and the others for hadron colliders. The former is suited for the study of EW-EDMs, 
tt'j and ttZ, while the latter for Chromo-EDM, ttg. Both of these results are summarized in 
Section O. 



1.3.2 Our contributions 

Although many studies have done for e~^e~ colliders with ^/s ^ 2mt and for hadron colliders, 
there was no one for e~^e~ colliders near the threshold to date. Thus we study this case. 
We concentrate on the anomalous EDMs in the production process (tt-gauge bosons) here, 
and leave tbW-EDM for a future study. Among the three flavor-diagonal EDMs, our prime 
concern is on Chromo Electric Dipole Moment (CEDM): tig-EDM. The reason is as follows. 
As shall be shown in Section |4.6|, two Electro Weak EDMs, tt'j- and ttZ-EDM, directly modify 



This is because the CKM matrix elements for the third generation are almost diagonal: \Vti,\ ~ 1. See 
Section for more quantitative statements. 

16 Provided th at th e initial state is a CV eigenstate. 



1^ See Section 4.5.2. 
See Section |a.4| 
See Sections 4.2 and 4.4 for the references. 



1.4 Organization of the thesis 



7 



the tt production vertex, and their interactions are proportional to the relative momentum 
\p\ ~ y/mtE between ti. This means open top region {^/s ^ 2mt) is more appropriate to 
study them. On the other hand CEDM modifies the rescattering of ti, which can is suppressed 
when y/s ^ 2mt. Thus in open top region at lepton colliders, it can be studied only by real 
gluon emissions. We find that the threshold region is more sensitive to CEDM than the open 
top region. As for hadron colliders, they are sensitive to the existence of a non-zero CEDM, 
since ti pairs are copiously produced in gluon fusions. In fact our result for CEDM is that the 
sensitivity of e^e~ colliders at the threshold are worse than that of hadron colliders. However 
precise measurement of CEDM may be difficult in hadron colliders. Moreover, while detailed 
detector simulation seems to be indispensable for the serious study of the sensitivity at hadron 
colliders, it has not been done so far0. On the other hand, a full detector simulation with 
realistic experimental setup is in progress [p.25|| based on our study presented here. Thanks to 
the clean environment of e^e~ collisions, their first result for the sensitivity is consistent with 
our naive estimate. 



1.4 Organization of the thesis 

In the thesis I studyQ the physics of ti threshold in e^e~ collisions for the study of both in the 
SM and beyond the SM dynamics. 

Chapters ^ and are devoted to the SM dynamics. Total production cross section cxtot 
is studied in Chapter |^. By reorganizing perturbation series, we obtain better convergence 
of (Ttot to NNLO than ever obtained. With this result, theoretical uncertainty to measure 
rrit becomes smaller than experimental uncertainties. The chapter also includes a review of 
previous studies, of theoretical set-up such as the Green function method to calculate cxtot- A 
derivation of non-relativistic Hamiltonian to NNLO is also given. Chapter ^ deals with the 
top-quark momentum distribution. Coulomb rescattering between ti is treated to NNLO. 

On the other hand. Chapter ^ is devoted to the contributions from new physics, especially 
to new sources of CV violation. It is explained that, the CV violation in top-quark sector 
within the SM is so small that if it is discovered it would immediately imply new physics. CV 
violation in top-quark sector is parameterized by Electric Dipole Moment (EDM) interactions 
of top quark, since it is the lowest dimension operator that violates CV symmetry. Many studies 
have been done on EDMs of top quark, but they are all in the open-top region -^/s ^ 2mt. 
We analyze them in the threshold region a/s ~ 2mt for the first time. Due to the multiple 
Coulomb rescattering near the threshold, differential production cross section is sensitive to t-i- 
gluon anomalous EDM even when ti is produced in e^e~ collisions. We find that the sensitivity 
of future e^e~ Linear Colliders (LC) is worse than that at hadron colliders by factor (9(1/10) 
or less. Sensitivities for ti-^- and ttZ-EDMs are also studied, and are found to be comparable 
to those in the open-top region at LC. 

Summary and discussions are given in Chapter 5. 



At least to our knowledge. 

This work is based on the collaborations with Y. Sumino and A. Ota and with M. Jezabek and 
Y. Sumino |6|l. 



Chapter 2 

Total Production Cross Section 



2.1 Overview 



As was explained in Chapter |1|, it was considered that the mass rrit of top quark t can be 
determined most accurately by the energy scan of ti total production cross section near the 
threshold in e~^e~ collisions. However there emerge the fear that this way of determination 
suffer from a large theoretical uncertainty, due to the large NNLO (= 0{l/c^)) correction to 
the ti production cross section crtot(s), where -y/i denotes the CM energy. 

There the Green function G 



In the thesis, we follow closely the procedure used in |27 



for a ti system is rewritten in a form convenient for numerical calculations. By implementing 
"renormalon cancellation" and "log(g//x) resummation" , we improve the convergence of atot 
obtained in 



This chapter is organized as follows. In Section |2.2| , it is explained that how an non- 
relativistic (color-singlet) ti pair is described in terms of the Green function of the pair. In 
particular, the total production cross section is given by the optical theorem, which is explained 
in Section |2.2.1| . In Section |2.2.2| , order counting is explained. Also shown is that the Coulomb 
potential is a part of leading Hamiltonian. At NNLO, there arise several difficulties, which 
is given in Section |2.2.3| . In Section p.2.4| , the Hamiltonian to NNLO is shown. The Hamil- 
tonian is looked at more closely in Section |2.3| . We obtain the Green function G by solving 
the Schrodinger equation numerically. For this purpose, we express G in terms of reduced 
Green function G', which satisfies simpler Schrodinger equation than G itself. This is done in 
Section |2]^. At that stage, we are equipped to reproduce the previously obtained results for 
the total cross section, including NNLO ones. They are collected in Section |2.(j|, and exhibit 



large correction at NNLO. Hereafter, we try to improve the convergence of the perturbative ex- 
pansion. It is known that (certain) perturbative expansions involving the pole mass of a quark 
are less converging than those involving MS mass, or "short-distance" mass. It is suspected 
that this is because pole mass suffer from non-perturbative dynamics of QCD. Section 12.7] is 
devoted to this issue. There we adopt potential-subtracted mass scheme, which is proposed 
in 



59 



This improves the convergence of binding energies. In Section |2.8| , certain kind of loga- 
rithm in the Coulombic potential is summed up using renormalization group. Tentatively, the 
convergence of R ratio is best with these two prescriptions. The section contains the estimate 
of theoretical uncertainties for this case; in terms of top quark mass, it is (Amj)**^ ~ 0.1 GeV. 
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2.2 Non-relativistic description of tt system 
2.2.1 Optical Theorem 

In terms of an amplitude Ai, the unitarity [Eq. ( |A.55| )] of S-matrix is expressed by optical 
theorem, which relates the total cross section of a two-body scattering to the imaginary part 
of the forward scattering amplitude: 

o-tot(^i, h anything) = i lmM{ki, k2 h, ^2) , 

sfJi 

where /3j = 2|A;cm|/\/s, and I^cmI = l^il = 1^21 at the CM-frame of two particles with momenta 
ki and ^2- If we restrict the intermediate states appropriately on the RHS, we have the total 
cross section for ki,k2 such-and-such. For the process {pe)e~ (pe) — 7* — ti, the total 
cross section of this process can be obtained from the imaginary part of the amplitude for 
e+(pe)e"(pe) — > 7* — > ^ 7* ^ e+(pe)e"(pe)- Let {-ieQtf ■ i^^^iq) be the ti contribution 
to the vacuum polarization of a photon with momentum q^^: 

'dSe^^- (0|T[/(x)/(0)]|0) 

(0|T[/(g)/(-g)]|0) 

d^p i i 



{27:)^ {p - g/2)2 - + ie + q/2)^ - m? + ie 
X tr [(^-//2 + m)7'^(^+//2 + m)7'^] , 

where j^{x) = i{x)-yH{x), and the last expression corresponds to one- loop. By contracting 
indices, we have 

= ^ (0| T[/(g)J,(-g)] |0) = ^9,.Il^''{q) • 

Thus for unpolarized e~^e~ , 

tM{pe,Pe -^■■■^Pe,Pe) = HcQtf " ^{q'g''' - gV)n(g') X 

X ^ X] tr [u{pe)u{pe)j^,v{pe)v{pe)-fu] 
spins 

= tiAnafQlQ',U{q') . 

Here we neglect electron mass. Note that the g'^g^-part of U'^'^ do not contribute, since e~^e~- 
current is conserved: v{pe)^u{pe) = 0, where = {pe +Pe)'^- Thus, with s = g^, 

R{s) ^ -u^^^-^r^tt) ^2.1) 

4.77- r 

= I2vrlmn(g2) = _ (o| T[nq)U-q)] |0) 
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Figure 2.1: Lowest-order (=0(l/c)) contribution of gluon exchange 
between tt to the total cross section. 




Figure 2.2: Ladder approximation of vertex function F. This is suf- 
ficient for Coulomb gluon, which is instantaneous. 



where 



Ana' 
~37 



0.817pb 



2 X 175 GeV 



(2.2) 



is "point cross section", the total cross section for e'^e^ —>■ ^* ^ f^^f^^ with ^/s ^ 2m^. 

As will be shown in Section [1.6.2| , Coulomb rescatteringQ between t and i can be taken into 
account^ by modifying (one of the) ti'j vertex as [Figure |2A[| 



where 



Gi\p\,E) 



YG{\p\,E) 



Pi 

mt 



da 



'^{p\G\x' = 0) , G = 



An " ' ' " H -{E + iTt) ' 

and if is a Non-Relativistic Hamiltonian. Thus to the lowest order of 1/c-expansion 

d^j9 i i 



iW^^q) = {-l)Nc 



27r 4 /E ^ Tt 
2 ^ 2 



— — + 1 — 

2mt \2 ^ 2 



X 



PL 

2mt 



X tr 



1+7^ 



1-7° 



reseat " 



-r 



t2{g^''-g^'g''')Nc 



d^p 



G{\plE) 



(27r)3 

i2{g^'' -g^'g^')NcG{x = Q,E) 



or 



n(g^ 



2N, 



c. 



G{x = 0,E) , Im -z (0| T[r{q)U-q)] |0) = 6NcG{x = 0, E) 



^ We have in mind the rescattering due to the exchange of Coulomb gluons. 

^ Only the effect of "soft" gluons can be taken into account by the Green function G. See Section 2.2.5 
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where G{x, E) = {x \ G \x' = 0). Thus up to short distance corrections and relativistic correc- 
tions, R ratio near the threshold is given as follows: 

R{s) = -NcQl -\mG{x = 0, . (2.3) 

/ s m'^c 

Note that the velocity (3 oi t ai ti CM frame is 



o,, l!!i = i_^^ = J_. (2.4) 

V s s 7^ 

If one expands as amplitude Ai into eigenstates of orbital angular momentum L, Ai = 
^^A^^^), then lA^'^'^^l oc P^. Thus S wave {L = 0) is the leading contribution near the 
threshold. We sometimes use S-wave projected notation of Green function G{r, r') defined in 
Section |A.2.4: 



G{r = 0, r' = 0) = G{x = 0,E) . (2.5) 
Here the energy dependence is implicit. 



2.2.2 Order counting 

The non-relativistic expression of a quantity is obtained by expanding it with respect to 1/c. 
Noting that 

^ Airhc \c J ' 

it is useful to introduce a dimensionful (but order unity) coupling = a^^{fis)c for order 
counting; here fig is a "soft scale" , which means the typical scale for the ti potential. Thus the 
QCD correction as is the same order to the relativistic correction /3, where /? = \/l — 4mf / s 
is the velocity of t and t at their CM-frame: 

n ^ / N asif^s) as , . , . 

P = - , asii^s) = — , — ^ = — = C I . 2.6 

C C jj V 

In fact as we shall see near Eq. ( p.l2| ), Pt ~ \jTt/nit ^ 0.1 ~ near the threshold. One can 
clearly see that one need to sum up as{fis)/P to all order. The Leading Order (LO) is C(l), 
the Next-to-Leading Order (NLO) is 0(l/c), and the Next-to-Next-to-Leading Order (NNLO) 
is 0(l/c2). 

The following relations may be useful to count the order of 1/c: 

dim[mrc] = dim[rp] = dim[/i] = dim[l] , dim[l/r] = dim[p] . 

Since dim[as] = dim[c] or dim[as/r] = dim[i?], the Coulomb potential is the same order to the 
kinetic term. Thus it is a part of the leading-order Hamiltonian Hq. 



2.2 Non-relativistic description of tt system 
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2.2.3 Complications at 0{l/c^) 



There are several complications at (9(l/c^) and beyond. Some are already in freely-propagating 
ti pairs, some are in the ti potential V, and some are especially in the Non-Abelian part ^na 
of V. 

For the free ti propagation, (i-1) the Hamiltonian contains p'^ term besides term, which is 
the standard kinetic term, (i-2) The wave-function normalization of the production ti current 
jtt depends on the momentum transfer g^, or in other words, anomalous dimension of ju is 
non-zero, (i-3) At 0{l/c) and below, the decay width Ft of t quark can be incorporated by 
replacing H ^ H - iV ||,|T§. At 0{l/c^) and higher, this is no longer justified. These three 
issues are discussed just below. 

For the ti potential V , there are log(g^/ n1) corrections p3|-|45[ 
to 0{l/c) aside from a finite term. However to 0{l/(?) 



This is the only correction 
besides the log corrections, (ii-1) 



there are also corrections called Breit-Fermi potential Vbf • These corrections are common 
to the potential due to both Abelian (7) and non-Abelian gauge {g) boson exchange, (ii-2) 



There is also a correction [52| especially due to non-Abelian nature of a gluon, which is absent 



at 0{l/c) and below. These two issues are discussed in Section 2.3 



Current normalization , C2 

For i?-ratio, the lowest in both f3 and as expansion is 0{(3^, a^). There are several corrections 
beyond the lowest order. Only relativistic corrections are considered hereQ For = and to 
all-order in /?, the production cross section is 



(2.7) 



where j3 is velocity of t at tt CM frame [Eq. (|2.4| )]. The Green function formalism [Eq. 
should reproduce this result. The Non- Relativistic Hamiltonian Hj^^^f"^ for free ti is 



(NR) 
free 



4ei - 2mc2 



PI 

m 



+ 



where H^^l is defined in Eq. ( |2.17| ). A form factor depends on the momentum-transfer in 
general, and so does the wave-function normalization for ti vector current j^{q)- A Taylor 
expansion, or derivative expansion in coordinate space, of R{s) reads 



R{s) = Q',-lm[-z (0| T[r{q)U-<l)] |0) 

2 s 'm?c 1 2m?c^ 



where C| 



(cur) 



2 -niiu -an J ^(cur) , ^(cur)'^r "I" 

0{as) and C^'"' 



Im G(r, r') 



r,r'^ro 



const. + (as); the constant of 6*2™'^'' is determined below 



(2i 



in this section. The reason for introducing the cutoff parameter ro(7^ 0) shall be explained also 



^ Others are considered in Section 



2.5 
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below in this section. If the potential between tt is neglected, a part of the expression above 
can be rewritten as follows: 



Im 



^(cur) ^(cur) ~l~ ■^^r' 



m 



+ 



\r') 



Im 



1 



\r') 



+ Im (r| 



P 



m 

4 



1 



m 



m 



(2.9) 



Here we used 

{/\r + /^r'){r\G\r') = -{r\{p^G + Gp'')\r') . (2.10) 
First, consider the case = 0. With the explicit analytic calculations^ of the Green functions 



in Section B.4.1, the /3-dependence of R in Green function method is 



while relativistic QFT says [Eq. (|3)] it is (3{l - (3"^ /?>). Note that u = p + 0{P^) is defined in 



Section B.4. Thus we obtain 



1 
3 



(2.11) 



As shall be seen in Section p.2.2| , soft gluon corrections are suppressed by as{fis) = C'(l/c), the 
lowest matching coefficient is sufficient for the 0{l/c^) calculations, since c!^^^^ is a coefficient 
of (9(l/c^) correction. 



Treatment of the finite width F^ 

Concentrating on freely-propagating ti, the Green function^ at the origin is 

1^ 



mr 



^^(0,0) = ^Re 
An 



rrit \c 



(2.12) 



Note that \jElmt = P + 0{p^). Thus effectively p ~ \/T\/rnt ~ 0.1 near the threshold. By 
using 



VE + iTt 



,/WTf^ + E 



1/2 



+ i 



^/WTfl-E 



1/2 



X sgn[Ft] , 



The relevants are ImGo and ImG2. 



^ This is derived in Section B.4.1. The relevant one is ImGn. 
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in general, we have 



+ ■ (2-13) 

It is known tliat tlie effect of finite decay width Tf can be taken into account by the 
replacement E ^ E + iTt to LO 0] and NLO [|T§]. However to NNLO with non-zero i t, 
imaginary part of the Green function is no longer finite when one takes the limit r = r' ^ 0, 
which is required by the optical theorem: 

= (-cr- 4) 1.1) , 

for r > r'. This means that the effect of the finite width Tt is not properly treated by the 
replacement E ^ E + iTt at C(l/c^). Here following ||2^, we regulate this singularity by 
putting r = r' — ^ ro 7^ 0, as shown in Eq. (|2.8|) . 



^/W+fl + E 



2mt 



2.2.4 Non-relativistic NNLO Hamiltonian H 

The NNLO HamiltonianQ H, which is relevant to us, is 

p2 Cptts 



H = Ho + V,ir) + Uip,r) , Ho 



[2.U) 



where Hq is the leading-order Hamiltonian, Vi{r) is radiative corrections to the Coulomb 
potential^. 



VAr) 



CpCls 



Atxc 



{2/?oln(/iV) +ai} 



+ (^)' {/3o' (41n2(/iV) + + m + 2/5oai) ln(/iV) +02} 
and f/(p, r) is a momentum-dependent potential, which includes only 0{l/(?) terms: 



P 



^ CaCfo^I 



CpO-s 



r 



1 + 5(3)(r) 



+ 



SCpds L S Cpag 1 



:s-r) 



2mlc^ 2mlc^ 



Amfc^ 2mtr'^c^ 2m^c^ 



P^> + 



3m? 



:2.i5) 



^ When we say simply "Hamiltonian" H, it means "Non-Relativistic Hamiltonian" H'-^^\ 
^ We may sometimes use the word "Coulombic potential" Vq for the Coulomb potential with radiative 
corrections: 



Coefficients 02 etc. are defined in Section 2.3.4 
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where { , } means anticommutator, and the last expression is for the S-wave state, L = 
0, 5^ = 2,r*r-? = 5*^^/3. The first term in U{p, r) is the relativistic correction to the free 
Hamiltonian, and the second term VnaI'") is due to Non-Abelian nature of a gluon. The rest 
is Breit-Fermi potential Vbf(P7 f)- In the next section, we shall see each of them more closely. 

We need to calculate the Green function G{r,r') for the Hamiltonian H. Because of the 
logarithmic corrections to Coulomb potential, G{r, r') cannot be obtained analytically^. For 
numerical calculations, there are two obvious complications for the H above. One is p'^, which 
makes the Schrodinger equation a fourth-rank differential equation, and the other is the 5- 
function. These two difficulties can be circumvented by rewriting H appropriately. We do this 
in Section 



2.2.5 Soft corrections and hard corrections 

It is known that the matrix element is singular near the threshold due to "soft" gluons, 

if one calculate to a fixed order of a^. This is called "threshold singularity". Here the "soft" 
gluon means that its momentum is ~ asirit. However it is also known that this singularity 
is absent once the terms of are summed over, and this summation is achieved in 

Non-Relativistic Green function G with Coulomb potential^ V . Thus Green function method 
[Eq. (p.8|) ] is convenient for the threshold physics. 

On the other hand, if the gluon, exchanged between tt, is "hard", which means its momen- 
tum is ~ m^, it is not treated by the non- relativistic potential, but by the matching coefficients 
of ti production current. It is called "hard vertex correction" which is the same as usual 
radiative correction to the vector vertex: 

r ^ t'^ (l - 2Cf^^ + 0{al)^ (2.16) 

for each vertex. Here fih — mt is a "hard scale" which is a typical scale for short distance 
physics. We limit ourselves to the contributions from the "soft" gluons for the moment. The 
matching coefficients, which represent the contributions from the "hard" gluons, are determined 



in Section 2.5 



2.3 Non-relativistic tt potential V to NNLO 
2.3.1 Normalization of wave functions 

In QFT, fields of spin-0 are dim-1, while those of spin-1/2 are dim-3/2. Thus if one uses the 
convention]^ that the dimensions of the creation-annihilation operators are the same for spin-0 
fields and spin-1/2 fields, the dimension of c-number wave functions for spin-1/2 fields is higher 
by 1/2 than that for spin-0. Since the dimensions of states with different spins are the same|^ 



Actually, the Green function without Vi can be obtained analytically. See Section 2.4.2. 
^ Note that V represents gluons of its momentum is ~ asirit. The relevant scale for these soft gluons (or for 
the potential) is denoted by mus- 

This is the convention that is usually used in (modern) relativistic calculations. 

Since {ip\ip) is probability, \^JJ) is dimensionless. While the dimensions of bases are determined by orth- 
normality relations. For momentum, (p\p') = {2t:)^S^'^\p — p') means dim[|p)] = dim[p^'^/^]. For spins, |s) is 
dimensionless in usual convention. 



2.3 Non-relativistic tt potential V to NNLO 
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in NR-Quantum Mechanics, one has to rescale the wave functions for spin-1/2 fields by hand. 
We may proceed as follows. Free Relativistic Hamiltonian i^fj-ee P^^^ 

H^^l = v/(mc2)2 + (pc)2 X 2 



(2.17) 



.2 , P P 



2mc^ + — 



4 / 1 

+ - 



Thus, the propagator for one-particle state in a non-Lorentz-covariant QM is 

i i 



2m 8m^c^ 



On the other hand, in Lorentz-covariant QFT, it is 



1 



p2 — + ie 2p° 



1 



1 



p2 _|_ ^2 _|_ pO _[_ /pi! _|_ _ 



or 



+ m 



+ m) i 

2_ 



i{fi — m) i 
-m? + ie pO - + ^2 _^ 2p0 



Thus the non-covariant normalization of spinors, which is used for non-relativistic calculations 
hereafter, is 



u{p, s)u{p, s) = , "^vip, s)v{p 

zp 



— m 
2p0 



or in Dirac representation, 

s) = 



p^ + m 



( 



2pO 



X 
pa 



\ 



-X 



"^{Pi ■s) = u'^{p-, s) = Cu{p, 



\p^ + m / 

where C is charge conjugation matrix defined in Section [A. 3. 2 



2.3.2 Breit-Fermi potential Vbf 

The potential V (in momentum space) between tt can be obtained by calculating the matrix 
element Ai for the elastic scattering of tt, t{p)i{p) t{p')i{p'), in QFT and comparing 
this with the formula in Born approximation, which is proportional to V{k), where k is the 
momentum transfer. Working in the CM frame of tt, the system can be dealt with one particle 
QM [Section [XXTj . Let us define^ k = p' - p, K = {p' + p)/2, oi p' = K + \k, p = K - \k. 

Note that k is not relative momentum but momentum transfer; likewise K is not total momentum. 
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Using p° = mc^ + p"^/ (2m) + (9(l/c^) etc., 

r 1 1.2 

m(p')7°m(p) = Xt 



1 -k^ + 2icrf(kxK) 1 



1 2K' - ie'^^y 



+ 



2m ' V C'^ 



The corresponding expressions for 

v{p)-f^v{p') = -v^{p'){YT'v''{p) = +u{p')-f^'u{p) 

— k, K — > —If, (Tt (Tt, xi'^ ~^ X'p ■ Thus in Coulomb 



can be obtained by replacing k 
gauge0 [Section |A.6.2|| , 



iM = i-CF){-tgsf u{p'h\{p) ■ u{p')j\{p) ■ + u{p')Yu{p) ■ u{p')-i^u{p) ■ D'^ 
where Vbf is the Breit-Fermi potential |[41|| : 

VBFip',p) = {p'\VBF\p) 



p'\p^ 



(p'-kf , ik-p] 



y 2m?c^ \ k^ fcV 2m2c2 \ {k^ {k^ 



+ 



3i 



kxp' ^ ^ kxp ^ 



k' 



k' 



2m? 



:s-k) 



where k = p' — p, S = St + St and St = cr(/2, St = (7't/2. There are several remarks. 
Since the overall phase of M. is unphysical, relation between iM. and V is determined so as 
to reproduce the Coulomb potential at the leading order. We rewrote K with p, p', fc in a 
symmetric way: 



AK'^ = 2p'^ + 2p2 _ fc2 



2\2 



4 [k-KY = 2 {p'-ky + [k-pY - {k') 

kxK = ^{—p'xk + kxp) . 

This makes the hermiticity of Vbf transparent. Due to CV invariance of QCD, Vbf is symmetric 
under St ^ St [Equations and [A.52| |. They are conveniently expressed in terms of the 
total spin S: 



SfS-t = -(s'-- 
2 V 2 



[Sfk){St-k) = ;:{{S-kY-^ 



The Vbf depends not only on k but also on K. This means the potential is momentum 
dependent [Section |A.2.3|] . This can be seen clearly in the coordinate-space representation 
given below. 

In fact Vbf depends on the choice of gauge, and is not hcrniitian in general. However Coulomb gauge 
gives hermitian result. 



2.3 Non-relativistic tt potential V to NNLO 
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The Breit-Fermi potential Vbf{i^) in coordinate-space can be obtained by Fourier transform: 

where r (r') is the relative coordinate for the initial (final) state. With 

d^p' d^p p' r' — p r 

= d^kd^p = k-r' + p-{r' - r) 

= d^p'd^k = p'-{r' - r) + k-r 

= d^Kd^k, = K-{r' + r)/2 + kir' -r) , 



and 



d^k 1^ k'y 1 k' 1 k'y 



(27r)3 [ ' fc2 ' |fc|' (fc2)2 

. .m/ N 1 f.ii nr'r^\ 477 ..,3^, ,21 zr^ 1 1 / .. rVA 



we have 



r 2m?c^ \ r ' J 2m^c^ r'^ m^c^ \ 3 y 



where { , } means anti-commutator, and p and L = rxp are differential operators that 
correspond to linear and orbital-angular momentum, respectively. Note that L and r commutes 
[Section |A.2.2|| . In this form, one can clearly see that Vbf(^") is hermitian. The following 



relations may be useful to compare with the formulas in other literatures: 

^ ^ (Srr)(Spr) 1 / JS-r) 



I- , p4 + ^ + 47r5(3)(r) = — (p^ + -r-{r-p) p^ 
[r ) r-^ r \ r y 



2.3.3 Non-Abelian effect V^a 



There is also a potential |Q proportional to Ca, which is zero for Abelian gauge group (like 
QED) [Figure IJ: 

One can see that the potential in momentum space is non-analytic, and in coordinate space it 
is a strong attractive potential near the origin. 
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I 
I 
I 

1 



Figure 2.3: Feynman diagram that contributes to Fna- 



2.3.4 Radiative corrections Vi for Coulomb potential 

There are also radiative corrections to Coulomb potential. We call this potential Coulombic 
potential Vq: 



^2.20) 



where the coupling in V^-scheme is [^^-|45[| 



+ 



1 + 



An 



/5o In — + ai 



4:71 



k 



(2.21) 



2\ 2 



k' 



/3o\n—] - (2/3oai + In — + as + ■ 



a'f'ik'] 



1 + 



a 



MS ,'^2^ 



47r 



An 



a2 + ■ 



for /i^ = fc^ 



Here the first expression is for fixed order series, and the second for Renormalization-Group 
(RG) improved series, jj? = k"^. In Section |2.2.4| we wrote the log corrections as Vi. 



Vc{k') = -C,^^ + V,{k'). 

The potential Vc(r) for coordinate space can be obtained analytically only for the fixed- 
order series: 



Vrir) = -Cf 



MS 



(2.22) 



where 

«r(i/r;/i) = «r(/i^) 



1 + \ {2/^0 ln(/ir e^) + ai} 
An 



(2.23) 



+ 



We can see 



An 



a^{l/r; /i) ~ aj(fc^ = /i^) , r' = r e' 



n 



/?o' 41n'(/ire^) + — + 2(/5i + 2/5o«i) ln(/ir e^^) + as 



,y^/.2 



'2. ,,2n 
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aside from 7r^/3. Here 7e = 0.5772 ■ ■ ■ is Euler-Mascheroni constant, and originates from the 
Fourier transform: 

d^A; ln(fc2)e^^-^ "1 r, / 2^ o i 

Coefficients oi and 02 are coefficients for 0{l/c) and (9(l/c^) corrections, respectively: 

31^ 20 31 10 

ai = — TpUf = Uf = 4.77778 ■ ■ • , 

9 9-^39^ 

/4343 ^ 2 22 \ ^2 A798 56 \ ^ ^ 

= 456.749 - 66.3542 nj+ 1.23457 = 155.842 ■■ • . 

Group theoretic factors Cp etc. are given in Section |A.1.1| . In the first calculation of 0{a1) 
correction [^], the 47r^ in C\ term of 02 was incorrectly calculated to be 67r^. With this, 

af'*^) = 634.402 - 66.3542 + 1.23457 = 333.495 ■ • ■ . 

Thus the correct value |^ of 02 is less than half of a!'2^'^\ 

Since 02 is a coefficient of 0{l/(?) correction, one can expect better convergence with 
the correction to 02. However for fixed order calculation, \\i{k/ ^g) term disturbs this naive 
expectation. While for RG improved calculation, the size of 0{l/c^) correction is diminished 
by 2. 

Coefficients Pq and [3i determine the RG flow of QCD coupling ^^(/i). These flrst two 
coefficients are independent to renormalization scheme: 



dlnjj? 47r \ At: 

das _ _Po_ 2 _ 3 _ 4 
^di2~ 27r"^ 87r2"^ 3271^"^ 



where 



Po = - ^Tprif = 11- \s = 7.66667- • • , (2.25) 

34 20 / 19 \ 

A = yCl - -CATprif - ACpTpHf = 2 ( 51 - -rifj = 38.6667- ■ ■ . 

At the leading order of running, one have 

, 2, «s(M2) Att 

as(/i ) = - 



where = exp (^ p^a'^M^) ) ■ identify ^^(/i)/?! and (as(/i)/7r) ln(/i^/A;^) to be the 

expansion parameters, then the coefficients are 

^ = 1.91667--- , — = 2.41667--- , 

4 16 
ai ao 
— = 1.19444 ■ ■ ■ , — = 9.74014 ■ ■ ■ . 

4 ' 16 

One can see that the finite term of 0{af) correction 02/16 is large compared to the other. 
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RG improved Coulombic potential 

As can be seen in Eqs. ( p.21| ) and (|2.23|) , leading log's, ^(agln)", in Vc can be summed up 
either in momentum space (/x = k) or in coordinate space (/i = 1/r'). However as we saw 
above, there is an extra term vr^/S in term for the coordinate space potential. Thus it 
seems that RG improvement in momentum space is better. The potential in coordinate space 
is obtained by numerical Fourier transform. There is another advantage^ for this prescription, 
which is related to the "renormalon cancellation", where the infrared part of Vc{k) plays an 
important role. However let us postpone this subject until Section p?7| . 

Scale dependence of av{q,fis) and a^^{fis) are shown in Figure for fig = 20 GeV, 
75 GeV, and q. One can see that perturbative convergence is improved by the RG prescription 
Us = Q over the whole relevant momentum scale. Since the coupling becomes stronger as higher 
orders are taken into account, it is expected that the binding energy becomes larger with higher 
order corrections. Scaling behavior of MS coupling is also shown for comparison. It is similar 
but higher order corrections are much smaller. In fact one can hardly distinguish 2-loop and 
3-loop running. This is because 02 (and Oi) is large. See the relation of V-scheme coupling to 
MS-scheme coupling [Eq. (|2.21|) ]. Also shown is the straight line q/{Cpmt/2). The solution 



q/{CFmt/2) = aviq), or 

fis = CFav{fis)m/2 

gives Bohr momentum pb = Us- We can see that — 20 GeV and ay(pB) — 0.16. 



As shall be explained in Section p.6.1| , the normalization of the total cross section is de- 
termined by the strength of the attractive force (W/dr between ti. This is shown in Fig- 
ure |2.5| . One can see that the perturbative convergence is worse for /i^ = 20 GeV ^ than 
for /is = 75 GeV or = q. 



2.4 Reduced Green function G' 



2.4.1 G in terms of G' 



As was explained above, the Hamiltonian H to NNLO (= 0{l/c^)) [Eq. (|2.14|) ] is not convenient 
for numerical calculation. Thus following |]2^, we rewrite the Green function G for H in this 
section. 

We are interested in a Green function G{r-^r') that is projected to the ^-wave. In operator 
language, it is defined as follows [Section |A.2.4|| : 



G 



H-uj 



E + iT 



t , 



(2.26) 



G(r,r') = {r\G\r') , G(r,r') ^ J ^^G(r, r') . 
Here H is defined in Eq. ( p. 141) . Using operator identities derived in Section |A.2.2 



[p^ ipr 



47r5(3)(r) + 



2L^ 



r 



It was argued in |(46| that large theoretical uncertainty remains even after the RG improvement of Vc ■ This 
claim was based on a large discrepancy between results of renormalization-group improvements in momentum 
space and in coordinate space. However in view of "renormalon cancellation" , momentum space is better. 
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Figure 2.4: The momentum-space couplings ay vs. momentum trans- 
fer q at LO (dot-dashed), NLO (dashed), and NNLO (sohd). Two are 
fixed order (^^ = 20 GeV, 75 GeV), and one is RG improved {^g = ?)• 
Intersection of ay and the straight line q/ {Cprnt/'l), which is also 
shown, gives Bohr momentum p-Q and the coupling ay (pb) at the scale 
Pb [Section |B.3(|. Running of MS coupling is also shown for reference. 
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or 

r2 



5^''\r)^[Ho,ipr] ^ + 



the momentum dependent 0{l/c^) potential U{p, r) can be rewritten as follows: 



^ \Ho,-) + — — [Hq , ipr 

AuitC^ { r ) 12 mtc^ 



where the last expression is for the 5'-wave. We concentrate on the ^-wave hereafter. To 
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summarize, 



U = hHl + {//o , Oe} + [Ho ,Oo\+W , (2.27) 
^ = , Oe = —. 5- - , Oo = -tt: T iPr , 



Arritc^ AmtC^ r 12mtC^ 

Clal ^2 , Ca\ I K 



Note that Oq is anti-hermitian, and thus U is hermitian. Using this expression for U, we can 
rewrite G up to 0{l/(^): 

G 



Ho + Vi + U-cj 

Hq — UJ Hq — UJ Hq — UJ Hq — UJ Hq — 00 Hq — u 

-b+il-boO-OE- Oo)-rT^{l -bu-OE + Oo) 

ijQ — U) n.Q — OJ iIq — LO iiQ — UJ ijQ — OJ 

-b+{l-bu;-OE-Oo) „ ^, (l-bu-OE + Oo) , 

H0 + V1 + bcu'^ + 2ujOe + W — CO 



or 



(2.28) 



G^G' -{bu; + OE,G'}-[Oo,G']-b + o(^-^^ , 
G' = -^—^ , H' = Ho + Vi + 2uOe + W , S = lu - buj^ . 
Note that H' does not contain neither nor S{r): 

H'^P--^ + V,{r)-^^, (2.29) 



where 



a, = a J 1 + — , E = {E^iY)[\- 



and 



2 rxK? ) \ \m(9- 



_Cl(2 , Ca\ fa.^^ 



2.5926 = 0.06637 x 



V0.16/ 

Thus numerical calculation of Green function is much easier for G' than G itself. In terms of 
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c-number functions, they are related as 
G{r,r') ■ 



E + iT 3Cf /o n 2 
1 + ^^ — ^ + 



2mc2 ' 4 \ c / V ma.^r ma.,r 
IIC^ ^a,\2 /Id Id 

c 



12 



-r + 



ma,r dr ma^r' dr' 



G'{r, r') 



Amc^ Aurr' ^ ^ ^ ) 



2mc2 



ma,r ma^r 



12 



G"(r, r') 



+ 



4mc2 Anrr 



—6(r - r') + o(r,r'\ 



(2.31) 



where we used the relation Eq. (|B.9|) for the Coulomb Green function. Note that the difference 
is higher order. Thus, 



ImG(ro,ro) 



1+ (- 

c 



-C 



F _|_ lie*!, 



SmasTQ 12 



Im 



(2.32) 



2.4.2 Green function for Coulomb plus potentials 

In this section, we consider the Hamiltonian, 



m 



K 

9 9 ' 



^2.33) 



because (i) the Green function G (projected to S-wave) can be obtained analytically; and (ii) 
the Green function G' for H' in Eq. ( p.29|) without log corrections Vi{r) can be obtained by 
replacing as ^ as, E S. The results for pure Coulomb potential, which we shall use several 
times, are also obtained by expanding in k. 
The Schrodinger equation for G{f, r') is 



E 



m 



K 



r mr^ 

where r = |r|. We are interested in the Green function G{r,r') projected to S'-wave: 



E 



-19^ 
mr 5r2 



K 



mr^ 



G(r, r') 



-1 



Airrr 



-6{r — r') . 



To simplify the equation, we introduce the following notations: 



G(r, r') 



g{r, r') 



z = CpasfTir 



V 



AE 



(2.34) 



Here 2; is r measured in unit of (twice) the Bohr radius pb, and l/i/^ is in the unit of the 
Bohr energy Eb [Eq. (|B.7|) ]. Thus by dividing by {CpasYm, we have 



^ + + ^ + ] 9{z z') 



-1 



4:7iCF(y. 



-5{z - z') 
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where g{z,z') is just g{r,r') witliout rescaling: 

g{z,z') = g[r = z/{CFasm),r' = z' /{CFOism)) , 

wliile tlie normalization of (5-function changes according to tt6{z) = Im [l/{z — ie)]. Solutions 
with (5-function source can be obtained by the solutions without (5-functions: 

''^ = ATiclas W i9<i^)9>{mz' - ^) + 9<{z')g^iz)e{z - z')] , (2.35) 



where is the WronskianP^ 



W^W{g^,g^;z) 



9>{z) g<{z) 



9>{z) g'^{z) - g<{z) g'^{z) 



and the function gy{z) {g^{z)) is the solution to the homogeneous differential equation: 

which is "well-behaved" at z — oo {z ^ 0). Since we use the value of the Green function near 
the origin, it is important to know its behavior around there. For k, > 0, which is the case, the 
potential is dominated by l/z"^ near the origin. Assuming a solution of power behavior, 

\dz^ ^ z^J - ' 



one obtains 



d = ^ = d± , (i+ + (i_ = 1 , dj^d_ = k . (2.37) 



This may indicate that the wave function is collapsed into the origin when the attractive force 
is as strong as k > 1/4. However we don't have to worry about this because k ~ 0.06 for 
our case. For small k, (i+ ~ 1 — k, (i_ ~ n. Thus in order to obtain the results for Coulomb 
potential, it is convenient to write d^{= 1 — d_) in terms of d^. For both choice of d±, g{z)/z 
diverges at the origin, but it's milder for d^. Thus g<:{z) ~ z'^+ near the origin. For Coulomb 
(k = 0), two independent solutions are 

z -^z^ + 0{z^) , l-zlnz + 0{z'^lnz) . 

In order for g{z)/z to be finite at the origin, g<{z) z — 

Further calculations are done in Section |B.2| . Especially, results for pure Coulomb potential, 
which we frequently refer to, is collected in Section p.3| . 

One can easily show that the Wronskian W for the two solutions of a 2nd-rank homogeneous ordinary 
differential equation without Ist-derivative term, is independent of z: 

dW 



Az ' 



Thus one can calculate W wherever he/she wants. 
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2.4.3 Explicit expression for G' 



The optical theorem shows that all we have to do is to calculate the Green function G, or G", 
near the origin. The reduced Green function G' is defined in Eqs. ( |2.28| ) and ( |2.29| ). Thus we 
numerically solve the Schrodinger equation 



m r 



9{r) 



0, 



(2.38) 



or, equivalently, 



z) = 



where 



2 {CFasfm ^, 



AS 



H z = fire' 



and 



^0= 1 + 



2m 



There are two independent solutions for this differential equation. Assuming their form near 
the origin as 



g±{z) = z"^ [1 + z {Lt \^\fi'z) + Lf \n{fi'z) + L±) + 0[z'')] 



(2.39) 



we have 



± 



LI 



2d± ' 
2d± + 1 
2dl ' 



2d. 



Ml + 2d± + 1 ^ 2d± + l„ 1 



44 



+ 



With these two solutions, solutions (nearly) regular for 2; ^ 00 and z ^ 0, respectively, are 



9>{r) = g-{z) + Bg+{z) , ^<(r) = g+{z) 



2.5 Matching of NRQCD with QCD 
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With the normahzation in Eq. ( |2.39| ), the Wronskian W with respect to z is the same as 
Eq. (^): W = VI - 4k. Thus we have 



Att Att 
G'{r, r') 



1 9>{z)g<{z') 



ni^cA-nCFCisW rr' 



[z > z' 



VI - Ak [ zz' 

Cpds 



g+{z)g^{z')_ ^ ^ g+{z)g+{z') 



zz 



+ Bz-^"- [l + 2z \xi\fi'z) + ln(/i'^) + L^^^^ + ■ ■ ■ } 



+ (vanish when z 



0). 



5z 



-2d_ 



i=± 



(2.40) 



where — >■ means r = r' = ro — >■ 0. 



2.5 Matching of NRQCD with QCD 



Having taken into account the "soft" gluon contributions, next we consider the effect of "hard" 
gluons. It is determined by matching the NRQCD calculation to QCD one. Short distance 
coefficients Cf^^\ 02^"^^ for ti production current, defined in Eq. (|2.8|), are determined in the 
same section to the lowest order of a^. Higher orders can be determined by matching the result 
for Green function method with that for usual perturbative QCD. The Green function method, 
or NRQCD, can be applied when /5 -C 1, while pQCD can be applied when as/ (3 ^ 1. Thus 
both formalism is valid when -C /5 -C 1. For this energy region. Green functions can be 
calculated perturbatively. While relativistic pQCD calculation was done in |H0|. Both of these 



results are summarized in Section B.4 



and C, 



(cur) 



Actually for cf^ 



By demanding these two results to coincide, we obtain 
we need only to the leading order, which was obtained 



already in free-propagating limit [Eq. ( |2.11| )]. With 02^^"* = 1/3 



^(cur) ^ ^(cur) + K' I ^^^^ 



2m? c 



CI 



(cur) 



CI 



(cur) 



2 

c. 



2 / -2Cf 



Im 



1 - 



11 
— ( 
12 



r,r'— >ro 

E + iT 



Im 



G{ro,ro 
E + iT 



1 + 



G"(ro,ro) 



(2.4i: 



where the equalities are up to 0{l/c ). Here we used Eqs. (|2.1CI| ) and ( |2.26|) , or 



^^-^-{E + ^T) + o('- 
m r \c 



G{r,r') 



A. 



m 



ImG(r, r') 



Im 



CEO'S 



Anrr 



-6{r — r') , or 



+ E + iT) G(r, r' 



+ 
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The same relation holds also for G'{r,r') since their difference is (9(l/c^). With the definition 



(cur) _j_ / 



IT 



(smja^ro ' 12 



we have 27 



'-Cf + Ca In 



+ 2/3oln 



19.17 + O.eilllln. + 0.251199 uh 



+ 38.3818 In 



(ref) 




4 



+ ( 22 - -nH In I — 



mt 



where 



a 



,(ref) 



39 



18.3077 + 38.3818 In 
e2-TE 2.07433 



+ 15.3333 In ( — 



2m+c 



rritc 



— -Cs + ttM -ln2- — 
S3 ^ I o 



-1.5215 



2 



151 13 



179 8 



-5.71378- ■• 



^ 9 



0.502398- •• . 



(2.42) 



(2.43) 



(2.44) 
(2.45) 



36 2 

= 1.22222--- 

44 4^ 

y " 9' 

Note that the velocity dependence of the matching coefficient is determined by single number: 
= 1/3. Thus it is highly non-trivial that the matching is surely possible. 
Let us consider the ro dependence of R ratio for the case = 0. Since only io and Lq , 
among £i and Lf, have imaginary part, we see from Eq. ( p.40|) , 



rn?c 



Im G'(ro, ro) = {im [B] z'^'- + Im [L+ + L,\ | 

vl — 4k > 



v/r^4^ 

Here we omit the terms that vanish when zq 0. Since k is a coefficient of (9(l/c^) interaction, 
we need only the terms with the first (or less) order in k. Thus 



ImC 



-2d- 



-2d_ In 20 



2rf_ In(zo) + 



1 - 2fi; ln(;zo) - 2k^ ln(zo)(l - ln(;zo)) + ■ 



2.6 Previously obtained results for R ratio 
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where Zq = CF<ys{f^s)''TT'tro- While vq dependence of matching coefficients [Eq. ( |2.43| )] is 

C ~ 1 + 2Kln(ro/rJ'''^^) . 

Thus we can see that In tq singularity is cancelled in R ratio to the order we are concerning 
However for finite width Ft, the expression for R diverges with l/rg and In(ro). This is due to 
the incomplete treatment of the width Fj. Here we follow the treatment of [^]; that is, expand 
the expression with respect to tq as detail as possible, and drop the terms that vanish when 
To ^ 0. We choose tq = r^'^^^ 

2.6 Previously obtained results for R ratio 

In this section, we show R ratio for freely-propagating tt with and without finite decay width 
Ft, Leading Order (LO) Coulomb rescattering with and without F^, Next-to-Leading Order 
(NLO = 0{l/c)) Coulomb rescattering with Ff, and Next-to- Next-to-Leading Order (NNLO 
= (9(l/c^)) Coulomb rescattering with F^. All of these are for fixed order asifJ^s) calculations]^, 
and are obtained before our work, except those with a|'^™. 

2.6.1 Analytic results 

Analytic formulas for R ratios [Eq. ( |2.1| )] are available for the following cases: 
3 / E 

R = -NcQf\ — , for free; without F^ , 

2 M nit 

R = ^NcQh ^^' + ^' + ^ , for free; with F, , 

3 2 T^Cpas ^ fmt 



R = ^NcQi ^_^J^ , ZE = TTCpasJ^ , for LO; without F^ , (2.46) 



where E = ^/s — 2mt. The second one is given in Eq. ( |2.13| ), and the last one is in Eq. ( |B.8| ). 
These results are shown in Figure p.6| . We can see several points. Because of non-zero width Tt, 



the production cross section is non-zero also for E < 0. More prominently, attractive Coulomb 
interaction between ti enhance ti production cross section by many times. This is because the 
ti wave function |?/^(0)p at the origin is enhanced by the attractive potential. 

Wave function ^(0) at the origin and potential V{r) 

The magnitude of a wave function at the origin is related to (the expectation value of) the 



attractive force d\^(r)/dr |^^. We are interested only in S wave. Higher orbital angular 

; to lp 



momentum states may not contribute to '?/'(0). Thus we define 

1 yij-) 



An r 



The K in G' is evaluated at the scale fig, while the k in C is evaluated at fi^- However, the difference 
between as{^J,s)'^ and as{iihY is higher order, at least formally. 



17 



LO and NLO corrections with RG improvement were also calculated; these are summarized in Figure 2.10 
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Figure 2.6: R ratios for the situations that analytic formulas are 
available. We used mt = 175 GeV, Tt = 1.43 GeV, as{^is) = 0.16, 
which is for /i^ — 20 GeV. Shown in vertical lines are the first four 
resonances for Coulomb potential: En = —E-q/v? where E-q = p\/mt 
and pb = Cpas{fis)iT^t/'^- Note that the self-consistent solution for 
PB = Cpasins = PB)'mt/2 is pB ^ 20 GeV. 



for S wave, which means 



1 y"ir) 



With this and the Schrodinger equation 

A 



2/i 



+ Vir) 



Att r 



where fj, = 1711/2 is the reduced mass, the following quantity can be evaluated in two ways: 

-1 



y"y' 



Thus we have 



-1 



dr{y 



2\/ 







-27r {ijj + rip 







2vr|V^(0)|' 



while 



2/i Jd\[E -V{r)]^^ = Jdr[E -V{r)]{y^y = fi JdrV (r) y'^ = 12 (V) . 



\m\' = ^{V'{r)) . 



(2.47) 



On the other hand, constant shift of the potential, V{r) V{r) + 6V, can be compensated 
by the shift of the energy, E ^ E + 6V. Thus it is related to the shift of the mass mt, since 
E = ^ — Ivrit. This issue is discussed intensively in Section in connection with renormalon 
ambiguity of both pole mass mF°'^^ of a quark and potential V{f) in coordinate space. 



2.6.2 Numerical results 

Analytic formulas for R ratio cannot be obtained with non-zero width Fj even for LO. For LO 
and NLO, the expression for R ratio is 

m = ^NcQ't IniG(r = 0,r' = 0) , 
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where the {S wave projected) Green function G{r, r') is the solution for 



{E + iTt) - 



2d 

dr^ r dr 



Cpas 



+ Vi(r: 



G{r,r') = 6{r-r') 

Anrr' 



with E = ^/s — 2m]i^°^'^\ Note that Vi{r) is zero for LO. For NLO, only 0{l/c) terms are 
used in Vi(r). The same notice is applied for the short-distance current renormalization C[^^'^\ 
It can read from Eqs. ( p. 421) and ( p.43| ); it is unity for LO, and (1 — 2Ciras(/i/i)/7r)^ = (1 — 
ACFasil^h)/Tr) for NLO [Eq. (CT)]. 

The situation becomes more complicated for NNLO (and higher). The last expression is 
the one we actually used: 



2 s m'^c 

2 s m'^c 



^(cur) ^(cur) ~l~ 



lmG'(r, r) 



1 + 



CpCi 



IT 



X Im 
3 



NcQi 



2 

X Im 



1 + 



4m 



1 + 



E + iT 
6mc^ 



1 + 



E + iT 
Qmc^ 



G"(ro,ro 



(2.48) 



CpCLs 



v/r^4^ 



where zq = CpOistTirQ, JI'zq = yuroe'^, and ro = r^^^^ = e'^ /{2mt). Here we used Eqs. 
( |2.41| ), ( |2.42| ) and ( |2.40| ). Reduced Green function G"(r,r') is the solution of 



-1 
rrit 



d? 2d 
dr^ r dr 



CpCts 



+ Vi(r) 



K 



-1 



G"(r, r') = -6(r — r') 



where 



, , 3 E + iF 
a.o = a» 1 + - 



2 m(? 

The coefficient B is determined as 



E^iV 
4mc^ 



B 



9-[r) 



9+{r) 



where g±{r) are the solutions of homogeneous differential equation for G'{r, r') with the bound- 
ary condition 



g^{z) = z"^ [I + z {Lt \Ti\li'z) + Lf Hli'z) + L±) + 0{z^)] . 



Our input is 



mt = 175 GeV , F^ = L43 GeV , af^{mz) = 0.118 . 
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Figure 2.7: i?-ratios for e+e" ^ 7* ^ at LO (dot-dashed), NLO 
(dashed), and NNLO (sohd) as functions of the energy measured from 
twice the pole mass, E = y^ — 2mpoie- Arrows indicate dislocations of 
the maximum point of R as the 0{l/c) and C(l/c^) corrections are 
included, respectively. We put mpoie = mt = 175 GeV, Tt = 1.43 GeV, 
and as{mz) = 0.118. Dotted lines show NNLO i?-ratios calculated 
with an old value of 02 [Q| , which is one of the coefficients in the two- 
loop perturbative QCD potential Vi{r). Figure (a) is for fig = 75 GeV 
and (b) is for fis = 20 GeV. 




Figure 2.8: i?-ratios for e+e —>■ ^* —>■ tt at NNLO for several values 
of tq: vq = a/2 (dashed), vq = a (solid), and vq = 2a (dot-dashed), 
where a = r^'^^'' = /2mt. Figure (a) is for /i^ = 75 GeV, and 

(b) is for = 20 GeV. Other notations and parameters are same as 
in Figure |2.7| . 



2.6 Previously obtained results for R ratio 



35 



There are also another inputs. One is the "hard scale" yU/^ = m^, which is for in matching 
coefficients, and the "soft scale" fig = 20 GeV, 75 GeV, which is for in potential. The other 
is regulator ro = r^^^^ for the incomplete treatment of Tt- 

Shown in Figure |3 is the i?-ratios for LO, NLO (= 0(l/c)) and NNLO (= 0{l/c^)), where 
the last one is calculated in [P5|-p5| for fig = 75 GeV. Our numerical results are consistent with 
them. As was already pointed in p6|-p8|, 0{l/c) and (9(l/c^) corrections for R ratio are not 



small at all. In fact this is the motivation for us to study this issue: how well can we improve 
the perturbative convergence of R. There are three possible sources of this bad convergence: 
(i) after the calculations in |p6H28[|, it was found |^5| that the previous calculation of a 



coefficient 02 in Vi, which is one of the 0{l/c ) corrections, turn out to wrong. In fact the 



original value of 02 is so huge that it was suspected to be wrong before the calculation in ^5 . 
We see the effect of this modification just below in this section, (ii) One can see that the soft- 
scale (fis) dependence of R is also large. This suggests the necessity of Renormalization-Group 
(RG) improvement. We do this in Section |2.8| . (iii) Yet another possible source is renormalon 
ambiguity of pole mass and potential. We pursue this issue in Section p7| . 

Before studying suspects (ii) and (iii) above, let us look at more closely these fixed order 
results for now. In Figures p.7| , the difference of those two figures are the choice of the soft scale 
fig, which is written simply as fi. One choice is /ig = 20 GeV ~ pb and the other is 75 GeV. 
The latter is chosen somewhat arbitrary. Naively, the natural scale for Coulomb potential is 
Bohr momentum — 20 GeV [Eq. (p.T) )]. However one can see that the convergence of the 
normalization is better for fig = 75 GeV than that for /ig = 20 GeV, although the convergence 
of the peak position (~ mass of the IS* resonance) is better for fig = 20 GeV than that for 
Us = 75 GeV. This maybe because the normalization is determined by the wave function 
at the origin, which is much shorter than 1/pb] |V'(0)P(oc {V'{r))). More explicitly, from 



Figure |2]^ one can see that the convergence of the attractive force due to Coulombic potential 
Vc is worse for fi = 20 GeV. There is another attractive force due to Vna for NNLO. This 
also increases the normalization^. On the other hand, the position of 15* peak is determined 
by Bohr momentum p^- However corrections to the potential Vi(r) are minimized neaiQ 
1/r ~ fig = fig^""^- Interplay of these factors makes the estimate of the optimal scale for 



fig complicated. We consider this issue once more in Section |2.8.1| , where the theoretical 
ambiguities of our calculations are examined. 



Also shown in the same figures are the effects of correcting 02 [0,^, which is a coefficient 



of NNLO correction to Coulomb potential. Numerically, the difference between 02°^'^'* ~ 333 
and Og'^*^™'' ~ 156 is large, but the effect of correcting it to R turns out to be small. This is 
because log term in Vi(r) [Eq. ( |2.15| )] are also large. In the following section, we give the results 



for RG improved potential, where large log's are summed. For those cases, NNLO corrections 
diminishes by factor two when the correct value for 02 is used. 

Shown in Figures p.8| are the ro dependence of R ratio. Both vertical and horizontal 



directions are magnified by two, compared to Figures The value of ro is varied by factors 



2 and 1/2. As was explained in Section 2.4.3| , the variation in R is due to not only 0{l/c^) 



corrections but also (9(rj/mj) terms in Eq. ( p.48|) . The size of this variation is a measure of 
one of the uncertainties of our theoretical prediction. We can see it is rather small compared 
to the 0{l/c^) correction itself. Especially, side-ward shift is very small. 



Short distance vertex corrections Ci etc. also affects the overall normalization. For example, with 
afS(mt) = 0.107, it is 1 - ACpaf^inh = mt)/n = 0.818 to NLO. 



19 



For reference, e"*"^ = 1.781 ■ 
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2.7 Renormalon ambiguity and PS mass 



2.7.1 Perturbative convergence 

It is known that the perturbative expansion in QFT is an asymptotic expansion. The ambiguity 
due to this is called "renormalon ambiguity" [^5| -[58[|. This is the subject of this section. 

For some quantities Q, explicit calculations show that perturbative convergence of Q is 
much better when Q are expressed in terms of MS-mass than in pole-mass One 
example is the leading fermionic correction (Ap)j to the p-parameter |49[| : 



^^on-shell ^ N^GpMf 



8^2: 



X 2.9 



TT 



TT 



MS 



8V27r2 



c- 



TT 



(2^-0.19) 



5.01n^ + 15, 



TT 



15 



+ 2.0£-4.0 



where0 Mt 



m 



pole 



mt{fi) = m^^(p), and i = ln[p^/m((p)^] = for p = pt = rrit^fit)- One 



can see that with = fit, the coefficient of a^/vr is smaller by factor 15 for MS-mass, and 
factor 4 in the case (as/vr)^. A similar behavior can also be seen in the QCD corrections to the 



interactions H, W~^W H, and ZZH ||50|. It is suspected that these worse convergence 
may be a signal of renormalon contributions. 

Likewise, as we shall see below, the pole mass m^"^^ of a quark itself suffers from "renormalon 
ambiguity" of O(Aqcd) even if Tt > Aqcd pTH^- For example, the invariant mass rrij^y of 
the jet and W, which are the decay products of t, distributes around m^. However since b is 
not color- singlet, the m^vK-distribution do not directly represent m)'°^'^. 

On the other hand, there is also "renormalon ambiguity" for quark-antiquark potential 
V{r) in coordinate space, and it can be shown [0,|6^ that the severest ambiguity cancel in 
the combination 2m^°'^ -|- V{r). We shall see this just below. 



2.7.2 IR gluon contributions to self-energy E and potential V{r) 
The mass m of a particle / is renormalized due to the self-energy S of / by 6m = mSi, since 



where 



^ — m — — m — mYii — {yi — m)Yi2 

1 i 
~ l-S2y-m(l + Si) ' 



E(j9, m) = mSi(p^, m) + (^— m)S2(p^, m) 



20 



See also i5C 
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Let us consider the one-loop gluon contribution to S of a quark. Since gluon is strongly coupled 
in IR region, we concentrate on those region. For a static quark ~ (m, 0) , 



-zLiR = [-tgs) I I / 77r-Ti7 7 — TTT^ rT~7o 



(27r)4 ' (p + A;)2-m2 + ie ' -|fc|2 + ie 

k = soft 

1+7°^ 2 /" d^'^ 1 

(27r)4(A;0 + ie)(|fc|2-2e) 



-Cf9. 



k = soft 

1 + 7° f d^k ATTCpas 



2 2 J (27r)3 

fc = soft 

I+7OI f d^k 47rC^a,(|A;|2) 



-0 



2 2 J (27r)3 |A;|2 

fe = soft 



where we used /c* ~ in the denominator, and -^o^ = Pp' ~ iirS^k^). Fermion bubbles in 
the gluon propagation alter the | A;|-dependence with the coefficient Uf, the number of (light) 
fermion flavor. By naively non-abelianizing nj to /5o, we have the last expression. Thus we 
have S2 = Si — S2, or 



'^'^IR = TnUi = m ■ 2S2 

1 f d^k AnCpasilkl^] 



2 J (27r)3 |fc| 

k = soft 



where Smin means the contribution from IR gluons. This expression is ill-defined due to the 
IR-pole of the QCD coupling Q;s(/i). On the other hand, the QCD-potential V{r) in coordinate 
space is also ill-defined due to just the same reason: 



(27r)3 \k\ 

soft 



However we can see that |^9|,Q the "severest IR renormalon pole" is canceled in the combina- 
tion 2mpo\e + V{r). The precise meaning of this statement is explained in the next section. The 
origin of the difference of the sign can be understood as follows. Consider an Abelian gauge 
group for simplicity. For the self-energy S, both end of a "gluon" couple to the "quarks" of the 
same charge; thus S oc H-Cira^, where Cp ^ (Q is Abelian charge). For the potential V , 
one end of a "gluon" couples to a "quark" while the other to a "anti-quark"; thus V oc —CpOis- 
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2.7.3 Borel sum and IR renormalon pole 

IR renormalon in mpoie is discussed in IR gluon contribution for a quark mass is 



6mM - 2 J = — I 



, "^"'(2^^ (2.49) 

n=0 ^ ^ 

= — /x/ dte-*/"^(^) (2.50) 

47r 

We can see several points. With the expression in the second line]^, which is written in 
series, the coefficients grows with n\, thus the series is asymptotic one, which means it does not 
converge at all. The procedure that leads to the final expression is called "Borel resummation" . 
One can easily check that this expression is the same as the above one at least formally, by 
expanding the denominator and integrating each term, which gives simply F-function or n\. 
The integration parameter t is called the "Borel parameter" conjugate to ^^(/u). It is sometimes 
convenient to use the rescaled Borel parameter u = ^ conjugate to ■ The complex plane 

of the Borel parameter t or m is called "Borel plane". With the Borel-sum representation of 
5miR, there is a pole on the path of integration, since /3o > 0. The pole in the Borel plane that 
originates from the asymptotic behavior of a perturbative expansion is called a renormalon 
pole. We can see that the pole mass has a IR renormalon pole at m = ^. The following 
expressions may be useful to follow some of the calculations here: 

^ 47r fc2 n=0 ^ 

A, 

where A = Aqcd- The most RHS of the second equation may be called Borel representation 
for the running coupling p7[] . 



Short distance mass and long distance mass 

IR renormalon pole disturbs pole mass mpoie but not MS mass ?n^. This is because MS 
scheme subtracts only the pole that originates from UV divergence. Such a mass scheme may 
be called "short distance mass". On the other hand, a mass scheme sensitive to IR physics 
may be called "long distance mass" . For example, the peak of the invariant mass rrijw of the 
jet and W that are decay products of t, may be a kind of long distance mass of top quark. 
Those mass schemes suffer from IR renormalon ambiguity. 

Here we used J da;^ln^^ = 2"n!. 
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Ambiguity due to IR renormalon pole 

There are several ways to estimate the ambiguity of mpoie due to IR physics. We show three 
ways here. It may be useful to see that both the Borel-summed expression and the asymptotic- 
series one give almost the same result. The symbol A denotes ambiguity. 

(i) In the Borel-summed expression Eq. ( p.50| ), there is a pole on the path of integration. 
There may be two choices to step aside this singularity; one is to around upper side of the 
pole, and the other is lower side. The ambiguity may be estimated by the difference of these 
two |5T|: 



A(5miR(/i)) ~ ilm 



vr ' Jc 1-2 



Iht 

4n 



Note that the power 1/2 of A^//i^ reflects the position of the renormalon pole u = 1/2. Be- 
cause of the exponential dumping factor exp[—t/as{fi)], the severest ambiguity is due to the 
renormalon pole that is closest to the origin. For the case of mpoie, it is u = 1/2. However in 
the combination 2mpoie + V{r), this severest pole is cancelled each other. The next severest 
pole, which may not be cancelled in the combination above, may be m = 1. The ambiguity due 
to this renormalon pole is 

/A^V A 
A {6mm) =A-— forM=l, 

thus suppressed by factor A//x ~ 1/175 for /i ~ mt. Precise determination of the residue may 
be difficult. 

(ii) Let us consider a asymptotic series J^^q 1131 • "^^^ ^^^^ approximation may be 
obtained when one sums up until the ratio of a term and the next term becomes greater than 
unity: [x"'~^^{n + l)!]/[x"n!] = xn < 1, or n < ricrit = '^/x. The ambiguity may be estimated 
by the magnitude to the last term: 



ME 



vn=0 



x"n! 1 ~ x"n!| ~ — 



For the case we deal. 



1 47r ^ ^ 

^crit ~~ — 6 ' 

2/?oas(/i) 



Note that if the severest IR renormalon pole is not m = 1/2 but u = 1, then ricrit = 47r/ {(SqUs 
Thus from Eq. (p.49|) we have 



C 



A {6mm) ~ — ■ V27frwe 



vr 



CFas{fi) 271 

Z /o ( N ^^QCD 
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(iii) There is also an another way to estimate. It gives small imaginary part directly to the 
denominator of the integrand of momentum integration 



A ((5miR(/x)) 



Im [5miR(/i)] 
-1 AnCp 



Im 



Cf^tt 



In || + 



/3o 



-A 



QCD 



For an asymptotic slavery theory /5o < such as QED, the IR renormalon pole is not on 
the path of integration, thus it may not be big problem. In terms of an asymptotic series, each 
term is sign altering, thus the singularity is less than the case for /3o > 0. There are also UV 
renormalon poles in QCD (and in QED as well). However, they are situated at the negative 
region of the real axis (for /3o > 0), thus cause no ambiguity at least with regard to the Borel 
resum. 



2.7.4 Potential-Subtracted mass mps(/i/) 

Arguments by now indicate that one should attempt to extract the MS mass m-y^ directly, 
not the pole mass mpoie, from experiments. However one crucial point is that MS mass is 
not calculated to the order we need. This is because the binding energy (~ —E-q) is already 
0{ma1) to the leading order [Eq. ( p.7|) ]. Thus to NNLO, one needs 0{mal) term, which 
is N'^LO correction to the quark mass. Recently N'^LO correction to the quark mass was 



calculated ||6^. But one more higher order may be very hard to calculate. One way to 
circumvent this dilemma is to introduce another "short-distance" mass scheme. There may be 
a host of schemes, but one scheme which is convenient to our problem is Potential-Subtracted 



(PS) mass mps proposed in ||5^. The argument goes as follows. As we showed above, the 
combination 2mpoie + V{r) is free from the (severest) IR renormalon pole, but each of mpoie 
and V{r) are not. The problematic part is the IR part of \^(|fc|). The idea is to add that part 
to the mass. Then the new potential is free from (severest) IR renormalon ambiguity. Thus so 
is the new mass scheme, since the sum is so: 



where[3 



2mpoie + Vc{r) = 2mps(/x/) + Vb(r; ///) , (2.51) 



mpsifi'f) = mpoic + ^m{fif) , Vc(r; /i/) = Vc(r) - 2Am(/i/) 



and 



Am(M^)sl f -^Vciq) (2.52) 

l+£{..-,„(.nJ-2 



2 J (27r) 

-CpCls 



TT 



2 

s 

+ (£)^{a.-(2a.3„.ft)(.ng-2).,3j(l„=g-4..|. 



Note that our Am{iif) is related to a corresponding quantity in [5^ by Am{fif) = —5ni{^f), and is 
negative. 
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Am(/i/) 


/i/ = 3 GeV 
/i, = 20 GeV /i, = 75 GeV 


/i/ = 10 GeV 
/i, = 20 GeV fXs = 75 GeV 


LO 
NLO 
NNLO 


-0.1922 GeV -0.1546 GeV 
-0.3123 GeV -0.2589 GeV 
-0.4003 GeV -0.3389 GeV 


-0.6408 GeV -0.5153 GeV 
-0.8945 GeV -0.7707 GeV 
-1.035 GeV -0.9255 GeV 



Table 2.1: The difference Am{fif) = mps{n j) — rupoic of PS mass and 
pole mass, defined in Eq. (|2.52D . Although we use ^/ = 3 GeV in the 
study, we also show Am{pf = 10 GeV) for reference. 



where Vc{q) is the Fourier transform of the Coulombic potential Vc{r) defined in Eq. ( p.20|) , 
and as = a^^(/i<j). The scale /x/ for factorization should be larger than Aqcd in order to 
remove the IR ambiguity. On the other hand it seems that fij should not be as large as Bohr 
momentum in order not to disturb the Coulombic potential near Bohr radius. Following 
the original paper ||5^ for PS mass, we choose fJ^f = 3 GeV. Explicit values are given in 
Table |2.1| . For reference, Am{fif = 10 GeV) is also shown. Note that in the definition above, 
we dropped the r-dependence of Am(/i/). Thus the equality in Eq. ( ^.51| ) is up to 0(AQQj-,r^) 
at least perturbativelyp^. There may be renormalon ambiguity also in this r-dependent term. 
This ambiguity can not be cancelled in the combination 2mpoic + V{r) since the pole mass is 
r-independent. However the ambiguity in this term may not be serious, because the typical 
scale r for (would-be) ti bound state is 1/pb — 1/(20 GeV), which means the ambiguity is of 
order A^^j^r ~ 0.05 GeV. It is known that V{r) oc r outside the perturbative region. However 
as we shall see in Section |2.8.1| , this ambiguity is negligible in our case due to the large mass 
and decay width of the top quark |^, |TT], ^ , n | . Thus the next-severest renormalon ambiguity 
above may also be negligible. 

Recently, the relation between MS- and pole-mass was calculated to C(«^) |Q. With this, 
PS-mass can be related more accurately to MS-mass than to pole-massQ. 



2.7.5 Results for R ratio with mps(/i/) 

Following the arguments above, it is PS mass (or a short-distance mass, in general) that is 
well determined from experiments. Thus we fix PS mass but not pole mass. With the choice 
fif = 3 GeV, the difference Am(/ij) ~ —0.4 GeV of pole mass mpoie and PS mass mpg is 
tiny compared to mpoie — 175 GeV. Thus the Schrodinger equation barely changes except the 
origin of non-relativistic energy E, since it is now E = ^/s — 2mps. Note that, as explained 



in Section 2.6.1, constant shift of potential V simply means the shift of the origin of E. 



Shown in Figure ^]9| are almost the same as those in Figure |2.7| except we fixed mps(/i/ 



3 GeV) = 175 GeV rather than mpoie- Since |Am(/ij)| is larger for higher order, the position 
of IS" peak moves more for higher order. We can see the sideward convergence becomes better 
by the difference |(Am)NNLO - {Amf^^\ ~ 0.1 GeV and \{Amf^^ - {Am)^^\ ~ 0.1 GeV, 
respectively. On the other hand the normalization barely changes, as is anticipated. In the 
next section we study RG improvement of the potential, which is expected to be effective to 
this problem. 



There is no Oi^Aq^j^r) term due to rotational symmetry. However I'm not sure if this is true also non- 
perturbatively. 

See the reference above for details. 
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Figure 2.9: i?-ratios for e+e~ j* ^ tt at LO (dot-dashed), NLO 
(dashed), and NNLO (sohd) as functions of the energy measured from 
twice the potential-subtracted mass, \/s — 2mps. We set fJ.f = 3 GeV 
and mpsifJ-f) = 175 GeV. Figure (a) is for /x<j = 75 GeV, and (b) 
is for Us = 20 GeV. Other notations and parameters are same as 
in Figure ^.7| . 



2.8 Results for R ratio with RG improved Vc and mps 



As mentioned in Section |2.6.2| , large renormalization-scale /ig dependence of R ratio indicates 



the necessity of log(/is) resummation. As first step[^, we sum up log(yLis) in the Coulombic 
potential Vc(r). Thus our results here still depends on the soft scale /i^, which is in Vbf and 



V^A- As we saw in Figure |2.4| , couplings converge much better for RG-improved case than for 
fixed order, all over the relevant momentum scale Cpdsf^t ^ V''" < "^t- Thus it is expected 
that the convergence is improved by the use of V^^'^'^ (r) that is RG improved (/i^ = q) in the 
momentum space and Fourier transformed to the coordinate space. 

PS mass mps{iJ,f) with RG improved potential is ambiguous, since V^^^^ diverges at IR 
region. Here we define as follows: 



|91>M/ 

|gl<M/ 



Our final result for i?-ratio in this paper is shown in Figure 2.10. For comparison, /2-ratios both 



with and without the RG improvement for Vq are given. We can see that convergence of both 
the normalization and the peak position is improved slightly. The latter is listed in Table |2.2| 
as the "binding energies" of the IS" resonance state 2mps(/i/) — Mis. They are determined by 
reducing the decay width F^ tentatively. 

Why is the convergence not improved much better? Let us drop the corrections other than 



those for Vc{r), for the moment. The R ratio for this case is shown in Figure |2.11|. The relevant 



A full resummation of logarithms up to NNLO requires a significant modification of the formula Eq. ( 2.48|) ; 
we will study its incorporation in our future work. 
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Figure 2.10: Our final result for i?-ratios in this paper for e"'"e~ — > 
7* ti at LO, NLO, and NNLO. Solid lines show those with 
renormalization-group improved Coulombic potentials, {r;iif). 
Dashed lines are those with fixed-order Coulombic potentials Vc(r). 
Arrows indicate dislocations of the maximum point of R as the 
©(l/c) and 0{l/(?^ corrections are included, respectively. We set 
lif = 3 GeV, mps(/x/) = 175 GeV, jis = 75 GeV, Tt = 1.43 GeV, and 
as{mz) = 0.118. 



(fixed-order) (RG-improved) 
At, = 20 GeV = 75GeV /x, = 20GeY /x, = 75GeV 



LO 

NLO 

NNLO 



1.390 GcV 
1.716 GeV 
2.062 GeV 



0.838 GeV 
1.453 GeV 
1.817 GeV 



1.573 GeV 
1.861 GeV 
2.136 GeV 



1.573 GeV 
1.861 GeV 
2.058 GeV 



Table 2.2: "Binding energies" of the IS* resonance state defined as 
2mps(/U/)-Mi5 at LO, NLO, and NNLO calculated with Vc{r) (fixed- 

order) and with Vq {r;^Xf) (RG-improved). These are determined 
by reducing the decay width Vt tentatively. We set jjLf = 3 GeV, 
mpsifif) = 175 GeV, Ft = 1.43 GeV, and as{mz) = 0.118. 
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Figure 2.11: i?-ratios for e^e^ j* ^ tt calculated with a Hamil- 
tonian H = p'^/mt + Vo(r), where Vo(r) includes only the Coulombic 
part of the tt potential. Other corrections (hard vertex corrections, 
relativistic corrections, etc.) are not included. Solid and dashed 
lines, respectively, show ii-ratios with (Vb(r) = {r;iif)) and 

without (Vo(r) = V^^^'^\r; ^f), fi^ = 75 GeV) a renormalization- 
group improvement of the Coulombic potential. We set fJ-f = 3 GeV, 
"^Ps(/U/) = 175 GeV, Tt = 1.43 GeV, and as{mz) = 0.118. 




fixed-Older (|I = 75GeV) 



formula is 



R^s) = ^'^^f' lmG(0,0) 
mf 



with 



-1 
rrit 



dr^ r dr 



+ Vo{r)-uj\G{r,r') 



ATirr 



both for fixed order Vo(r) = V( 



(fixed) 



r; 



jj,f) and RG-improved Vo(r) 



We can 



see clearly that the convergence is much improved by the \og{iis) resummations. Thus it is 
other corrections that disturb the perturbative convergence of R ratio. Among them the 1/r^ 
potential V^x{r), which is a strong attractive potential near the origin, is suspected to be 
the worst. It remains as our future task to gain better understandings of these residual large 
corrections. 

Shown in Figure 12?!^ is a^^(m^)-dependence of i? ratio. The peak position moves 0.12 GeVx 
2 with the variation of Q^^{mz) = 0.119 ± 0.002, which is the current uncertainty M. 



2.8.1 Theoretical uncertainties for NNLO calculation of R ratio 

There are several uncertainties in our NNLO calculation of R ratio. In this section we study 
the size of them and estimate the theoretical uncertainty (Am^)*^ of top mass determination. 
Since the perturbative convergence of R ratio is the best with the RG-improved potential and 
the Potential- Subtracted mass, we study the uncertainties for this case. 

Leading log(g/yUs)'s of Vc{q) can be summed over not only with the choice fis = 1, but also 
= 2g etc. Their difference is related to the summation of the next-to-leading log's (and 
higher). Shown in Figure p.l3| are the uncertainties due to this freedom. One can see that it 
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Figure 2.12: Dependence of R ratio with the variation of a^^(m^): 
0.117, 0.119, and 0.121 for the left, and 0.108, 0.118, and 0.128 for the 
right. For both figures, IR part A: < A = 3 GeV of Coulombic potential 
Vc(A;) is not taken into account for the Fourier transform to Vc(r). For 
the figure on the right, twice the Ani(iJ, f = 20 GeV) — Am (A = 3 GeV) 
is further subtracted from Vc{r) and added to mps- No width is 
included for as{mz) = 0. We set Hs = 75 GeV, m,ps(/i/) = 175 GeV, 
Ff = 1.43 GeV. 



is fairly small for both NLO and NNLO. The dependence is large for LO, since there is no log 
that compensates the log dependence of the coupling. 

The soft-scale dependence still remains in the potential, since we RG improve only the 
Coulombic part. It is shown in the left-half of Figure p.l4| . A physical quantity Q do not depend 
on the choice of renormalization scale /x provided Q is calculated to all order. If perturbative 
expansion is terminated to certain order, Q depends on /i since the discarded higher order 
terms depend on /x. Thus it is expected that the truncated perturbative expansion is a "good" 
approximation when dependence is small, and the magnitude of dependence around there 
is a measure of (uncalculated) higher order corrections. With this view, one can see that the 
choice /is ~ pb — 20 GeV is not "good" for R ratio. The right-half of Figure |2.14| shows cut-off 
To dependence, which is absent for NLO and LO, and also for NNLO if one treats the effect of 
top decay width Ff properly. One can see that the dependence on tq is small for tq ~ rQ^^\ 
Thus one can see both /i^ and tq dependence is small compared to the difference of NNLO 
R ratio and NLO R ratio. Besides the shifts are only horizontal, thus do not disturb the 
determination of top mass mf^(/ij). 



More quantitatively, the variation of the peak position of cr, 
= 20 GeV-150 GeV and /i,/g = 1/2-2. 



NNLO 
tot 



is 0.08 GeV x 2 both for 



Estimate of (Am^) 



th 



Now let us estimate the theoretical uncertainty of top quark mass mf (/i/) with /i/ = 3 GeV 
based on our NNLO calculations. Since both /i^- and ro-dependence are fairly small in the side- 



ward direction, we can forget about them. RG-improve prescription dependence [Figure |2.13 
may serve as a measure of theoretical uncertainty. One can see that the sideward NLO cor- 
rection to LO i?-ratio is similar to the uncertainty of LO i?-ratio due to the choice of ^s/q- 
Likewise, NNLO correction of NLO R is also similar to the uncertainty of NLO, but it is a little 
smaller. This may be because there are sources of NNLO correction other than those log's; that 
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Figure 2.13: Dependence of R ratio on the prescription to sum-over 
leading log(/ic,./g)'s of Vdq)- Three choices are given for each with 
NNLO (top), NLO (bottom-left), and LO (bottom-right) corrections: 
lis = q/'2., q, 2q. We set /x/ = 3 GeV, nipsifif) = 175 GeV, Tt = 
1.43 GeV, and as{mz) = 0.118. For NNLO /x^ = 75 GeV. 



I* 0.8 




-2 

2 mps [GeV] 



1.6 
1.4 
1.2 
1.0 
« 0.8 
0.6 
0.4 
0.2 
0.0 



NNLO, RG improved V^r), n = 75GeV, Hf= 3GeV 




-2 2 

2 njps [GeV] 



Figure 2.14: (Left) The soft-scale (ns) dependence of R ratio with 
NNLO corrections, iig is written as fj, for simplicity. The normalization 
changes monotonically in the choice of shown in tlie figure: fig = 
20 GeV for highest and fig = 150 GeV for tlic lowest. (Right) The cut- 



off (ro) dependence of R ratio with NNLO corrections. 



IS written 



as a for simplicity. In both figures, we set fif 
175 GeV, Tt = 1.43 GeV, and as{mz) = 0.118. 



3 GeV, mps(/x/) 
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is, Vbf and Vna- Thus we estimate the NNLO theoretical ambiguity to be {ArritY^ ~ 0.1 GeV. 
Here we are conservative, since there are corrections aside from log's, such as VVa- 

As for N'^LO = 0{l/c^) corrections, those terms for Vbf is ~ Cpa's/im^c^f^'^), since mcr 
is dimensionless. Here only the dimension is considered. Thus translations like ~ p^, 
~ 5^^\r), 1/r ~ d/dr are possible. Likewise, 0{l/c^) correction to Vna may be ~ 
CACpal/iTn'^r^c^), since may be 0{l/c^). Another significant modification is that real 
gluon emission is possible to 0{l/c^). 

In our calculations, leading log's in Vq are summed with RG. Recently next-to-leading log's 
corrections of N'^LO, ln^(l/a<(), to R ratio is calculated in One can read that the size of 
N^LO corrections to the peak position is ~ 2 x 0.1 GeV, which is consistent with the estimate 
above. 



Non-perturbative effect 

Non-perturbative corrections modify QCD potential at r ~ 1/(1 GeV) and far. Thus their 
effects can be estimated by the dependence of Green function G with respect to the variation 
of the potential around r ~ 1/(1 GeV); and it turns out to be small. For r > r', 

G'(ry)- ^ 9>ir)9<{r') 



AirCpasW rr' 

where g>{r)/r and g<{r)/r are the solutions of (reduced) Schrodinger equation [Eq. ( p.38| )] 
with the boundary conditions limr-^oo fi'> (^) = and \im.r-,og<:{r) ^ [CFdsfTit^Y^ ■, respectively. 
Since only G'{r,r')\r=r'~Q is relevant for R ratio, it is g>{r) near the origin that determines R 
ratio. It may be affected by non-perturbative effect (r > 1/(1 GeV)) through the boundary 
condition. However we can see in Figure p.l6| that gy{r) dumps well within perturbative region 
r/r-Q < 20 with r ~ 1/(1 GeV). Thus non-perturbative effect (to the QCD potential) do not 
alter the Green function for ti. We also show (7<(r) in Figures p.l7| and |2.18| , for reference. 



Effect of Z exchange 

To (9(l/c^), axial- vector current also contributes since it is suppressed by C(l/c). Note 
that the "cross term" j^jy of axial- vector and vector (jy) current do not contribute to total 
cross section Section ^.3.2|. The effect of Z exchange between e^e~ current and ti current is 



incorporated by rescaling the vector-vector coupling and by adding the effect of axial- vector- 
axial- vector coupling. With the notation defined in Eq. (|A.33|) , the former changes only overall 
normalization by the ratio -|- |[a^f*]p to Q^. The effective couplings [f^f*] etc. are 

defined in Section A-.5| . Numerically it is 0.480/0.444 = 1.080 on the threshold. Especially, 



theoretical uncertainty (Am^)*'^ for determination of rrit is not altered by this inclusion, since 
it is energy independent. While in it is shown that the effect of i\iA is small and flat with 
respect to E near the threshold^, since it is suppressed by compared to jyiv- Note that 
only the leading order calculation of j\jA is sufficient for the NNLO calculation of R ratio. 
Thus R ratio near the threshold is not altered when the effect of Z exchange is taken into 
account, and so is (Amt)*'^. 



Note that the correspoding coupUng is |[i;'^a*]p + |[a'^a*]p = 0.1430 on the threshold. 
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Figure 2.15: The function g>{r)/r is a solution of the (reduced) 
Schrodinger equation to NNLO [Eq. ( 2.3^ )] that dumps at r — > oo. 
CM energy is changed from E = —5 GeV to OGeV. They are plotted 
vs. r/re = PBf, where l/re = Pb = CFasmt/2 = 17.94 GeV for 
tts = a^^{^s = 20 GeV) = 0.1534. Perturbative potential description 
is valid within the hard scale ~ ps/"^* = 0.1023 and the Aqcd scale 
~PB/(lGeV) = 17.94. We set Hs = 20 GeV, fif = 3GeV, mps(^/) = 
175 GeV, Ft = 1.43 GeV, and a^(mz) = 0.118. 
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Figure 2.16: The same as Figure 2.15 . CM energy is chaged from 
E = 1 GeV to 5GeV. The last two figures are both for E = 5GeV; 
only the range of r differs. 
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Figure 2.17: The function g<{r)/r is a solution of the (reduced) 
Schrodinger equation to NNLO [Eq. ( |2.38 )] that dumps at r ^ 0. 
Other situations are just the same as Figure 2.15. 
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Figure 2.19: R ratio in m-poie scheme 
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(left) and ra-y^ scheme. See 



Effect of using nxy^ instead of mpoie 



Shown in Figure p.l9| is the effect of using different mass schemes. In the left figure, pole-mass 
is fixed to mpoie = 174.788 GeV; while in the right one, MS-mass is fixed to fnim) = 165 GeV, 
which corresponds to mpoie = 174.788 GeV at NNNLO. In both figures, a^^{mz) = 0.119 is 
used; this corresponds to a^^(165 GeV) = 0.109. Note that the used potential is one in fixed- 
order; also PS-mass scheme is not used. By letting — > 0, one can investigates the position 
of IS* peak. These positions are also shown in the figures by vertical lines. Numerically, they 
are (348.321 GeV, 347.504 GeV, 346.986 GeV) for (LO, NLO, NNLO) when mpoie is fixed, while 
(344.408 GeV, 346.892 GeV, 347.340 GeV) when m(m) is fixed. 



Chapter 3 

Top-Quark Momentum Distribution 



Due to the large decay width Fj, the momentum \p^\ of top quark is not determined by the 
CM energy y/s = 2mt + E, but distributes around \>ymt{E + 1 GeV + iTt)\, where 1 GeV ~ 
binding energy. Momentum distribution da /dp of top quark in e"^e~ ti near the threshold 
is calculated in |13[ to LO, and in []T5[ to NLO. In this section we calculate a part of NNLO 
correction; that is, only gluon exchange between ti is considered. Gluon exchange between t 
and 6 etc. are not considered. 



3.1 Derivation of the distribution 



Differential cross section of tt pair production can be written as follows: 



da(e+e- ^ tt ^ bW+bW 

E 



Nr ^ (2mt)4 



2s 



X,Y='y,Z 



[s — mx){s — itly) 
1 



pj — mj + irritTt pj — mj + irritTt 



d<^i{bW+bW~ 



where 



{Exvr = tr 



AM 



;i-Pe+Pe-)-(Pe 



9w_ 
V2 



tr 



i + -p^ -Pt + m 1 - 75 . 1 - 75 



2m 



X 



2m 



7 



7^ 



X 



-^^^-^^^^1-75^^,1-75 
^ ' Y 7 



X 



2mt ' 2mt 2 
{Pw)a{pw)i3 



X 



-9af3 + 



m 



w 



-9pa + 



2 

{Pw)p{Pw)a 



m 
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and 



AM 



9x [v'^'j 



eX 11 n^-^^ii^. 



a 



x = (|PtrM-(^ + *ro)x 



This can be obtain by usual Feynman rules with the modification of the ti production vertex 
due to the Coulomb rescattering. The formula above can be reduced more by approximating 
the momenta of top and anti-top by their on-shell valuesQ. By decomposing the phase space 
of hWnW- into tt 



d^i{hW+hW- 



d*^p 
(2^ 



d$2(t* ^ hW+) d$2(t* ^ hW- 



where pt = g'/2 + j9 and pt = q/2 — p {q is the total momentum), both bW~^ and bW can 
integrated out: 



{Pw)a{pw)l3 



m 



w 



= -75). 
mt 

As for the remaining phase space d'*p/(27r)^, the time-component can be integrated with the 
formulas 



1 



27r z — {x + iy) z + {x + iy) 2{x + iy) ' 

■J 1 T 2 .2 

dz 
2^ 



1 1 


2 




2; — (x + iy) 2; + (x + iy) 




2(x + iy) 



since we put 
we have 



ypl + mf everywhere aside from at the denominators of propagators. Thus 
d(T(e+e- ti) 



1 N, 



c 



(Pp 



E 



2r, 



(27r)3 2s ^^^^^^ (s-m^)(.-m|-) 



1 



E - p'^/mt + iVt 



X tr 



2m/ 



2mt 



XY)fiu 



The momentum distributions of the decay products of W^s depend on the polarization of 
top-quark. For such calculations, the following relations may be useful: 



'-I- m ^ +m_^+ml — |*75 „ 

-Gr = G 



2m 



2m 



2m 



^ — m^p — m 
2m 2m 



2m 



m 1 zl^c 



^ Except at the denominators of the top and anti-top propagators. 
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where 



1 - P-s 



C = tr 



m 1 — 
"2m 2 



1 - P-s 



C = tr 



i'— m 1 — 
"2m 2 



There relations hold when = rn^, P-p = s-p = 0, and the likewise for p. These relations 
can be easily verified at the rest frame of p. Note that although G is a general 4x4 matrix, 
only 2x2 sub-matrix contributes because of the projection operators ^ + m/(2m) etc. Those 
four degrees of freedom are expressed by one scalar C and one 4- vector P with one constraint 
P-p = 0. 

For the case of 7-mediated tt production, the distribution da /dp can be calculated in usual 
way with the vertices 

A'^ = eQ^r , - (\Pt\Vm -{E + iP,)) eQtG{E, p,)^ . 

At the leading order, the distribution can written as 



1 da 



4V 3 2 4^2 An 

'J^c^t ~^ * 



s m^c 



G{r = 0,p) 



apt dp (27r)3 2 

For Pt — > 0, this is zero unless G{r = 0,p) diverges. This is consistent with 

1 da 



= S{p - ^/^) , for Ft = 0. 



3.2 Unitarity relation 

Momentum integration of the differential cross section da jdp should coincide to the total cross 
section a^^t- This is relation is called unitarity. It is important to confirm oneself that unitarity 
holds both theoretically and numerically within the accuracy considered. 
For a Green function 



G^ 

its imaginary part can expressed as 

ImG^ -(G-Gt) =Gt 



(G-^) 



1^ 



G-\ ^ ^(G-i)t - G- 



2% ' 2i 2i 
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Thus the relation 
da 
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G{\p\,E) =lmG{x = 0,E) . 



G{\p\,E) 
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da 
dp 



aptA^cQ? opTt 



irmr 



g{\pIe) 
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is consistent with the expression for the R ratio in Eq. ( p.3| ). 

Including the higher order corrections for tt current, the unitarity relation above becomes 



Im 



G(r,r') 



^(cur) . ^(cur) "I" ^r' 

/^(cur) _^ ^(cur) + ^r' 



(2vr) 

d^P /^(cur) ^(cur) + A,,./ 
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Thus the momentum distribution is expected to be expressed as follows: 

G{r,p)G*{r',p) 



s m^c 



apt dp (27r)3 2 

With the relation 

A 
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2m2c^ 



r,r'— >ro 



— G'(r,p) = - 



Gpcis 



G{r,p) - 



sin(pr) 
pr 



+ - 



which is just the Schrodinger equation with the Coulomb potential, the derivatives can be 
removed: 



x-y(cur) ^(cur) A^ -|- ZX^.' 



^(cur) _ 



Gf /a,\2 1 



G{r,p)G{r',p) 
E 
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3 V c / masTQ Sm'^c^ 
1 sin(pro) 
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3mc2 pro 



Re 



Giro,p) 



Here the equality is up to 0{l/c^). We want to express the unitary relation in terms of the 
reduced Green function G": 



G = G' + 6G' + 0{ ^ 



where SG' = - {bcu + Oe , G'} - [Oq , G'] - b . 



Thus with 



G{r, p) = G\r, p) + 5G'{r, p) + 0[- 



where G{r,p) = (r| G \p), we have 



G{r,p) = G'{r,p) 



+ 2 Re 



6G'{r,p) G'{r',p)* 



+ 



Differential operators can be rewritten in terms of c-numbers as we did for the total cross 
section: 



1 d 

ma.r dr 
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Thus with the definitions 
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we have, up to 0{l/c^), 
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Here we used Eq. (|2.42| ). Up to the terms that vanish when ro ^ 
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Figure 3.1: Top quark momentum distributions at LO (dot-dashed), 
NLO (dashed), and NNLO (solid) for /x, = 20 GeV (left) and fis = q 
(right). For each curve, we set the cm. energy on the IS resonance 
state, ^/s = Mis- 



3.3 Results for da /dp 



To summarize, the momentum distribution of top quark near the threshold can be written as 
follows: 



where 

/(P, ^o) 



1 + ^) |G"(p,ro)P + ^Cp«.(/x.)2Re ^^/.(p, tq) G"(p, tq)* 



^CpasifisY Re G'ip^{p,ro)G'{p,roy 



+ ^^^Re 



G"(p,ro) 



Qrrit pro 

Here p = \p^\. We changed a notation of momentum-space Green functions a little bit: 

G{p,ro) = jdh e*^^-^"G(r,ro) , 

1 



Gi/rip, ro) 



dV e^P-"^ 



as{iJ,)mtr 



G{r, ro) 



G.p.(p,ro) = /dV e'P-^^ ^^G{r,ro) . 

The unitarity relation holds numerically, to the degree we are concerning. 
Our choice for the soft scale is /x^ = 20 GeV since a relevant scale around the distribution 
peak is the scale of Bohr momentum pb- 



Shown in Figures |3.1| and |3.2| are the top quark momentum distributions. We normalize 
the distribution to unity at each distribution peak, because overall normalization do not have 
extra information other than the total cross section R has. This distribution depends on CM 



energy. The Figure |3]^ is on IS resonance and Figure ^]2] is 4 GeV above the IS". We use the 
notation AE = \/s — Mis- One can see that corrections are fairly small compared to those for 
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Figure 3.2: Top quark momentum distributions at LO (dot-dashed), 
NLO (dashed), and NNLO (sohd) for /x^ = 20 GeV (left) and mus = q 
(right). For each curve, we set the cm. energy at 4GeVabove the IS 
resonance mass. 



total cross section R. Using RG improved Coulombic potential also do not change qualitative 
feature. Both 0{l/c) and (9(l/c^) corrections are similar qualitatively, while quantitatively, 
the latter is larger. The most prominent correction is the shifts of peak momenta Ppcak for 



AE = 4 GeV Figure 3.2. However this is mainly due to the shift of IS* resonance, which 



is already visible in R ratio. That is, for energies AE > 1-2 GeV, the peak momentum of 
the distribution tends to be determined only by kinematics, Ppeak — |a/s — 4mj . Note that 
binding energy becomes deeper with higher order correction Table p.2| . Correspondingly, the 
peak momentum in Figure p]2| lower with higher corrections. On the other hand, this kind 



of correction is small on the 15* resonance (AE = 0, Figure |3.1|) . Quantitatively, 0{l/c) and 
C(l/c^) corrections shift Ppeak is 5ppeak/Ppeak = —0.8% and +2.5%, respectively, for fixed order, 
while +0.5% and +2.2% for RG improved. Behavior in intermediate energy is as follows. At 
AE = the corrections are positive ~ +few%; between AE = and AE = 1-2 GeV, 
the corrections decrease and change sign from + few % to — few %; at higher energies, AE > 
1-2 GeV, the corrections stay negative, but their magnitude |5ppeak/Ppeak| decrease with energy. 

Aside from the correction mentioned above, the next prominent correction is for higher 
momentum region. This is because higher order corrections in 1/c-expansion are less suppressed 
for higher energy and momentum. 



Other corrections 

Although the momentum distributions above are calculated with (9(l/c^) corrections, they are 
not complete. This is mainly because we do not treat the decay process of t and i properly. In 
Section |2.2.3| we mentioned that treating the decay width F^ with the replacement E —>■ E + iV 
is valid for 0{l/c) and below. In fact, this is only for total cross section R. Final-State 
Interactions (FSI) between the decay products of t and i modify the momentum distribution 
of t already at C(l/c). In the analysis above, we dropped this effect. The effect of FSI to 
0{l/c^) is not calculated yet. Shown in Figure |3]^ are 0{l/c) momentum distributions with 
and without FSI. We can see that the effect of NLO FSI is somewhat larger than that of 
NNLO potential modifications etc. Quantitatively, the final-state interactions reduce the peak 
momentum about 5% almost independent to the energy [|^. Note that the energy dependence 
of the corrections are different for NLO FSI and NNLO potential modifications etc. 
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Figure 3.3: Top quark momentum distributions at NLO with the 
renormahzation group improvement for the Coulomb part of the po- 
tential. The cm. energy is set on the 15 resonance state (left) and 
4 GeV above it (right). The solid (dashed) line is calculated with 
(without) the 0{l/c) final-state interaction corrections. 



Chapter 4 



Anomalous Electric Dipole Moments of 
Top Quark 



In the early stage of studies on e^e — tt near the threshold, the effects of a Higgs exchange 



between tt were also examined |]TT|. It showed that the mass thh and the coupling g^^jj of 
the Higgs can be probed, if the Higgs is light. Many of extensions of the Standard Model 
(SM), including Supersymmetric SM (SSM), contain extra Higgs multiplets. They in general 
introduce complex couplings that are unremovable by field redefinitions, or "physical", thus 
provide sources of CP violation other than those in CKM matrix. Not only in extra Higgs 
multiplets, but also in other extensions, there are plenty of sources for CP violation in general, 
once we expand the particle content of a model. For example in SSM, mass matrices for 
squarks and sleptons can have complex phases. These extra sources for CP violation induce 
EDM for various particles to one-loop level, while the leading SM contribution is three-loop. 
The magnitude of EDM of a fermion is proportional to its mass, in general. Thus we study 
the EDMs of top quark in this section. Because EDMs are CP violating, their effect can be 
seen in the difference of the polarization of t and t, which is zero if all the relevant interaction 
is CP conserving. 

In Section ^Tl] it is explained that a non-zero value of (intrinsic) EDM violates both P and 
T. Field theoretic description of EDM is also given. In Section ^]2| expected values of EDMs 
in various models are summarized. 



4.1 Electric Dipole Moment (EDM) 

Total angular momentum is a sum of orbital angular momentum and (intrinsic) spin. Just the 
same situation holds for Electric Dipole Moment (EDM) and Magnetic Dipole Moment 
(MDM) 

f^E = f^EDM + D, flB = f^MDM + M. (4.1) 

Here we mean intrinsic dipole moments by tJ^EDM Mmdm; "orbital" ones by D and M. 
As is well known, the transformation properties of electric field E etc. under, say, parity can 
be defined so that the classical electromagnetism is invariant under any combinations of Parity 
P, Charge conjugation C, and Time reversal T. However this is not the case for quantum field 
theories, as we shall see below. 
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P 3 


E B 


D M s X p 


p 


+ - 


- + 


- + + - - 


c 








T 


+ - 


+ - 


+ - - + - 



Table 4.1: Transformation properties under Parity P etc. of the quan- 
tities related to classical electromagnetism; p is charge density, j cur- 
rent density, E electric field, B magnetic field, D "orbital" part of 
electric dipole moment, M "orbital" part of magnetic dipole moment. 
These are fields, s is spin (or angular momentum, in general), x is 
coordinate, p is linear momentum. P etc. mean the eigenvalues under 
P etc. For P and T, arguments of fields are changed from {x,t) to 
{—x^t) and (a;,— t), respectively. 



The EDM and MDM ^tg are defined by the following interactions: 

V = -fiE-E - tXj, B . (4.2) 

Within the classical electromagnetism, only the "orbital" contributions {D,M) to the dipole 
moments can be calculated: 

D{t) = jd^xxp{x,t) , M{t) = ^ jd^xxxj{x,t) , 

where 

p(a;,t) = ^g,5(%-s,(t)) , j{x,t) = J2(1^^^S^''\^ - Xi{t)) , 

i i 

are charge density and current density, respectively. From these definitions, we can see the 
transformation properties under Parity and Time reversal. For example, 

p{x,t) - {-Xi{t))) = +p{-x,t) , 

i 

because (5-function is an even function, and then 



Dwi/d=..(-.)p(-..o^-z,w. 



Here we changed the integration variable x ^ x' = —x. The transformation properties of 
other quantities are given in Table |4.1| . 

Within the framework of QFT, spin s is the only vector available for a particle at rest, thus 

Medm ! Mmdm // ^ • 



Thus ^j,-^Y)M'B is (P,T) = (+,+) while ij,-^y)m'E is (P,T) = (—,—), which means non-zero value 
of (intrinsic) EDM ^tEDM is allowed only when both P and T (or CP) are violated. 
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As is well known, wlien tlie interaction of fermion-fermion-gauge boson has the form of 

2m 

MDM is expressed as 

where m is the mass of the fermion, Fi{q'^) and F2{q'^) are form factor^, which depend on the 
momentum of the gauge boson (flows into the fermion current), and g is Lande (/-factor: 
{g — 2)/2 = a/(27r) + ■ ■ ■ . Thus anomalous MDM is expressed by the second term of the 
following effective Lagrangian: 

C = -eF,{q')i/jY^A^ - -^F2(g2)^a^-^a^A, . 

2m 

From this we can obtain the corresponding term for EDM by replacing B hj E. This can be 
done by using dual field strength: 

FfMu = (Foi, —eijkFjk) = (E, B) , —^ t^^paF'"' = {—B, E) . 
Our convention for the totally antisymmetric tensor is 69123 = +1. Thus with 



we have 



r'^ = ^^F,{q') + ^a'^-'q. [F,iq') + ^75 diq' 
2m 



or the following effective Lagrangian for EDM: 



C = -eF,iq^)i,^^M>. - ^^i'ta'-'^.ijdpA, 

2m 

~ -eF,iq'')i,YMf. + ^-^{i^il,id^ij)A, , (4.3) 

2m 



where d{q^) is a dimensionless form factor for anomalous EDM: 

ed <7 



The "~" in the last step means they coincide with the use of the eqs. of motion [Section |A.4 
The operator \^ia^'''-i^^l}F^,y ~ -{^i-i^id^^)A^ has (P,T) = (-,-) [Tables Q The 
anomalous EDM interactions with Z and gluon g are defined in the same way. Since the 
interactions are flavor-diagonal, its existence immediately means CP violation. 

In fact, as shall be shown in Section [A.4| , the EDM-interaction above, which is dim-5, is the 
lowest dimension operator for t-t-gauge boson that violates CP invariance. Thus CP-violating 
effect is taken into account by the EDM to the leading order. The form factors (i(g^) depend on 

^ Fi{q'^ = 0) is the electric charge of the relevant particle. 
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the momentum transfer in general. However to the lowest order of derivative expansion, it is 
constant. To our knowledge, all studies to date concerning sensitivities of collider experiments 
to the EDMs assumed that the momentum dependence of d{q^) can be neglected. Also, as 
we shall see in the following section, the g^-dependence is shown to be moderate, by explicit 
calculations of EDMs based on specific models like Minimum Supersymmetric Standard Model 
(MSSM). Thus we also follow this assumption. The magnitude of EDMs are 

tt7-EDM = ^ = dt^ X 1.13 X 10"^^ecm , 
mt 

ttZ-EDM = = rf;^ X 1.13 X 10^1 V cm , (4.4) 

nit 



ttg-YiDM = = dtn X 1.13 x 10"^% cm 
mt 



Note that gz/e = {g/ cos 6'vi/)/(fi' sin 6*14/) = l/(cos^i4^ sin ^14^) = 2.372 ■ ■ ■ . 

As can easily be checked, the term dipia^'^'-j^ipF^i, is hermitian only if the EDM-coupling d 
is real. However as explained in Section |A.7.5| , an effective Lagrangian need not be hermitian; 
its anti-hermitian part simply parameterizes an absorptive part T — T^, where T = [S — 
is a T-matrix defined from the S'-matrix. Thus we assume (i's are complex. 



4.2 Predictions of EDMs in various models 

In all models with CP violation, EDMs of top quark may not be zero. For example in the 
Standard Model, tt7-EDM is estimated to be ~ 10~^°ecm []69|,|70|. Those for tiZ- and tig- 



EDM are obtained by changing the unit to gz cm and gs cm, respectively]^. Compared with 
the "(9(l)-coupling" result [Eq. ( |4.4| )] e/nit ~ lO^^^ecm, it is very small. This is because 
in the SM, the lowest contribution arises from three-loop diagrams, which is proportional to 



Since the operator ipia^'^'-^^ip responsible for EDM is chirality hipping, one need to pick up 
a mass of a fermion at least once. Large mass may be advantageous to obtain large EDM, but 
it may also be a suppression factor when the relevant fermion is in a loop. Thus in various 
models, EDM of a fermion is typically proportional to its mass. This makes the EDM of top 
quark very interesting. 

There are several models that induce EDM at 1-loop. These include models with extra- 



Higgs multiples and supersymmetry, which are natural candidates for beyond the SM. See |67 
for reviews on early studies. Here we summarize studies made after these reviews^. The 
conclusion of this section is as follows: the size of EDMs are expected to be 

tig-YiDU ~ lO^^^-lO^^V cm , 

and the same for other EDMs with appropriate changes of the unit , if CP- violating parameters 
in the models are 0{1). 

In a model with two Higgs doublets, it is estimated by W. Bernreuther, T. Schroder and 
T. N. Pham (1992) |7T1 that 



tt7-EDM ~ ^ = 3.1 X lO^^^ecm { 



100 GeV 



^ These can be obtained by assuming that the resuh for duui in ||7C| can appHed also for dfi^ etc. 
3 See also [Esl. 
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where the extra mf comes from the Higgs-fermion couphngs. Numerical calculations show the 
real part and the imaginary part of the Electro Weak-EDMs are the similar order each other. 
For tt7-EDM it is ~ 2 x IQ-^^e cm for = 150 GeV, mn = 200 GeV, y/s = 400 GeV; and 
ttZ-EDM ~ 0.2 X ti'y-EDM. They assumed CP-violating parameters in the model to be of 
order unity. EDMs induced by neutral Higgs boson exchange are considered also in A. Soni 
and R. M. Xu (1992) They estimated 

tt7-EDM ^ 4 X 10"^°e cm , 

ttg-EBM ~ 7^ X tt7-EDM ~ 6 x 10"^°5(,cm , 
Qte 

for niH = '2mt and rrit = 125 GeV. CP- violating parameters in the model are assumed to 
be of order unity. They argued that the effects due to charged-Higgs-boson exchange may 
be suppressed by a factor m^/M^, where is a linear combination of mf and m|^. Thus 
the charged Higgs boson contribution is likely to be subleading to the neutral Higgs boson 
exchange. They also calculated the dependence of the EDM form factors, which turn out 
to be moderate. The tt'j- and ttZ-EDM induced by Higgs exchange are calculated also in 
II, 11,11 • 



In the MSSM, one can introduce (unremovable) complex couplings with only one genera- 



tion |7|]. Induced EDMs in MSSM is considered in [|73[. They find that both tt-f- and ttZ-EBM 
can reach 10~^°ecm. Exchanged particles are gluino, charginos, neutralinos and squarks. In 
most cases gluino g gives the leading contribution. However when > 500 GeV, they find 
that chargino contribution can be larger than that from gluino due to the large Yukawa coupling 
of top quark. They assume the unification of gauge couplings and gaugino masses but others. 
In particular, they do not assume unification of the scalar mass parameters and the trilinear 
scalar coupling parameter Ag of the different generations. Their reference set of parameters 
are such that squark are lighter than 400 GeV, and complex phases in squark mixing matrices 
and in higgsino mixing parameter /i are of order unity. They vary one parameter while others 
are fixed to the reference values. The ti'j- and tiZ-EDMs in Supersymmetric Standard Model 
(SSM) was studied also in lTl|,|ri9l. The tbW-EBM in SSM was studied in pTl , [TTq] . 

There are also other models. For example in the model |^ with SU(3)c x SU(3)l x U(1)x, 



contribution of charged and neutral Higgs exchange at one-loop is estimated to 10~^^ecm and 
10~^^gs cm for the values of relative phases of the vev's such that CP violation is maximal. See 
also ||77H80[| for works related to top-quark EDM, |8ll] for EDMs of 6-quark etc. in supersym- 



metric models, and p2[ for EDMs of 6-quark etc. in leptoquark model. 



4.3 CP violating observables 

At the Lagrangian level, CP-symmetry is violated if the field-dependent phases ricrjp in Eq. ( |A.44| ) 
cannot be tuned so that the all term in C transform with the same phaseQ. Since kinetic terms 
are invariant, or "even", under any combinations of V, C, and T, a CP-odd interaction is syn- 
onymous with CP-violating one. We are concerned with CP-odd observables, which needs the 
interference of CP-even and CP-odd amplitudes: A^cp-cvcn-^cp-odd- Transformation properties 
of or £ are discussed in Section |A.7| . Here we apply it to |A1p, or da. As we shall see 
in Section [A.7.5| , observables are classified into CPT-even and -odd ones. The latter can be 



non-zero if and only if there is a finite absorptive part, which is defined as the anti-hermitian 



One usually redefines the fields to include {rjcrip) ^1"^ . See Section A. 7. 6 
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E-E E + E p-p 


p + p pxp 
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s + s 
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- + + 
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+ 
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- + - 
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+ 


+ 
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+ + - 
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+ 
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+ 
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- + - 
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+ 




+ 




+ 


t 


+ + - 
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+ 






+ 


+ 



Table 4.2: Transformation properties of the energies, the momenta 
and the spins of particle-antiparticle system. In their CM frame, E — 
E = 0, E+E = 2E, p—p = 2p, p+p = 0, pxp = and p p = — 



part T — of T-matrix, S* = 1 + iT. Sources for an absorptive part are, for example, decay 
width and final- and initial-state rescattering. For the case we are considering, e^e~ ti near 
the threshold, there is plenty of Coulomb rescattering between tt. Thus we can probe both 
CPT-even and -odd observables simultaneously. 

A word for notation may be in order. P is a field-theoretic operator acting on creation- 
annihilation operators [Eq. ( |A.44| )]; P is an operation for Ai defined by 7-matrices, which is 
c-number [Eq. ( |A.46 )]; P is an eigenvalue with respect to P or P [Tables [4.1| , A.3 |; P is an 
operation for quantities such as spin s or momentum p [Eq. ( [4.5D , Table [4.2| |. However there 
may be no confusion even if these distinctions are not made, aside from Time Reversal. 



4.3.1 P, C and T for \M\^ 

Basic discrete symmetry operations in QFT are V, C and T. These operate on creation- 



annihilation operators, and are summarized in Section |A.7.1| . We can define corresponding 
operations P, C and T for (a part of) an amplitude Ai. These are represented not by operators 
of QFT, but by combinations of 7 matrices. These operations are summarized in Section |A. 7. 4 . 
We can also define related operations P, C and T for the each term of |A^P, or da. For 
definiteness, we limit ourselves to the case we are considering: e~(pg, Se) e+(pg, Sg) — * 7*/Z* —>■ 



t{pt, St) t{pt, St), followed by t 
operations are defined by 



bW~^ {t — > bW ) etc. For particle-antiparticle pair, these 



p ■■ 


■ {P,s,p,s)- 


-i-p, 


s,-p,s) 




c ■■ 


■ {P,s,p,s)- 


(p, s 


P, s) 


= {-P, s, -p, s 


T ■ ■ 


■ {P,s,p,s)- 


-i-p, 


-s, -p, -s) 




CP ■ ■ 


■ s)- 


-i-p, 


s, -p, s) 


= (P, s,p, s) 


CPT ■■ 


■ iP,s,p,s)- 




s,p, -s) 


= {-P,-s,-p, 



(4.5) 



The equalities are for the case p + p = 0. Some combinations are given in Table [4.2| , whose 
contents can also be organized as follows: 



(CP, CPT) 
(CP, CPT) 
(CP, CPT) 
(CP, CPT) 



:+,+) 



s- s , p + p , 

sxs , pxp , 

s-s , p p , 

s + s , p - p ■ 



E-E 
E + E 



(4.6) 
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Candidates of the particle-antiparticle pair are (e^,e+), {t,t), and {b,b). Tlie last two 

originate from the decays of t and i. 

Note that the operation T do not interchange initial and final states. The tilde ~ on T 
emphasis this point. This operation T is sometimes called "naive T"0. It may be difficult, if 
not impossible, to built the field theoretic operator T for T. For example, the operator with 
the same properties as T but without the anti-unitarity nature does not do the job. This is 
because the transformation for creation-annihilation operators Ta(p, s)T^ = a{—p, —s) etc. 
with an unitary operator T, do not lead to the same transformation for ip and as T, since 
ip (and ip) contains u and v, which are c-numbers. 

With the knowledge of Section |A.7.4| , it is sensible to operate P and/or C to a part of 
a process when all relevant Feynman diagrams are topologically the same. In our case, we 
can apply P etc. separately to the final or the initial current. On the other hand, one should 
operate T to the whole process. 

Symmetry property of |A^p is directly related to that of the expectation values of observ- 
ables, which we now turn to. 



4.3.2 Expectation values 

Since an expectation value is obtained by integration over the phase space, not all terms in 
|ca/Mp contribute to it. The transformation properties under Parity etc. can be used to 
determine a relevant combination of interactions. 
Consider for example, an expectation value 



{{Pe - Pe)iPt - Pt)) = 4 {Pe-Pt) = — f PeP 

C^tot J 



* d<l> 



which is P-odd with respect to both the final and the initial currents. Thus only a part of da 
contributes to the expectation value; that is, which have the same transformation property to 
the relevant observable. There are four combinations for the SM vertices: 



c.c. 
c.c. 



Here [v^v*^] etc. are (energy dependent) couplings that combine both 7 and Z contributions. 
These are defined in Eq. ( |A.33| ). Note that P of both initial and final current should be oppo- 
site. Thus only the combination 03 = [v'^f*]* [a*^a*] -|- [a'^f*]* [f^a*] contributes to the forward- 
backward asymmetry (after one sum over spins)[|. The effective couplings [v'^v^] etc. are defined 
in Section [A.5| , and the combinations 03 etc. are defined in Eq. ( [4.26|) . 



Another example is the total cross section atot- After summing over spins, momenta are 
the only vectors. Thus each term in should be composed of even number of momenta. 

From Eq. ( [A.51D , we know that the product of the Parity of the final-currents should be even. 



which means even number of momenta of final state. Thus the number of momenta of initial 
state should also be even, which means the product of the Parity eigenvalues of the initial- 
current should be even. To summarize, an interference of amplitudes M.1M2 contributes to 

^In my opinion, "naive time reversal" seems to mean "just flip initial and final states" ; just the opposite. 

^ Before summing the spin Sf of i, (sfP^) (p^-Pf.) oc cos6te contributes to "spin-direction-specified FB asym- 
metry". This observable is P-odd w.r.t. the final tt current but P-even w.r.t. the initial e~e~^ current. Thus 
the combination 04 = [v'^v*]* [v'^a*] + [a^v*]* [a'^a*] also contributes to this asymmetry. Note that non-zero 
polarization of initial e'^e~ effectively mixes the vector vertex and the axial- vector vertex of e^e^. 
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the total cross section if and only if the Parities of the final-currents in A^i and Ai2 are the 
same, and so are initial-currents. 

CP transformation interchanges the spins of a particle and its anti-particle (Sf, Sf) (s^, St) 
at the CM frame of them. Thus the interference (A^cp=±-^c*p=±) amplitudes with the 

same CP contributes in the same sign to the polarizationsQ of a particle and its anti-particle. 
On the other hand, the interference (A^cp=±-^cp=t) amplitudes with the opposite CP 

contributes in the opposite sign to the polarizations of a particle and its anti-particle. We see 
the example latter. 



4.3.3 Several CP violating observables 



From Table or Eq. (gj), we can see that {st =F St)-{p^xp^) is (CP, CPT) = (=F, ±) and 
{st ^ St)-{pt — Pe (PePt) /IPeP) IS (CP, CPT) = (=F, =f). Here the CP operation can be limited 
to the final current. Within the Standard Model, where CP-odd contributions to both tt and 
e~^e~ vertices are highly suppressed, we see that 

{stiPe^Pt)) - {SfiPe^Pt)) = , 

{stiPt - Pe (Pe-Pt) /IPeH) " {SfiPt " Pe {PePt) /IPeT)) = , 

because they are CP-odd combinations. The former is suppressed in the open top region, 
since {st + St)-{PeXPt) is CPT-odd, which means the expectation value is proportional to the 
absorptive part of the amplitude, which come from EW and QCD radiative corrections. It is 
known 0] that the polarization Pn (~ (sf(pgXpJ)) of top quark normal to the e~^e~-ti plane 
is (9(10~^) for the SM. However for the SM, the difference Pn — Pn is far smaller than the 
current experimental reach, since it is CP-odd. 

The following are some of CP-odd combinations of observables^ (at the CM frame of e+e"): 

s// - Sji (x{st- st)-p^ CPT = - 

- OC {St - -St)iPt - Pe (Pe-Pt) /IPeH CPT = - 

sn - sn oc (st - st)-(peXpJ CPT = -t- 

Ai. = Ej-Ej if = t,b) CPT = - 

M = PeiPf^Pf) CPT = + (4.7) 

4=PeiPf + Pf) CPT = - 

Tis = 2ipj - pj), {PfXpj)^ CPT = + 

QL = 2(P/ + Pf)3 {Pf - Pf)3 - U\Pf\' - \Pff) CPT = - 

Note that p^ = {p^ — Pe)/2 and p^ = {p^ — Pt)/'^- Also, kinematical cuts should respect 
CP-invariance. 

Expectation values of CP-odd observables are not the only CP-violating observables. Event 
rate asymmetries also do the job. However, these two kinds of observables may be related 
closely. For example, Forward-Backward asymmetry A-pB of top quark is similar to the corre- 
lation (PePt) oc {cos9te)- Morc precisely, Ape = (sgn(pg-pj) /o"tot- Several kinds of event rate 



^See Section |4.5| for polarization. 

^ We follow the notations of ||2^, pT5| , [T22| 



4.4 Studies of EDMs at open top region 



69 



asymmetries are considered in literatures, 
for tLtL and tntR PHHg 



One is the difference of production cross sections 



CPT = - 



where ii = t^. These are chirahty-fiipping, and thus are small for relativistic top quarks. 
Note that from Eq. ( [4.5|) , o-itLth) ^ CFiiRtRj under CPT, while the other two chirality-non- 
flipping's, aitiiR) and aitRtLj, are self-conjugate under CP. As we shall see in Section ^.5.3| , 
the charged leptons produced in t (or more directly, W^) decays are emitted preferentially 
in the direction of top-quark spin. Thus ^+ produced in ti decay is less energetic than that 
in Ir decay. Thus the event asymmetry above is similar to A^^ defined above, which is also 
CPT-odd. Another event rate asymmetry is up-down asymmetry Aup_down 
rate difference between the events with 
below the e^e~-tt plane: 



which is the 

* above the reaction plane and the events with £^ 



A 



up-down 



(/) 



iVup(/^) 



A'^downl/^y 



iVup(/±) 



+ Ardown(/±) 



CPT = + , 



where / = ^"'",6, and N^pif"^) = N^ipif) + N-ap{f) etc. This quantity is "similar" to the corre- 
lation {{p£+ + P £-) ■ {Pe^ Pt)) 1 which is CPT-even. Yet another one is the particle-antiparticle 
asymmetry [p^ , |117|| : 



Nif) - N{f) 



CPT 



where / = i~^,b. Since both momentum and spin are summed]^ CP-odd means CPT-odd. 
Still another one is |^ 



AMf) + AMI) oc {N^if) - NM) + imf) - NbU)) 



CPT = + 



where / = £'^, b. This quantity is (CP, CPT) = (— , +), because p ^ — p for CP and p ^ +p 
for CPT. 

For pp collisions, there may be non-zero CP-odd asymmetry that is not related to CP- 
violating interactions, since the initial state is not a CP eigenstate. However such an effect is 
estimated to be small |[7B|, |103 . 



4.4 Studies of EDMs at open top region 

There has been many studies on EDMs of top quark. These exploit "open" t's; that is, 
relativistic t and t: /5 ~ 1. Although our concern is at the threshold region (3 < 0.1, we 
summarize results obtained from the analyses in open top region, for comparison. 

4.4.1 Present bounds for EDMs 

There are several EDMs for top-quark; ti'y, ttZ, ttg, and tbW. Experimental bounds on these 
EDMs at present are summarized here. Conclusion is that "(9(1)" couplings are allowed: 

ttg-EDM < lO'^^gs cm etc. 

^ References in this section are not meant to be extensive but only examples. 
In this case, both CP and CPT just interchange particle and antiparticle. 



70 



Chapter 4 Anomalous Electric Dipole Moments of Top Quark 



P. Haberl, O. Nachtmann and A. Wilch (1996) Q studied ttie present bound on ttg-EDM 
and -MDM from the pp — ^ tiX total production rate at Tevatron ^/s = 1.8 TeV. They found 
the coupling of order unity is allowed, which means ttg-EDM = IO^^^Qs cm is allowed, and so 
is for MDM. 

K. Cheung and D. Silverman (1997) also studied the present bounds on qqg-EDM and 
-MDM, and are both calculated to be 0.89 x lO^^^ggCm. They used the transverse momentum 
distribution of prompt photon produced in qg — > 
1.8 TeV, 0.1 fb"^ (Run 1). 

EW-EDMs are constrained indirectly from B - 
allowed. 

A. De Riijula, M. B. Gavela, O. Pene and F. J. 
of quarks from neutron 7-EDM: |tt7-EDM| < 2.0 x lO'^^ecm, \tiZ-EDM\ < 8.7 x lO'^^^^^cm, 
and |tt5f-EDM| < 5.8 x lO^'^^g^ cm. They used = 100 GeV; this may change the limit above. 
As they themselves pointed out, it is doubtful whether their calculation is applicable to top 
quark. Thus we don't take these limits seriously. 

In [^, CP conserving anomalous EW couplings of top quark is constrained from LEP/SLC 
data. 

181. 



Ql, 99 Ql: ^ 91 at Tevatron ^/s ~ 
[Q; it shows that (9(() 1) coupling is 



Vegas 



1991) H studied bounds on EDMs 

-18, 



For r, present bound from LEP is obtained in 



4.4.2 Expected future bounds on EDMs (Open top region) 

There are already several studies on the top quark EDMs. Some are for e+e~ colliders, others 
are for hadron colliders. We choose two representative studies for each of these, and report the 
results of their analyses. In e~^e~ collisions, ttj- and tiZ-EDM can be studied in the production 
process, while in hadron collisions, ttg-EDM can be studied. The summary of this section is, 

• Bounds of lO'^'^ecm on tt'-f- and tiZ-EDM are possible at a e+e~ colliders at ^/s = 
500 GeV. With the use of "optimal observable", this can be improved by factor 2. How- 
ever this observable may have larger systematic error than simpler observables. 

• Bound of 10~^^-10~^^ gs cm on tig-EDM is possible for hadron collisions at ^/s = 14 TeV. 



(i) Anomalous ElectroWeak vertices in e+e or 77 collisions 



i-0) Detector simulation ||1 13|1 was performed with JC = 50 fb ^, nit = 180 GeV, ^/s 



500 GeV, polarization of e~ beam = ±80%. They used the decay chain tt bbWW bbqq'iiy, 
where £ = e, /i. The relevant branching fraction is 8/27. They obtained overall efficiency of 
the analysis, including branching fractions, reconstruction efficiency, and acceptance, ~ 18%. 
Their resuhs is that the sensitivity for EW-EDM is 0(0. 1) with 90% CL. 

(i-1) D. Atwood and A. Soni (1992) They studied the sensitivities for tij- and ttZ- 
EDM, (and -MDM) of e^e^ —* ti with ^/s = 0.5 TeV. No specific model is assumed. Assuming 
that the top decay vertex tbW is the SM one, they claimed Icr-bounds are ~ 10~^^ecm for 
lO'* events. They introduce the notion of an "optimal observable" that have the maximum 
sensitivity to the relevant coupling. This observable turn out to be powerful as we shall see. 
However in general, this observable is rather complicated combination of momenta (and spins, 
maybe). Thus as themselves noted ||103|] , there may be larger systematic error than that for 



simple observables. The events they use are that with both W decay leptonically: bi^ugbi vg. 
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Re[tt7-EDM] Im [tt7-EDM] Re [ttZ-EDM] Im [ttZ-EDM] 



4 X 10' 


47 


1 X 10" 


-17 


1 X 10' 


-17 


4 X 10' 


-17 


4 X 10" 


48 


1 X 10" 


-18 


4 X 10- 


-18 


4 X 10" 


-18 


4 X 10" 


-19 


4 X 10" 


-19 


4 X 10" 


-19 


4 X 10" 


-19 



without Ew, with optimal 
with Ew, with simple 
with Ew, with optimal 



Table 4.3: D. Atwood and A. Soni (1992) Q. la sensitivities in 
units of ecm. "With -Evk" means "to use polarization of W" . Sim- 
ple observables are explained in the text. See the original paper for 
optimal observables. 



In this case, one can measure the polarization E!^ of W, statistically by the distribution of 



E\ 



where its "gauge" is fixed by an arbitrary lightlike vector ujq. As summarized in Table O , 
the use of E!^ improves the sensitivity about (9(10^). They propose several observables that 
are expectation values of certain operators. Explicit formulas of the operators for optimal 
observables are quite lengthy. They also studied simple operators. Those with best sensitivities 
are 

e^.p.P.^Q'^H+PH-^ for Re[tt7-EDM] , 

e^.^p^P^Q'^.H+PH^" for Re[ttZ-EDM] , 

H--Q^ for both Im[tt7-EDM] and Im[ttZ-EDM] , 

where P^, = -p^, = -p^-, = p^+ +p^~, = 2 {Ew+-Pt) Ew+ ±2 {Ew~-Pt) Ew-- 

(i-2) M. S. Baek, S. Y. Choi and C. S. Kim (1997) 13 . They studied the sensitivities 
for tt'y-EDM and ttZ-EDM of e~^e~ collisions and 77 collisions at NLC with ^/s = 0.5 TeV, 
J Cee = 10 fb~^ for polarized electrons and the twice for unpolarized electrons. No specific 
model is assumed. The real part of EDMs are bounded by CPT-even observables (when there 
is no absorptive part). With polarized electrons, the tightest bound is obtained through A\, 
defined in Eq. ( |4.7|) , in the inclusive top-quark decay mode: 

I Re(tt7-EDM)| < 1.4 x lO'^^ecm , | Re(ttZ-EDM)| < 2.3 x lO'^^ecm . 

Likewise, the imaginary part is bounded by CPT-odd observables. With polarized electrons, 
the tightest bound is obtained through A^, defined in Eq. ( |4.7|) , in the inclusive top-quark 
decay mode: 

I Im(tt7-EDM)| < 1.8 x lO^^^ecm , | Im(ttZ-EDM)| < 3.1 x lO^^^ecm . 

For polarized 77 colhsions, |Re(tt7-EDM)| < 1.8 x 10~^''ecm for y/s = 0.5 TeV, and 0.2 x 
10~^''ecm for 1.0 TeV. They argued that the sensitivity is much sensitive to ^/s than e~^e~ 
collisions. 

(i-3) These couplings are also studied in |9|^. In l9|,|9|], ti-f-EDM and tiZ -EDM are 



consider within two Higgs doublets model. QCD corrections to this process well above the 



threshold are considered in Studies in 77 colliders are given in p9| , |100|| ; sensitivity is 

~ 10-i^e cm with JC = 10 fb"^ and ^/s = 500 GeV. 
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(ii) Anomalous Chromomagnetic vertex in hadron collisions 



(ii-1) D. Atwood, A. Aeppli and A. Soni (1992) ||103|| . It was studied that the sensitivities 
of gg — > ti at SSC and LHC for the real part of tig-EDM. la sensitivity is 0.3 x 10~^^ gs cm 
with optimal observable /opt, 0.5 x lO"^^^'^ cm with /i, and 0.6 x lO"^^^^^ cm with /2. We do 
not reproduce their explicit form. Assumed experimental setup is mainly, SSC, y/s = 40 TeV, 
10fb~^, 10^ leptonic ti events. However, they say that the difference is only a few percent 
by changing SSC y/s = 40 TeV to LHC i/i = 16 TeV. Correlations between the t and i 
polarizations are used: /opt, /i, /2, /s- The /opt and /i are tailored to maximize the sensitivity 
for Re(tt(y'-EDM), and are fairly complicated combination of momenta and spins. As they 
pointed out, there may be larger systematic error than simpler observables such as /2 and /s: 

J'. 



{PerPe2PbrPb2, 



1/2 



/a turn out to be less sensitive than /2 by C(10). 

(ii-2) S. Y. Choi, C. S. Kim and Jake Lee (1997) |lOl. The sensitivity of gg ti (LHC) 



for ttgf-EDM was studied. Assumed setup is acceptance efficiency e = 10%, Bp = Bi = 20% 
for i = e, fi, y/s = 14 TeV, and JCpp = 10 fb^^. CP-odd energy and angular correlations are 
used for unpolarized proton beam, la sensitivity is 

I Re(%-EDM) I = 0.899 x 10'^^ gs cm • • • T^^ , 

I lm(tt^-EDM)| = 0.858 x lO^^^^fe cm ■ ■ ■ , 

e 

33 ' 



lm(tt^-EDM)| = 0.205 x lO^^^^f, cm ■ ■ ■ Q 



where written in the RHS are relevant observables defined in Eq. ( |4.7] ). There may be (at 
least) two more operators that are CP-odd: A{ and Ag, also defined in Eq. ( [4.7|) . However 
these correlation are zero, since the cross section is symmetric under {pg, pg) {pg, Pg) due to 
Bose symmetry: \Pg,Pg) = o){Pg)o){pg) |0) = + \pg,PgY For polarized proton beam, one can 
probe CP-violating effect without detailed reconstruction of the final states. Thus they used 
CP-odd rate asymmetry, 



a I + a_ 



where a± is the cross section for tt production in collision of an unpolarized proton to a polarized 
proton of helicity ±. They obtained | Im(tt5f-EDM) | < 10^^'^(?sCm. They found sensitivity 



depends crucially on the degree of polarization of gluons. Berger-Qiu parameterization [|105 
for polarized gluon distribution functions was used. 



(ii-3) In |7|, N{tLtL) - N{tRtB) was studied. In ]T06| Tevatron Run II (v^ = 2TeV) was 



considered. In 77 — * tt is also considered; a list of sensitivities is also given. In ||108|| , 



two-doublets Higgs models was considered. In ||107|| the process gW^ th was studied, and 
the contribution turned to be small; the reason is also given. 

(ii-4) Here are comments for the uncertainties in the theoretical analyses of sensitivity stud- 
ies at hadron colliders. Some of the existing sensitivity studies make use of some complicated 
combinations of momenta of leptons and b's; those include "optimal observable". Whether 
these analyses can be taken seriously was questioned; there are several reasons. Firstly the CM 
energy of a reaction cannot be measured in hadron collisions. In view of event reconstruction, 
more questionable is the efficiencies for the assignments of jets to the parent partons. A jet. 
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which is assigned to a parent parton, must be well separated from other jets in order for it 
to be used as a signal jet; in other words, not all hadrons are assigned to one of the signal 
jets, which means the jets used for the analysis. This is in contrast to e+e- collisions, where 
every hadrons are assigned to one of the jets. Generally certain amount of energy is not cor- 
rectly assigned but included into other (non-signal) jets. Usually the missing energy will be 
supplemented artificially by multiplying some factor to the observed energy. Under such cir- 
cumstances it should be examined carefully, incorporating realistic detector simulations, how 
well the required parton four-momenta can be reconstructed. 



(iii) Other situations 



Anomalous chromo-EDM vertex in e"^e~ collisions is studied in ||101| , |102|] . In ||101|| , gluon jet 
energy distribution for the process e+e~ tig is studied. They concluded that the sensitivity 
is (9(0 1) with ^/s = 500 GeV, JC = 50fb~^, identification-efficiency for top-quark pair- 
production event ~ 100%, and E^^^ = 25 GeV. The last one is the cut for the minimum 
gluon-jet energy, which is necessary because (1) the cross section is IR singular, and (2) and 
to avoid the contamination due to the gluon radiation from 6-quark; the cross section grows 
rapidly as ii^™™ is lowered. On the other hand, in [|102|| , CP-odd spin correlation of t and t is 
studied. They concluded that the sensitivity is 0(10) with ^/s = 500 GeV, JC = 50fb~^, and 
top-detection efficiency e = 0.1. 



Anomalous EW vertices in hadron collisions are studied in [[76| 



Anomalous decay vertex is studied in [|109H112| . Both production and decay anomalous 
vertices are studied simultaneously in ||113| - |119|1 for e'^e' collisions, and in [p.20|| for hadron 
collisions. CP violation in top sector is studied also in ||121|| . CP violation in other sectors 
(EDM of T etc.) are studied in p2HT2l. 



4.5 Polarizations 

In the previous section, we saw that the spin-direction of t and i may serve as a probe of 
EDMs. In this section, we see how the effects of polarization can be incorporated to actual 
calculations, and how the polarizations can be measured. 



4.5.1 Polarizations of e+e 

Polarization is a incoherent sum of states with definite spin direction. Thus this can be treated 
by a density matrix. For spin- 1/2 particle, polarization P is defined to twice the expectation 
value of spin S: P = 2(5'). Thus the density matrix p is 



2(i + Ep^-0 = 2U + ^P. 1-p/J 



because tr (p) = 1, tr (pSj) = {Si). Then for longitudinal polarization, we just sum "spin-up" 
and "spin-down" with the weight and respectively. Since electron mass is negligible, 
chirality is the same to helicity for electron, and the opposite for positron: {iPlY = ('^'^)_r- Thus 
the matrix element E^'^ ^ '^spms^sejge initial e+e~ current with incomplete polarizations 
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X 

1 I — 1 ^ 




0.4 / ,/■■ / / /l-j 
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+ _ / ■' : /' ■'' 

^ - ; / / r J ; 4 

Oh ; / /■ : / / \ 
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P(e") 

Figure 4.1: Contours of tiie measure of initial polarization: x = 
(Pe - Pe)/(1 - PePe). In the figure P(e-) = Pe and P(e+) = P^. 
From left-top to right-bottom, the contours are for x = 0.99, 0.9, 0.8, 
0.4, 0.0, —0.4, —0.8, —0.9, and —0.99. For unpolarized positron beam 
Pe = 0, X = -Pe- 



can be written as follows 

E 



Here Pe = 1 (Pe 
(positron), and 



fJ-U 



tr A,, 



X A,, 



1+ Pel + 75 ^ 1 



Pe 1 - 75 



l+Pel+75 ^ 1 



Pel 



75 



1-P.P, 



2 2 2 
tr [A^(l - X75¥e^.;^e] • 



X 



1) means the complete polarization with right-handed helicity for electron 



X 



1 - PePe 



(4.8) 



is a measure of initial polarization; x = means (Pe,Pe) = (=Fl,±l)- The contours of x 
are shown in Figure |4.1| . We can see that non-zero polarization effectively mixes vector- and 
axial-vector-vertices, or x 7^ means that the initial state is parity violating. Thus when 
X 7^ 0, an expectation value of parity-odd observable can be non-zero, even if all the relevant 
interactions are parity conserving. 

The effect of non-zero x can easily be recovered from the results with x = 0, in the following 
way. Since we don't assume anomalous vertices for initial e~^e~ current. 



tr[(t;^^7M - a^''w)a ' X^Mi^'"^*!^ ' «^"'*7.75)#e] 
^ tr[7;./^7.;^^ {{v'^'v'^* + a^^a^^*) + xiv'^'a^''* + a^^t;^^*)} 



+ 7MK7.#e75 {(t;^^a^^* + a^^t;^^*) + xiv^'^'v^^* + a^^a*^^*)}] . 

Thus we can recover terms proportional to x from those without x follows: every terms in 
|A^P contain two factors of or a*^; exchange each one of them with the rule v'^ ^ a*^ and 
divide by two. The product of couplings in Eqs. ( [4.26| ) and ( [4.27|) are defined such that each 
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75 



pair of (ai, 02), (03, 04) and (05, Og) are "conjugate" with the procedure above. Thus we can 
recover x dependence simply by replacing (ai, 02) with (ai + ^02, 02 + xai) and similarly for 
(03, 04) and (05, oe). 



4.5.2 Polarizations of tt 

Although the spin s of, say, t is surely a 3-vector, but one should be careful to its quantum 
nature. That is, two spins that are orthogonal in a sense of 3-vector, are not "orthogonal". 
For a spin- 1/2 state, 

\ia Aini2 2 ^ 1 + cos^ 

\{eA\z)\ =cos - = — - — , 

where 9 and (j) denote polar coordinates, and \z) = = 0). Thus, when an amplitude vanishes 
for some choice of spin s = —n, it means the spin-direction is +n, not perpendicular to — n. 

Let us look at them more closer. In quantum mechanics, the polarization P along some 
direction n = 1) may be calculated as follows: 



p=±i 

Here we used 

P 

1 = ^ |s = Pn) (s = Pn| , n-S \Pn) = - \Pn) . 
p=±i 

Note that Si means spin operator, while Sj means spin-direction: |s|^ = 1. 

With this experience in mind, next we consider the process e^e' — > t(s)t(s), where s (s) is 
a 3-vector for the spin-direction of t (i) at its rest frame. Let us write the spin dependence of 
the matrix element M. as follows: 



4 



+ -s s^M ^ 
4 



(4.9) 



As the notation indicates, Pj can be interpreted as the polarization vector for t. This can be 
seen as follows. Total cross section a is 



croc Jdn J2 \M{s\s^)\'^^ Jdnc 



Note that momentum 3- vectors are summed over all direction, while spin 3- vectors are summed 
only "up" and "down", because this is the way they form a complete set. The directions n,n 
of quantization-axes are arbitrary, here. While, the polarization vector P* of t for a fixed 
momentum configuration is 

{n'-2S') oc P\M{s' = Pn\ = Pn^)f = n'P' C . 

p=±i 
p=±i 
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The spin quantization-axis n for t is arbitrary. Dividing by a, we see that P* defined at Eq. ( |4.9| 
is surely a polarization vector for t, when one sum over the spin of i. Likewise, 



p=±i 
p=±i 

Thus we see 

which means the spins of t and t are correlated when M*-' ^ 0. This can be seen already in the 
definition Eq. (^) ; when M*-'' = 0, \M{s' = -F\s^)\'^ = irrespective of (when |Pp = 1). 
This means the spin-direction s is P, and is determined solely by the momentum configuration, 
because so is P. In general M*-' 7^ 0, because momenta and spins have independent degrees 
of freedom. However in some cases, the two spin-directions factorize. For example, for the 
process e~(L) e~^{R) — 7* — t{s) i{s) near the threshold. 



where f3t = '^\Pt\/^/s■ Here e~^{R) means positron of right-handed helicity, not chirality. This 
shows that the direction of spins of both t and t are anti-parallel to the direction of (left- 
handed) initial electron. This is due to the conservation of total angular momentum; on the 
threshold, there is no orbital angular momentum. 

The polarization vectors P,P should be composed of odd number of and Pg, since | p is a 
scalar under 3-dimensional rotation. This means the polarization vectors are odd under "over- 
all" parity transformation [Section |A.7.4|| , which means parity violating. These vectors are 
non-zero even for zero initial polarization x, because the coupling of Z to fermion-antifermion 
is parity violating. 

A |A4p for the definite momenta and spin-directions of t and t can be calculated with the 
usual Trace technique with spin-projection operators [Section [A.3.7 |. One can easily see that 



the general spin structure for |A^p is that given in Eq. (|4.9| ). In the next section, it is given 
that the simple example of a calculation with a spin-projection. 



4.5.3 Spin direction of top quark 

The direction of the spin of top quark can be measured statistically for the decay process 
t —>■ biui. This is because [Q at the rest frame of t, the charged lepton i is emitted preferentially 
to the direction of the spin of top quark. This can be seen as follows. Using MiFf 2 = — t'gF'^'u^ = 
— ■U2F'^f 1 for commuting spinors, and Fierz identity, 

= |M(pb)(l + -f5)v{pu) u{Pe)i^ - 75)u{pt) , 
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which means 



= ^Pi ■ {pt -mst) 

= ArritEiil + Sffii) at the rest frame of top. 



That is, roughly speaking, St // + and si // — {p^-). Thus the component n-St of the 

spin St of top quark parallel to certain direction n can be obtained statistically by (n-pi). Then 
we treat the spin (or polarization) of top quark as if it were a observable. 



4.6 Corrections due to Coulomb rescattering 
4.6.1 Relevant Lagrangian 

From now on, we concentrate to the ti production near the threshold. This also means to use 
e~^e~ collisions, since for hadron colliders ^/s of each subprocess is not fixed. As emphasized in 
Chapter the most prominent feature of threshold region is leading-order multiple Coulomb 
rescattering[^ , which means higher orders of are not suppressed compared to the propagation 
without gluon exchange. This drastically changes the behavior of ti production cross section, 
compared to the tree-level expectation. For example, the magnitude of the "imaginary part" 
(= absorptive part) of the amplitude is similar to the "real part". 
The relevant part of the Lagrangian is as follows: 



Zmt 

+ terms required by gauge invariance + ■ ■ ■ , 



where the couplings are defined in Eq. ([A.32|) . Feynman rule for, say, tt'j vertex is proportional 
to QtJ^ + dt^yl {2mt)a'^^''i^qy, where is a momentum flows into ti current. Note that 



Zm 2m 2m 



where the first "~" means they coincide with the use of eqs. of motion, and the second means 
to sandwich with {pt,Pt\ and |g) and to integrate over x [Section |A.4|| . We checked that both 



of these forms of the interaction give the same result, at least for the case we deal with. 



Since the Coulomb rescattering is indispensable for the threshold region, in some case its effect is implicit, 
and the word "rescattering" is saved for the other corrections; that is, gluon exchange between t and b for 
example. However here we do not care about those corrections involving decay products. 
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As we saw in Section p73| , tlie potential Vq between tt due to C^t = —gsitl^TH)G'^^ is 

Vc = —CFgs/\q\'^i where q is momentum transfer from i to t. Likewise, one can show that the 
potential Vcedm due to the itg-EDM interaction above is [Figure 



This potential is (P, C) 



Vcedm 



-CfQs -d. 



tg 



+ 



as it should be. 



iq- 



2 



2 



(4.10) 



t 



t 



Figure 4.2: Contribution of anomalous ttg-EDM to tt potential. 

Note that the magnitude of the EDM-interactions are proportional to the momentum trans- 
fer between ti. This can clearly be seen in the modified vertices Eq. ( |4.12 ) in the next section. 
Thus the effects of two EW-EDMs becomes larger for higher CM energy ^/s, because they 
affect directly the tt production vertex. While the effect of CEDM do not become larger for 
higher y/s, because|3 it affects the Coulomb rescattering between ti, whose relative signifi- 
cance becomes smaller for y/s. Thus the threshold region does not suit for the measurement 
of EW-EDMs. On the other hand, CEDM can be measured only near the threshold for lepton 
colliders. Note that even though the sensitivity to the existence of finite CEDM may be better 
for hadron colliders, the precise measurement of its magnitude may only be achieved at lepton 
colliders. 



4.6.2 Corrections to vertices 

In the next section, we show that the effect of Coulomb rescattering between t and t is to 
change the vertices for ti'y and tiZ. The result is as follows. To the tree-level (and without 
anomalous vertices), the electroweak ti vertices are 



v'^V^ - a*^r^ , where F^;, = 7^ 



J- A 



times —igx- See Section |A.5| for the definition of the couplings. Rescattering and anomalous 
vertices modify the vertices: 



where 



) -V iy-a + dtxi- EW-EDM ' 



(4.11) 



^ V 



J- A 



YG{E,\p\)-i^,^dtgD{E,\p\) 
mt 

2 



\P\ 

nit 



Yl5F{E,\p\) X 



nit 



— 

Tew-edm = ^75— — i^(^, \p\) X 
nit 



{E + tTt, 
\P\' 



{E + tVt 



(4.12) 



More concretely speaking, because the combination of the Green function D = F — G, which is relevant 
to CEDM, becomes zero for y/s ^ 2mt. 
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where p = {Pt — Pt)/'^ is the relative momentum of tt, and E = y/s — 2mt is non-relativistic 
CM energy; G{E, \p\) and E{E, \p\) are S- and P-wave Green functions, respectively 



G{E, \p\) = {p\G\x' = 0) , pF{E, \p\) = {p\pF\x' = 0) 



and 



and D{E, \p\ 
that 



UJ + Vn]G 



\P\ 



uj + Vc] pF = p 



(4.13) 



(4.14) 



rrit / \ nT-t 

G{E, \p\) — F{E, \p\). We can see from the equations for Green functions 



G{E,\p\) = F{E,\p\) = l/ 



\P\ 



nit 



{E + iVt 



for 1^ , 



and D{E, \p\) = 0. Thus in the weak potential limit, or for ^/s ^ 2mj where the effect of 
Coulomb rescattering is diminished, the vertices reduces the tree-level ones with anomalous 
EW-EDMs but QCD-EDM. 

We can see from Eq. ( [4.12|) that the effects of anomalous CP-violating interactions are 
suppressed]^ only by Pt = \Pt\/^t- Thus we take into account 0{/3t) corrections, and neglect 
the higher orders^. 



4.6.3 Sketch of the calculation 

Using Non-Relativistic form of the propagators, 

^^ 2^ 2 E/2 + kO -\q\y{2mt) + irt/2 



9- 



7 



+ o{- 



+ 



2' 2 E/2-k^ -\q\'^/{2mt)+iTt/2 \c 
Lippmann-Schwinger equation for the vector vertex is written as follows: 



ry{E,p)=Y + Gp{-tg, 



d^'k 



(27r)4 P/2 + F - |fc|V(2mt) + iTt/2 
i 



X 



E/2 - A;0 - |fc|2/(2mi) + iTt/2 





X 



A" 



+ 

+ 



d 



tg 



2mt 



C 



2m 



! c 



3 



\p — k\ 



(4.15) 



(4.16) 



13 



Note that the time-component of the vertex do not contribute to |A^P, since the initial e+e current 



is conserved q^j'^ 



(because the mass of e can be neglected). Thus at the CM- frame — (^/s, 0), the 



time-component of e+e current is zero j 



0. Propagators of gauge bosons are proportional to g^^, aside 



from g^gi/. Thus only the space-component of tt current contributes. 
We also neglect the 0{(3t) corrections from FSIs. 
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A similar equation holds also for the axial vertex. Since there is no p^-dependence on the 
right-hand side, consistency requires ry{E,p) = Ty{E,p). Thus we can trivially integrate 
and obtain 



ri (F «1 = V - /^'^ -'^f9s I 

7 J(2n)^\p-k\^\k\ymt-{E + irt) 



X 



X 



1 + 7' 



(4.17) 
1_/ 



{E, k) - ^iP-J^ {a^%rV(^, k) + T^E, fc)a^%5} 
zrrit c 

In view of the structure of gamma matrix, we decompose the vertex function ry{E,p) as 
follows: 

jU6'^TG{E,p) + rj^{E,p] 



+ ia'^^,ri(E,p)-i^,r^(E,p) 



1-7^ 



(4.18) 



where "^ij^^iE, p) = 0. By plugging this expression into the integral equation above, one 
obtains the integral equations for each (reduced) vertex functions rG{E, p) etc. One can see 
that r^(S,p) = 0{dtg), T%{E,p) = 0{d^l). Thus we forget about T%{E,p) hereafter. One 
can do just the same procedure for the axial vertex r^(£', p), and finds that anomalous ttg 
vertex we inserted do not alter the axial vertex to the order we consider now. We extract 
Lorentz structure of T'^p{E,p) etc. as follows: 



Tg{E, p) = f ^ -{E + zT,)) G{E, \p\) , 



rUE,p) = 



-{E + iTt 



—p 

rritc 



F{E,\p\) , 



P^(£;, p) = f ^ -{E + iV^)] -^dt,D(E, \p\) . (4.19) 
In terms of G{E, \p\), the Lippmann-Schwinger eqs. become NR Schrodinger equations: 



\P\ 

\P\ 

\pI 

rut 



{E + .r,) j G(E, IpI) + [Vc{\p - k\)G{E, \k\)\ = 1 , 



d^A; 



{E + ir,)j p'F{E,\p\) + 



(E + irt)jp'D{E,\p\) + 

fd^k 



Vc{\p-k\)k'F{E,\k\)\ ^p' 
Vc{\p-k\)k'D(E,\k\)] 



(27r)= 



Vc{\k-p\){p-kyG{E,\k\) 



We find that D = G- F. As for EW-EDMs, since 



2m, 



rritC 



we can put as 



rWEDM(^,P) = il^—K{E, \p\) - {E + iTt)\ 

mtc \mt J 



(4.20) 



(4.21) 
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By comparing witli 



7*75 ~ i(T*^75- 



rritC 



one can see tliat p'^K{E, \p\) satisfies tlie same Lippmann-Schwinger eq. as p^F{E, \p\), thus 
K{E,\p\) = FiE,\p\). 



Likewise, 



^\ ~ ^cr^^754 - -iE + tVt)^ Yl^FiE, \p\) . 



4.6.4 Polarization vector of top and antitop 

Including 0{Pt) contribution, the production cross section for ti pair can be written as follows: 
da 



where 

n \Pt\ a PePt 
Pt = , COS Ote = -. n r ' 

and X is a measure of the polarization of initial e~^e~ defined in Eq. ( |4.8|) . The coefficient Cfb 
is defined in Eq. ( [4.23| ). For unpolarized positron beam (Pe = 0), x = ~Pe, where Pg = — 1 for 
the electron beam that is completely polarized to left-handed helicity. Hard vertex corrections 
are not included, since we are working with non-renormalizable interactions. In the formula 
above we summed over spins of t and i. Before summing over the spins of tt, the production 
cross section for them can be written as follows: 



da{st,St) da 1 + P-St + P-St + {st)i{st)jQ 



d^Pt d^Pi 4 
As was explained in Section [4.5.2|, P and P can be interpreted as polarization vectors for t and 



t, respectively. Both SM and anomalous interactions contribute to the polarizations: 

P = PSM + 5P , P = PsM + . 



With the knowledge of Section |4.3.2| , the interference of an anomalous vertex, which is CP-odd, 
and the standard one, which is CP-even, contributes in the opposite sign to the polarizations 
of t and t. On the other hand, the interference among the standard interactions contribute in 
the same sign: 

PsM = PsM , ^P = -5P , (4.22) 

or 



PsM = (P + P)/2, <5P = (P-P)/2 
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Note that we are working to the leading order of the anomalous couplings. Hereafter we express 
those vector by components: 

P = P//n// + PLUi. + PnUn , St = s//n// + sj_nj_ + snHn , 

where the axes are determined by the beam direction and the production plane: 



n 



Pe 



UN 



Pe^Pt 



\Pe\ iPe^Ptl 

Now the polarization of t is as follow^: 



Psm)// = C// + Re ( C/}^ ) PtcosOte 



[Psm)l = Re [ C^- ] Pt sin 9te , 



(Psm)n = Im y^j^Qj f^t sinOte 

for the contributions from the standard interactions!^ , and 

SP// = , 



6P^ 



Im ( Bfdta?; ) + Im ( B^dt,^- ) + Im ( Bfdtz- 



D 



G 



'^G 
F 



Re ( Bf,dtg- ) + Re ( B^dt^- ) + Re ( B^dtz- 



G 



G 



PtsinOte ; 
A sin Ote 



for the contribution from the anomalous interactions. Coefficients Cjj etc. are defined below 
There is also correlation between St and Sj: 

(st)i(st)jQij = {st)i{st)jQij^-^ + {st)i{st)j 6Qij , 

(st)i(st)jQygj^ = S//S// + (s//sj_ + s±s//) Re (^Cn^^ PtsinOte 
+ (s//Sn + snS//) Im (^^^^^ f^t sin 9te , 



{st)i{st)j 6Q,ij = 
+ {s//s± - S^S//) 



Im ( B'^dtg^ ] + Im f BZdt^^ ] + Im f B^,dtz^ 



Re BU 



G 
D 
G 



G 
F 
G 



Re ( B7dtj^ ) + Re ( Bfdtz 



Pt sin 9te 
(3t sin 6te . 



We use notation similar to |2J]. The are two differences: (i) as and 04 of them are twice of ours; and (ii) 
Cat of them is —Cn of ours. 



16 



These resuhs coincide with those in 



4.6 Corrections due to Coulomb rescattering 



83 



We do not use these correlations in this work. 

Coefficients of the polarization are defined as follows: 



cf,{x) = 

Ci.{x) = 



«2 + X^l 
04 + XOg 



for the standard vertices, and 



<^N(,Xj — \ — (^FB 



(4.23) 



Biix) = -1 , 5&(x) = +q?(x) , 

05 + xae 



i?r(x) 
5r(x) 



«1 + X«2 

Q6 + XQ5 
ai + 



^n(x) (^t7 + ^n(x) dtz , 



(4.24) 



where 



Blix) 



1 



^f(x) 
5^(x) 



5n(x) 



oi + Xa2 
1 

Ol + X«2 
1 

Ol + X«2 
1 

Ol + X«2 
1 

Qi + Xa2 
1 



{ {[v%']*v'^ + [aV]* o^^^) + X {[v'vf a'^ + ^t;*]* t;^^) } rf(s) 
(xbV]*^;^^ + [aV]*v"^) , 



(4.25) 



Ol + X«2 

for the anomalous vertices. Symbols ai^e are the combinations of EW couplings: 



a V 



02 = 2 Re ([^;V]* [aV]) , 



(4.26) 



for the standard vertices and the anomalous ttg^-vertex, and 

= [v%'Y v'U'^ + [a%']* a'^d'^ + [v%']* v'^d{s) d'^ + [a^^;*]* a'^d{s) d'^ 
= [^;^^;*]* + [v^']* v'^d{s) d'^ + [a'v']* a'^d{s) d'^ , 

Qg = [v^vf [a^d*] + [a%']* [v^d'] 

= [v^vf a^'^d*^ + [a^v'Y v^'^d''^ + [v^v'Y a^^d{s) d'^ + [a^v']* v^^d{s) d^^ 



(4.27) 
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{G*D) 




G*F 






(SM)* (SM) 


(SM)* (EW-EDM) 




tti (12 


CI3 (24 




Ptt 


+ (-) +(-) 






Pe-e+ 


+ 


+ 


+ 



Table 4.4: "Eigenvalues" of the products of couplings oi etc., defined 
in Eqs. ( 4.26 ) and ( 4.27 ) under Parity transformations P for the final 
{tt) and initial (e~e^) currents. See text for the details. 



for the anomalous ^7- and ttZ- vertices. Symbols [v'^a*] etc. are defined in Eq. ( |A.33| ), and d{s) 



is a ratio of Z-propagator to 7-propagator [Eq. ( |A.34| )]. These combinations ai^g are chosen so 



that the dependence on the initial polarization can be seen easily; each pair of (ai, 02), (as, 04), 
(05,05) is "conjugate" each other with respect to the exchange e^e~ couplings [Section [4.5.1|| . 

Table is "eigenvalues" of the products of couplings ai etc., defined in Eqs. ( |4.26|) and 
(|4.27|) under Parity transformations P for the final {tt) and initial {e~e^) currents. It read as 



follows. For example, the coupling = [v'^v^]* [v^d^] + [a'^f*]* [a^d^] is a coefficient of the cross 
term of SM contribution and EW-EDM contribution, A^sm-^ew-edm- Thus in |A^p, 05 is 
accompanied by G*F, since 

A^SM ~ G + F = C(l) + 0{Pt) ■ ■ ■ CP-even 

A^c-EDM --D = 0{Pt) ■ ■ ■ CP-odd 

-Mew-EDM ~ = 0(A) ■■■ CP-odd. 

Note that we are working to 0{f3t)- Under P^, f*-part of A^sm is transformed with the 
eigenvalue +1, while (i*-part of A^ew-edm is —1; thus in the term with 05 should be — 1 

under Pff. A complication occurs for oi and 02, which contain Vt, which is a coefficient of both 
^'^G, which is P-even, and i'^^{p^ /m)D, which is P-odd. The former is 0{(3^), while the others 
are 0{j3). Since now we are working to 0{(3), only the interference G*D contributes, as far 
as D is concerned. Thus can be Ptt = il? where +1 for jCp, while —1 for G*D. Thus 
denoted in parentheses for oi and 02 are for G*D. On the other hand, Vt in os-e should be SM 
contribution, not chromo-EDM contribution. This is can be seem from the order counting for 
the velocity /5; all anomalous contributions to M. are 0{(3), while SM contributions are 0{1). 
This table can be used as follows. For example in |A1p, a term with (P^, Pe-e+) = (+> +) 
should be accompanied by ailGp. Likewise, a term with (Ptf, Pe-e+) = (~) ~) should be 
a^G*F, which is SM contribution, or a2G*D, which is the interference between chromo-EDM 
contribution and SM one. 



Physical considerations 

Let us see if our results for da{st, St)/d^Pt is plausible. For this purpose, it is convenient to 
express it in terms of vectors p and s, not the components of them; see Section p.5| . We 
already saw in Eq. ( [4.22|) that CP-property of the terms in |A^p relates the s-dependence and 
s-dependence [See also Eq. ( ^^ ]. Here we shall see that Parity-property can be used to derive 
momentum-dependence. As was explained in Section |4.5.1| , non-zero polarization x ior initial 
e~^e~ current mixes the vector and axial vertices, thus it complicates symmetry considerations. 
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(P ^ PA 


(bJVlj [bM) 




(bJVlj [hiVv-tjDM.) 


\Pe? 


0{l) 


( + >+) 














Re [a3G*F] 






Pe 




(-5 +) 


I/^|2 






Pt , Pe (Pe-Pt) 




( + . -) 




aiG*D 


a5G*F 


PeXPt 


0(A) 




asG*F 


a2G*D 


aeG*F 


be). iPe)j 


0{\) 


( + >+) 


ai|G|2 






(Pe). iPt)j 








a2G*D 


aeG*F 


{Pe)i{Pe^Pt)j 


0{f3t) 


( + >-) 


a4G*F 


aiG*D 


a^G*F 



Table 4.5: Note that spin (s) is even under P. Thus P^f etc. is 
determined solely by momenta (pg,pj). On the other hand, CPT 
depends on s. 



However as explained in Section [4.5. 1| , the results for non-zero x can easily be recovered from 
those for X = by replacing couplings appropriately: ai ai + for example. Thus 
in this section, we consider the case when the spins of e^e" are summed with even weights: 
X = 0. Before calculating the traces for |A1p, it is easy to see that each momentum^ {Pt,Pe) 
can appear only twice or less in |A^P, and each spin- vector (s^, St) can appear only once or 
less. However since we are working to 0{(3t) or less, momenta of t and i cannot appear twice: 
\Pt\ = f^fytmt, 

Let us start with the unpolarized part dcTunpoi- As is shown in Eq. ( |A.53| ), CP-odd part 
of the (product of the) tt current {jj^)tt = jttjlt ^'^ '^'^^ contribute to spin-summed cross 
section. Since the anomalous interactions are CP-odd while the Standard ones are even, the 
interferences between them are CP-odd. This means those interferences (A^sm-^edm) do 
not contribute to ddunpoi; the interference among anomalous contributions (A^edm-^edm) is 
0{P^), and thus is beyond our concern here; only A^g^^A^sM contributes to ddunpoi- Since 
dcTunpoi is a scalar, (pg, pj-dependence of it is either IPgP; Pe'Pty or see Table ^75|. Among 



these, IpjP is 0{j3f), and thus is beyond our concern here. "Eigenvalues" under P^+e- etc. 
can be easily seen; for example, p^ Pt has odd number of p^ and odd number of p^; thus it 
has (Pe+e-, Ptt) = (—5 — )• Accompanied couplings and Green functions can also be determined 



from Table |4.4| ; for example, a term of (Pe+e-, Ptt) = (~5 ~) is accompanied by either a2G*D 
or azG*F. But the former do not contributes to ddunpoi, since it is A^sm-A^edm- Since Pe'Pt 
CPT-even, it is accompanied by Re (asG*^); note that a CPT-even (-odd) term is proportional 
to "real" ("imaginary") part [Section [A.7.5|| . Among dcTunpoh only the part with Ptt = +1 
contributes to the total cross section atot [Eq- ( |A.51| )]; thus Utot oc ai|Gp. On the other 



hand, Pu = —1 part contributes to the Forward-Backward asymmetry A^b = (sgn(pe'Pt)) /c"tot 
[Section ^.3.2|| . Thus the "couphng" of the term p^ Pt oc Pt cos 9te is Re (asG*^). This completes 



the analysis of duunpoi- 

Next we consider the terms proportional to the spin- vectors: PfSt + PfSt. Since the 
polarization vectors Pf and are 3- vectors, each term of them is proportional to either of p^, 
p^ or PeXPf The number of p^ should be zero or one, since we are working to 0{l3t)] note that 
= E'^, where E is beam energy. Thus the momentum dependence of P^ and P^ is either 



"'^'^We are working at the CM-frame of e+e 
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Pg, Pj, (p^-p^), or p^y-Pf. These[^ are listed in Table |i75| . Also shown is the order of velocity; 
only the interference among SM contributions are 0{1), and the others are 0{(5t)- Note that 
under CPT, "eigenvalues" of (p, s + s, s — s) are (— , — , +) [Table [4.2|| . From this, for example, 
one can see that the chromo-EDM contribution to the CP-odd term|3 (sj — st)-{p^-p^) in da is 
proportional to Iie{a2G*D), since the term is CPT-even. This completes the analysis of the 
term PfSt + PfSt in duunpoi, except the relative magnitude of each term, especially of p^ and 
Pe {PePt)- Our result for (5P shows that they are proportional to either 

{St - St)-{Pt - PeTZrw) - St)-(PgXpJ , 



IP, 



e I 



which are both perpendicular to the beam direction p^. It seems that we miss some symmetry. 
However for the effect of chromo-EDM, it can be understood as follows. For that case, ti is 
produced in (1,3) = (0,1) or J^^ = 1~ at first. Note that P = (-1)^+^ and C = (-1)^+^ 
for fermion-anti-fermion pair. Then they are modified by the anomalous tig-EDM potential 
[Eq. (|4A0|) ], which has (P, C) = {-,_+); thus the resuhing state is J^^ = l+~ or (L, S) = (1, 0). 
On the other hand, (spin-l|(5'j — S'j)|spin-0) is zero for i = z, and non-zero for i = x,y. Thus 
there is no contribution to parallel component. 

One also see the terms proportional to St*s/ can be understood in the same manner, except 
the anomalous contributions: 

\Pe\ 

We can see that neglecting 0{(3t) and higher, 



d3pj 



(Xl+C/i{s// + Slj) + S/jS/, 



Thus if s/i = 1, which means the spin of t is parallel to Pg, then s// ^ —1, which means the 
spin of i cannot be anti-parallel to Pg, or should be parallel to Pg [Section |4.5.2|] . The situation 
is more clear for the case x = il? or C// = 



7/ 

dajst, St) 



oc (Its//) (Its 



This factorization shows that when x = +1 [~1]) o^' when the helicity of initial (e^,e^) is 
(L, R) [{R, L)], the spins of both t and i are anti-parallel [parallel] to Pg. For other choice of x, 
the spins of t and i have no definite direction. Including 0{l3t) corrections, the factorization 
do not hold even for x = ±1. 



4.7 Results for EDMs 

Here comes our numerical results. We shall see that it is possible to disentangle that the effects 
of each three anomalous interactions dtg, dt^, and diz- 



Here [p^-p^) is the vector p^ with the coefficient p^ Pf 

Of course there is no such contribution to the CP-even term such as (st + St)-{p^-p^). 
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If all the couplings dtg etc. are real and the decay width Tz of Z is neglected, the effects of 
anomalous interactions for polarization are 



(5Px = ( - Bldtg + Bldt^ + Bfdtz) Im A sin 9te 



5Pk 



( - B'^dtg + Bldt^ + Bl^dtz) Re ( G ) ^ ^^c?*^ 



A sin 6'te 



This assumption is realistic, since the imaginary part is generated by physical threshold in loop 
diagram. From Figures [4.11| and |4.9| , we can see that Chromo-EDM dtg and EW-EDMs (ij^, diz 



can be separated through their energy dependence, since typical j3t become larger for larger 
energy. The Re [{F — G)/G\ (3t becomes smaller for larger energy, because F = G for the case 
= 0; note that the effect of rescattering becomes smaller as one leaves from the threshold. 
On the other hand, from Figure [4.3| , we can see that the two EW-EDMs, dt^ and dtz, can be 
separated through the response with respect to the initial polarization x- This is also pointed 



out some time ago jO^] in the context of open top analysis. Note also that from the definitions, 
Bs_ = ±5n for X = ±1- 
The relation 



Ppeak 



V {E + 1 GeV + iVt 



(4.28) 



agrees qualitatively with Figure 
with V ^ 0, the relation Ppeak = 



.9. Here 1 GeV ^ 2mt — Mis = binding energy. Note that 



^/mWTi^t) hoidsQ. 

For stable quark under LO Coulomb rescattering, G and F can be obtained analytically p4 



lim 
lim 

E^p^ /mt 



E- 



where = CFOis'mt/2 ~ 20 GeV. Thus 



F 
G 



1 



P=Ppcak 



Ppoak 



^mt {E + 1 GeV + ^F^) 



This agrees qualitatively well with Figure |4.1CI| . Note that there is no energy dependence of 
Re [F/G] if Tf = 0. Thus this may provide a way to measure the decay width of top quark. 



4.7.1 Sensitivity 

The polarization (P) projected to certain direction can be measured as 



(P) 



This can be obtained from G = l/{p1/mt ~ {E + iTt)) for ^ 0, and da /Apt (x |Gp 
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Figure 4.3: Coefficients for polarizations as functions of initial polar- 



ization X- III the figures, C, 



para 



perp 



Cj_ and Cn 



(Left-top) Coefficients for SM contributions. (Right-top) Coefficients 
for anomalous EDM couplings. Lines in the lower two figures the 
same to those in the upper two figures: (Left-bottom) Coefficients for 
the polarization parallel to rv±. (Right-bottom) Coefficients for the 
polarization parallel to iin. 
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Figure 4.4: Green functions 
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Figure 4.5: Green functions 
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Figure 4.6: Green functions 
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Thus the error may be estimated roughly as S{P) ~ N^s . Here N^g means the number of 
events used for the analysis: 

Nes = (^tt X X [BfBh) X e 

= 0.5 pb X 100 fb"^ X f ^ ■ ^) X 0.6 = 4 X 10^ events, 



which means = 1.5 x 10~^. We used the branching fraction Be (Bh) of to or 

/x"*" (hadron) is 2/9 (2/3). We assumed integrated luminosity JC = 50 fb^^, and the detection 
efficiency e = 0.6. 
For example, with 

IT 



B^ Re ( dtg — — — ) Pt sin 9te , 



statistical error for dtg is 



1 

d(cos6'te) 1 



^"^'^ ^ 5^ X Re (dtg^) A f\._ , 



d (cos 6*46) sin 6te 

1 



eff 



1 4 

X - X 1.5 X 10" 



ttg-EDM = —dtg < 10"^ V cm 



1 X 0.1 TT 

~ 0.2 ~ C(10%) . 
From this estimate, we have 

gs 

mt 

for Chromo-EDM. Since all coefficients {C±, -Bn, etc.) are similar in magnitude and so are 
the ratio of Green functions (Im (F/G), etc.), bounds for electroweak EDMs are just the same 
order: 

tt7-EDM = —dt^ < lO^^^ecm , ttZ-EDM = —dtz < lO^^^fe cm . 
rrit 

A realistic Monte Carlo detector simulation based on the calculation here is also in progress 
by experimentalists [p.25|| . They show that high detection efficiency is possible with simple 
event selection criteria and 6-tagging. Up to now, only lepton energy asymmetry {E^+ — E^- ) 
for dilepton+dijet event is obtained. Statistical error with 100 fb~^ is 

^(stat) ^^^^^ _ _ Q_g4g _ 

Their results for anomalous couplings, based on our calculations, are 

I Re [dt^] I < 1.48 , I Re [dtz] \ < 1-04 , | Re [dtg] \ < 3.93 

to 95% CL. In fact, lepton energy asymmetry is not a good observable for polarization. This is 
because the time component of spin 4- vector is suppressed by A ~ 1/10 compared to the space 
component. On the other hand event asymmetries are not suppressed in that way, since they 
can be obtained by integrating 3-momentum asymmetry. Moreover the branching fraction for 
single-lepton + 4-jets is larger than that for dilepton + dijet. Thus it is expected that the 
sensitivity for anomalous EDMs is improved by factor 10 or more. This is consistent with the 
naive estimate we made above. 



Chapter 5 

Summary and Discussions 



We study how the ti threshold region serves as a probe of physics at electroweak scale and 
beyond. 

In Chapter ^, we studied the convergence of the perturbative series of the total production 
cross section crtot for e~^e^ —>■ ti in the threshold region. Two efforts have been made to 
improve it. One is to use Potential-Subtracted mass mP(/i/), and the other is to resum leading 
log(g//i<i)'s of the Coulombic potential Vc{q) in momentum space using Renormalization Group 
(RG). The former is to circumvent the (severest) renormalon ambiguity inherent to the pole 
mass and the coordinate-space potential of an IR-slave interaction. With this prescription, 
the convergence of the position of IS" resonance becomes better. The latter is to reduce the 
renormalization scale [fig) dependence of R ratio, where fig is the scale for the ti potential 
V. Among various contributions to the potential V = Vq + Vbf + H^a, only the log(g//Xs)'s 
in the Coulombic potential Vc are summed up in this work. RG improvement of the whole 
F is a subject of a future study. With these two implementations, we estimate {AmtY^ ~ 
0.1 GeV, where (Amt)*'^ is the theoretical uncertainty for the determination of rrit using the 
process e'^e^ —>■ ti with ^/s ~ 2mt. This size is consistent with one of the N^LO corrections 
calculated recently ||32|. Lower bound to the CKM element \Vtb\ may also be obtained. Note 



that Tt oc \Vtb\ , and the present bound [|| for \Vtb\ is 0.06-0.9993 when the unitary of CKM 
matrix within the first three generations is not assumed. 

In Chapter ^, we calculated the momentum distributions da /dp of top quark near the 
threshold. Corrections to 0{l/c^) are included. However the Final-State Interactions!^ between 
t and the decay products of i etc. are not taken into account. The effect of the NNLO correction 
for (l/(Ttot) da/dp from the rescattering between t and i, which is calculated here, turns out 
to be rather small compared to the NNLO correction to cxtot- In fact the correction is smaller 
than the effect of the NLO correction from the FSI. Thus the NNLO correction from the FSI 
may contribute with the similar magnitude to the NNLO correction from ti rescattering. 

In Chapter ^, anomalous EDMs of top quark are studied. They can be measured by CP- 
odd observables. Chromo-EDM modifies ti rescattering, while Electro Weak-EDMs modify ti 
production vertex. Thus their effects can be distinguished from the energy dependences, since 
the latter is enhanced by raising CM energy, while the former is not. On the other hand 
two EW-EDMs, ti'-f-EDM and tiZ-EDM, can be distinguished from the dependences on the 
polarization of the initial e~ (and e"^). We estimate the sensitivity of ti production near 
the threshold in future e"'"e~ colliders to be ~ 10"^'' (7^ cm (e cm, Qz cm) for tig-EDM {ti'y-, tiZ- 
EDM) with the integrated luminosity ~ 100 fb~^. It is worth mentioning that a realistic Monte 



^We call this FSI for brevity. 
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Carlo detector simulation is in progress by experimentalists ||125|| , based on our theoretical 
calculations. Their first result is already available, and our naive estimate is consistent with 
it. We note that hadron colliders may be powerful to prove the existence of non-zero EDMs; 
however precise determination of their values would be difficult. In this respect, e~^e~ collider 
is much advantageous. 
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Appendix A 



A.l Notes on group theory 
A. 1.1 Group theoretic factors 

Generators of a Lie group satisfy the following commutation relation: 

Trace normalization Tr and quadratic Casimir Cr for a representation R is defined by 

tr [T^T^ = TrS""' , T^T^ = CrI ■ 

For SJJ{N = 3) 

1 A^^ - 1 4 

7> = -, Ta = Ca = N = 3, = _-_ = -, (A.l) 

where F and A mean fundamental and adjoint representations, respectively. 



Coefficients Pq etc. of renormalization group equations are given in Section |2.3.4 When 
necessary, we use 

n/ = 5 , nH = l , (A.2) 
where the former is the number of light quarks and the latter is the number of Higgs bosons. 

A. 1.2 Color charge summation for color-singlet state 

For a particle- antiparticle pair ipip'^ in color-singlet state, the product of color charge is —Cp'- 

- 1 

(T'^f^) = -Cp = ■ (A.3) 

This corresponds to Q^-{—Q^) = — for Abelian gauge group. Sum over the index a is 
implicit. In this section, we show this relation. The relevant Lagrangian is 
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where = -T"^, and = Cip'^ with C7^"^Ct = -7/^. Then 



^CFtT[l]^-CF , 



where 



for example. The indices i and j mean the color of initial and final ip. respectively. 

Although the calculation above is specific to fundamental representation F, generalization 
to arbitrary representation R is obvious; (T"T") = —Cr, where Cr is quadratic Casimir for 
the representation R. 

A. 2 Notes on non-relativistic quantum mechanics 
A. 2.1 Relative coordinate 

As is well known, a non-relativistic two-body system can be decomposed into relative and 
center-of-momentum motions: 

2 2 2 p2 

H — + — — + 

2mi 2m2 2// 2M ' 

where 

„ mi 1712 _ 

R= ■ ri + ■ r2 , P = pi + p2, 

mi + m2 mi + m2 

m2 mi 

r = ri-r2, p = ■ Pi ■ P2 , 



mi -|- m2 mi -|- m2 



and 



— = 1 , M = mi -|- m2 (mim2 = /xM) . 

II mi m2 

One can easily check that those coordinates and momenta satisfy the standard commutation 
relations [Ri , Pj] — i5ij etc. Thus they are conjugate also in quantum mechanical sense. 

A. 2. 2 Radial momentum 

The operator pr for radial momentum is defined by the Weyl-ordered product of r/r and p: 

1 /I 1 
Pr = -[ -rp + pr- 

2 \r r 
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Using 

dr dr^ ' dr 



d ,9 d (r''\ n — 1 . , 
r—^r—, — ( - 1 = — — with n = 3 , 



one obtains 



. _d_ _ld_ / •^2__^ 2 _}_d_( 2^ 

P'^ Qj. r dr ^ dr'^ r dr r dr^ dr \ dr 

With these relations, one can show that this satisfy [r , Pr] — i, and ('0i|pr'02) = (Pr'0i|'02) for 
the states of hm^^o '^V'(^) = 0. Note that (V'ilV'2) = / / r^dripKr) ijj2{r). On the other 
hand, the angular part of is given by orbital angular momentum L — rxp: 

,2 2 



where 



d ( d 

1} = r^p^ - r-ir-p) p + lir p = r^p^ + — r^— 



Since [r , L*] = 0, I? and r also commutes. The following relations may sometimes be useful: 

(r-») = r^p^r^p' , r (r p) p = r^r^p'p^ — (r-p) +ir p — —r'^—-^ . 

One must be careful for the singularity at the origin r = 0. For example, with — A — p^ = 
for angle-independent function, naive calculation gives A(l/r) = 0, while actually 

A- = -47iS^^\r) . 
r 

The contributions of 5-functions are easy to miss. They arise from differentiation of functions 
singular at the origin. It is safer to integrate by parts so that to avoid these subtle points. For 
example, with functions of well-behaving, 

(/|P?;|S) = (P?/|;I«) 



= -/dSipdr{i|l(r/)}i9 
= 47r/(0M0)- /dQ / r'drf-—g, 



'0 



or 



p^l = 47r<5(3)(r)-i|^. 
We use following relations in the text: 

\pl,ipr\ = 47r5(3)(r) , [p\ip,] = 47r5(3)(r) + 3 



2L^ 



[iPr, -] = -J ■ 

The first relation is highly anomalous; p^ and p^ do not commute. 
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A. 2. 3 Operators in momentum and coordinate space 

A momentum independent potential V is diagonal in coordinate space: 

{x\Vix)\y)^Vix)5^'\x-y) , 

where hat " means operator. Such an operator in momentum space depends only on momentum 
transfer: 



{p\Vix)\q) = Jd'x{p\x)V{x) {x\q) 
= V{p - q) . 



Using 



{x\p\) = -iV^{x\) , {\p\x) = +iV^{\x) , ^f{x-y) = —^f{x-y) , 



we have 



dy' 



{p\[p, V{x)] \q) = (p- q)- V{p - q) , 
{x\[p, V{x)] \y)={- iV,V{x))5^''\x - y) 



dx' 



More generally, 
{p\VG\) = 

{p\pVG\ ) 

{p\VpG\ ) 



d^g 



(27r)3 
(2^ 



V{p-q)G{q) , 
pV{p-q)G{q) , 
V{p-q)qG{q) , 



{p\[p, V]G\ ) = J-0j-siP - Q) ■ V{P - q) G{q) 



{x\VG\ )^V{x)G{x) , 
{x\pVG\) ^ -iV,{V{x)G{x)) , 
{x\VpG\) ^V{x){-iV,G{x)) , 
{x\[p,V]G\) = {-iV,V{x))G{x) , 



where 

G{p)^{p\G\) , G{x)^{x\G\) , 
Our convention for Fourier transform is 



|f,(^|p)(p|/)-/^e-7(p), 



= = Jd'x{p\x) {x\f)^ Jd 



^xe-'^-^f{x) . 



For (5-functions, 



5('\x-y) = {x\y) = J^^ 
{27rf 5^^\p - q) = {p\q) = Jd^xe-'^-^-^ 



A. 2 Notes on non-relativistic quantum mechanics 



103 



A. 2. 4 S-wave projection 

If a quantity / depends on a 3-vector, say, r, it can be decomposed into "radial part" and 
"spherical part" . For example, a plain wave of momentum p can be written as 
Sometimes we use the shorthand notation for S-wscve projection: 



G 



1 



H-uj 



, uj = E + irt , 



Gir,r') = {r\G\r') , G{r,r') = J ^^G{r,r') . 



An An 



(A.4) 



L = 0, 5^ = 2, r^r^ = d'^r^S. 

In general, partial wave expansion can be defined as follows: 



f{k,9) = f^{2e + l)Mk)Pe{cos9) , Mk) = P,{cos9) f{k,9) . 



The second relation holds because 



Here Pi{cos9) is Legendre polynomials: 



J^P,{cos9)PAcos9) 



2i+l 



3 1 

Po(cos 9)^1 , Pi(cos 9) = cos ^ , P2(cos 9) = - cos^ ^ - 7; ■ 
Also the following relation holds: 

1/(^.^)1' = E(2^+ 1)1/^(^)1' ■ 

An example may be in order. A plane wave is expanded as follows: 

00 

gifcrcose _ Yy^^ + 1) i^u{^r) Pe{cos9) , 



where ji{kr) is a spherical Bessel, which is a solution Ri{r) — ji{kr) of the radial Schrodinger 
equation: 



dr^ r dr 



+ 1) 



Re{r)+k''Re{r) = , 



where k^/{2m.) = Explicitly, 



. , , sin p ... sin p cos p ... / 3 1 \ . 3 

Jo(p) = , Ji{p)^— , J2(p)=^-- smp-— cosp 

P P P VP PJ P 
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Fourier transform goes as follows: 

oo 

fix', x) = J](2£ + 1) Mx', x) A(cos e,>,) , 
f{p',x)^ y"dV e-'^^-^'f{x',x) 

oo 

= J](2£+l)/,(p',x)P,(cosV) , 
where 9x'x is the angle between x' and x, and 

PCX) 

fe{p', x) = / Anx"dx' {-iYje{p'x') h{x' , x) , 
Jo 



a;) = ^ ^^^^^3^ i^jeip'x') flip', x) . 



Here we used 



The following are completeness relations for the angular part: 

d{cos9)Pe{cos9)Pe'{cos9) 



- J^(2£+ l)P^(cos^) P^(cos^') = 5(^)(cos^ - cos^') , 



2 y_i ^ ' ^ ' ^ ' 2£ + 1 ' 
1 
2 

— / d</. e^("')0 ^ 5^^, , 

^ oo 

^ ^ e^m(,^-'^')^5(i)(0_0') 



^=0 

-.27r 



27r 

m=— oo 



Since 



/ d X a; — x') f{x', cos 6*', 0') = f{x, cos ^, 0) , 
(5-function in the polar coordinate can be written 
5(3) {x-x')^ ^5(^) {x-x') 5(^) (cos 9x, - COS ^,,,) 5(^) (0, - 0,0 , for / = /(x, cos 9, 0) , 
= ~ ^^^^ i^^^ ~ "^o^ ^^'2) , for / = /(x, cos 61) , 

= ^2^^'^i^-^') Jorf^fix). 
Here 6*^2 is the angle between x and certain axis, and likewise for (p^- 
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A. 3 Formulas for Dirac spinors 



References [|126|| are useful. 



A. 3.1 Mode expansions 

The anticommutation rules for the creation- annihilation operators and bp, are 

{a;, aj} = {6;, &?} = {27r)'2Ep6^'\p-q)6^' , (A.5) 

and all the others are zero. Here (b^J) is the creation operator for fermion (anti-fermion) 
with the momentum p and the spin s. A quantized Dirac field ip{x) can be expanded in terms 
of them: 

Since the anticommutator relations above are invariant under the phase rotation a — ^> e**^ a (and 
the similar for b), the relative phase between u and v is not fixed, a priori. Our convention is 
explained in Section |A.3.3. 



A one-particle state is \p,s) = a^J |0). Note that {p,r\q,s) = 2Ep5^^\p — q) 6"^" is 
Lorentz invariant. 

Note that all flavor- dependence is carried by the creation- annihilation operators and 6*; 
c-number spinors u{p, s) and v{p, s) is flavor- independent. 

A. 3. 2 Several conjugations: F^, f , F'' and F^ 

There are several conjugations^ for 7-matrices F; that is, hermitian conjugation F^, Dirac 
conjugation F, Charge conjugation F^, and the conjugation related to Time Reversal F^: 

f = 7°F^7° , F" = Cr'^C^ , F^ = B^V'^B , (A.7) 

where the Charge conjugation matrix C and the Time Reversal matrix B are defined by the 
following properties: 

C-t = C-i , = , B^ = B-^ , B^ = -B , (A.8) 
C^/^T^t = _Y , 5V*5 = +7° , B^-f'*B = -Y ■ (A.9) 

The last two relations can be combined to the single relation B'^'j^'^B = +7^^. For a product 
of 7-matrices, one can show that F1F2 ■ ■ ■ F„ = f „ ■ ■ ■ r2ri and (F1F2 ■ ■ ■ F„)^ = F^ ■ ■ ■ FgF^ 
Some of these are collected in Table [A.l| . By using these conjugations, the operations V, C and 



T for a spinor bilinear, which are defined in Eqs. ( [A.48|) , ( |A.49D and ( |A.50|) , can be expressed 
as follows: 

7017° = (f )^ , CT'^C^ = r , B^T*B = (F'')t . (A.IO) 



^ We limit ourselves to the spinors in 4-dimension. 
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+ 
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+ 


+ 
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Table A.l: "Eigenvalues" of 7-matrices under the several conjuga- 
tions defined in the text: f = j^T^'^, = CV'^C^ and = B'<T'^B. 
The symbol (—1)'^ is defined to be +1 for = and to be —1 for 
IJ, = 1,2,3. The index /i on (—1)^ is always not summed; besides 
there is no distinction whether it is covariant or contravariant. The 
notations = (—1)^7^ and a^*^ = (— 1)^(— l)''cr^'^ may sometimes 
be convenient. 

One can see that "Parity" and "Time Reversal" do not reverse the order of 7-matrices. For 
successive applications of the conjugations, one can show the following relations: 

F=(f)% P = (f)^ 

Also one can show {CB)'-f'^{CBy = — 7^^, which implies CB = 75 up to phase. Note that 
BC = {CBy. 

A Dirac conjugated 7-matrix F appears when one takes the hermitian conjugation of a 
spinor bilinear: 

(V^iFV^2)^ = V'2fz^i , (A.ll) 
for both anti-commuting and commuting spinors. 



A. 3. 3 Antiparticle field ?/;^ 

The antiparticle field ip'^ of a Dirac field if) is defined by 

i)\x) = Cifj'^ix) , = -^^(x)C^ , (A.12) 



where C is the charge conjugation matrix defined in Eq. (|A.9|) . It coincides with the charge- 



conjugated field CipC"^ up to field- dependent phase. The following relations can be shown: 



where ^l/r = Pl/r^P and tpL/R = ^Pr/l, where Pl/r = (1 =F 75)/2. Our convention^ for the 
relative phase between c-number wave-functions u and v are 



V = u 



Cm"^ , v = -u^C^ , = -uC 



^ In general, v — exp** u'^, which means u — exp"' v'^ . 
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See the next section for details. With this convention, an anti-fermion field ^p'^{x) is expanded 
in terms of the creation-annihilation operators as follows: 

P g 

One can see that a and 6 in ^ are interchanged in Thus the role of particle and anti-particle 
is indeed interchanged. 

A spinor bilinear can be rewritten in terms of their anti-particle fields: 

ipiVip2 = +'4'2^'^'^i 5 anti-commuting spinors. (A. 15) 

For commuting spinors, 'Uirf2 = — fjF'^'u'j; = —U2T'^Vi, for example. The second equality holds 
irrespective of the phase convention between u and v, since v = e^^ means u = e*'' v'^. 



A. 3. 4 Dirac representation 

For an non-relativistic particle, the time component of the 4-momentum is much larger than 
the space component. Thus it is convenient to choose the representation of 7- matrices such 
that the time component is "special" . Dirac representation is one of them: 

) , a^''^'-[Y,Y] , (A.16) 

7' = W^7'^7^7''7'^ = ^7°7'7'7' = 

where a* are Pauli matrices: 

These satisfy {7^ , 7*^} = 2g'^'^, where g'^'^ = diag(+, — , — , — ); and aV-^ = 5^^ + ie'^^^a^, where 
gi23 = _|_x_ Our convention for the 4-dimensional totally- antisymmetric tensor is that eoi23 = 
Note that the phase of 75 is chosen such that the chirality (= eigenvalue of 75) coincides 
with the helicity for a massless particle (or -u); for a massless antiparticle (or f), the chirality 
is the opposite to the helicity [Section [A.3.5|| . Note also that (75)^ = +1. Explicit forms of the 
charge conjugation matrix C and the time-reversal matrix B in the Dirac representation are 

C ^ ^ . B^,W-i^"\^). (A.17) 

Overall phases are not fixed by the defining properties in Eq. (|A.9|) . Our convention is such 
that CB = 7^ which means BC = {CBy = 75* = 7J. 

Freely-propagating spinor states are specified by the 4-momentum p^^ = {E, p) and the 
(3 dimensional) spin direction s (|sp = 1) at the rest frame. Note that s is not the space 
component of the spin 4-vector s'^ that is obtained by boosting (0, s)'^ to the "laboratory" 
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frame. However when the distinction between s and is not important, we write simply s. 
The c-number wave-functions are as follows: 



u{p, s) = VE + 



m 



pa- ^ 

-X 



pa- ^ 

-X 



\E + m^ J 



, v{p, s) = u%p, s) = VE + m\ E + m | ^ (A.18) 



where u'^{p,s) = Cu'^{p,s). The second relation defines the relative phase of u and v. From 
this it holds that v'^{p, s) = u{p,s). One can see that the upper (lower) two components of u (v) 
are 0{1), while the lower (upper) are 0{(3). This is the reason why the Dirac representation 
is convenient for non-relativistic particles. The expressions above can be obtained by boosting 
the expressions at the rest frame of a particle: 

u{p,s) = V2ml'^ ] , v{p,s) = 




^0 

Note that the mass-dimension of a spinor field is [ip] = 3/2, while that of a scalar field is 

m = 1. 



Other representations 

Representations other than the Dirac representation can be obtained by unitary transforma- 
tion0 U. For example, so-called chiral (or Weyl) representation 7^ can be obtained from the 
Dirac representation 7^^ as follows: 



7w 



Accordingly ?/^w = Uipr, and ipw = i^oU"^ ■ One can see that if 70 is hermitian (anti-hermitian), 
then 7w is also hermitian (anti-hermitian). As for the charge conjugation matrix C and the 
time-reversal matrix B 

Cw = UCy^U^ , 5w = U*By,U^ , (A.20) 

in general. These follow from the fundamental requirements C7^"'"C^ = — 7^^ and B^'j^'^B = 7^^. 
Since U for the Dirac-to-Weyl transformation is real, if one choose Cd = ^7d7d fo'^ the Dirac 
representation, then Cw = '^7w7w ^^^^ Weyl representation. However in general, the 

representations of C and B in terms of 7-matrix depend on the representation of 7-matrix 
itself. On the other hand, since the fundamental requirement for the parity transformation is 
AY^A^ = Y {A = 7°), which means Aw = UAdW, the representation of "parity matrix" A is 
independent to the representation of 7-matrix. One can also see that the relations = —C 
and B'^ = -B are representation independent. On the other hand, a relation such as = —C 
is representation dependent. 



^ The fundamental relation {7^, j"} = 2g^^l remains intact by the transformation 7 UjU^^ with 
invertible matrix U that is not necessarily unitary. However if U is not unitary, ^'^4' changes when one changes 
the representation. It may be convenient if 7** is either hermitian or anti-hermitian. Since (7^)^ = 2g^^ [no 
summation], 7'^ should be hermitian (anti-hermitian) if g^^ > (< 0) [no summation]. 
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A. 3. 5 Pauli spin spinors and helicity h 

Spin s 



Pauli 2-component spinors ^"^^ X 
at the rest frame of a particle, is 



for the spin direction^ = (cos sin 6, sin sin 6, cos 9) 



I ^ \ 
' cos - > 



Note that x"^"'^ 
the direction Sr: 



VA 

sm - e^f 
2 



X 



-^^'KX' 



( - sin ^ e"^"^ \ 



V 



cos- 



(A.21) 



-X^ . The X* cind x ^ cire eigenstates of spin operator for 



sr-(tx 



SR-(rx 



-X 



Thus x^^ surely represent the states of spin = ±S/j, and (1 ± sr-(t)/2 are the spin-projection 
operators for 2-component spinors. 

Spin 4-vector is obtained by boosting (0, s^Y at the rest frame of the particle to the 
"laboratory" frame: 



2 



Sr + 



7 



7 + 1 



(A.22) 



where (3 = p/E, (3 = |/3|, 7 = I/a/I — Note that {p^,s^) (p^^^s^) is equivalent to 
{p, Sr) {—p, that is, the Parity and Time Reversal transformations. 

The following relations are for the process e~^{pe) e~(pe) — ^ t{pt, St) i{pt, St) with -y/i = 2mt'y, 
where 7 = 1/ -^/l — /S^. We work in the CM frame of tt (and e"'"e~). We neglect the electron 
mass, and dropped 0{j3'^) and higher terms. Also note that St and St are spin 4-vectors at 
the "laboratory frame", which is the CM frame of tt, while Sf and are spin 3-vectors at 
the rest frame of t and t, respectively. Our convention for the totally antisymmetric tensor is 
eoi23 = +1- Use Pe + Pe = Pt + Pt and SfPt = SfPt = to obtain other relations. 



SfPe 
SfPe 



+ StiPt - Pe) 
+ StiPt + Pe) 



SfPe 
SfPe 



-Sf{Pt 

-^tiPt 



Pe) 
Pe) 



et.upa{Pe)^{pe)%Pt)'isty = -2mtStip p,) 

e^,upa{Pe)^{Pe)''iPt)^{St)'' = -2mtSf {p ,x p^) 



iPe)'{Pe)''is,nS,r = -2m,p^ 



StXS 



{Pe 



npt)''{st)^{s,r 



-pupcr 



{Pe)'{PtnSty{St^'' 



+mt{StXSt)-{Pt - Pe) + {{SfPt) St + {SfPt) St)iPeXPt) 

nit 

+nit{stxst)-{pt + Pe) - —{{sfPt) St + isfPt) St)-{p^xp^) 

rrit 



^ The subscript r emphasis that the sji is not the space component of the spin 4-vector in Eq. ( A.22 ). 
Since we rarely refer only to the space component of s^, we sometimes wri te sj; simply s. In fact, the difference 
between s and the space component of is O[0^^ , as can be seen in Eq. ( A.22 ). When the distinction between 
s and s'^ is not important, we may simply write them s. 
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Helicity h 

Spin degrees of freedom can also be specified by lielicityQ h = SR-p/\p\. The state of h = +1 
(—1) is called right-handed (left-handed) helicity state, for both particle and anti-particle. Thus 
the spin for helicity eigenstates are s = ±p/\p\ = ±f3/\f3\ for h = ±. Since s-crx'^'^ = 
we have 

\P\ 

where = with s/j = +p/\p\. By noting x^^"^^ = ^X^? Dirac spinors of helicity 
eigenstates are 



u{p, h = ±) = 4E 



m 



±\P\ 



X^ 



-x^ 



= ±) = I ^ + m'^ I . (A.23) 



Here u(p, K) = u'^{p, h), since charge conjugation do not flip neither momentum nor spin. For 
a massless state, we have 

uip,h = ±) = VE(^^^^^ , vip,h = ±) = VE(^~^l^ . 

With 75 in Eq. ( |A.16|) , we have 

75 u{p, h) = +h u{p, h) , 75 v{p, h) = —h v{p, h) for m = 0, 

which means that our convention for the phase of 75 is such that it eigenvalue (= chirality) is the 
same as the helicity for u (~ fermion), while the opposite for v (~ anti-fermion), for a massless 
case. Note that u (v) can be used also for anti-fermion (fermion), since uiTv2 = —v^T^Vi, etc. 
[Section [CTl . 

The spin 4-vector s'^ for helicity eigenstate is 

±1 / \ ±7 ff3''\ 7 fl 



7+1 , 

Thus s'^lh = ±1) ±p^/m as /3 ^ 1 both for a particle and an anti-particle. 



A. 3. 6 Lorentz transformation for 7-matrices 

Parity operation etc. are defined to act on ip's, but in Section |A.7.4| we show that equivalent 



operations can be implemented as manipulations on 7-matrices. Likewise, although Lorentz 
transformation changes only fields and (x^,p^, s^), equivalent manipulations for 7-matrices, 
which are constant to Lorentz tr., are possible. Consider a Lorentz transformation 

pf" A'^^p'^ , ^ Al^p . 



The helicity of massless particles may be defined by the limiting procedure m — > 0. 



A. 3 For mulcts for Dirac spinors 
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Since ■^■^ is a Lorentz scalar, ip — >• ■^Ai ^. Since ipp^'^^ip is also a scalar, we have A^j^Ai = 
Iv, or 

This means 7- matrices are "transformed backward" by Ai~^ and Ai. For example, let us 
consider the Lorentz transformation = A'^^{m, 6)", or {A~^Y — (m, O)*^. Thus we have 

Ai-VAi = (A-i)>'^7. = m7°. 

A. 3. 7 Projection operators for energy A±{p) and spin S(s) 

Thus energy projection operators A±(p) for positive and negative energy at the rest frame are 

. _ /l ^ _ 1 + 7° _ m + TO7° boost m + / 
A_|_ — 

A_ = 

or 



OJ 2 2m 2m 

\_l-7°_m - m7° boost m - / 



1/ 2 2m 2m 



A.(p)^i^. (A.24, 



It has a following properties: 



A+(p) u{p, s) = u{p, s) , A_(p) s) = v{p, s) , 
A_(p) it(p, s) = A+(p) s) = , 
A±(p) A±(p) = A±(p) , A+(p) A_(p) = for p^ = m^, (A.25) 
A+(p)+A_(p) = l , 
T^^±{p) — ±?77, A±(p) for p^ — rri^. 

Use the relation (■uirt'2)^ = V2TU1 for u{p,s) etc. 

Likewise, for the spin-projection operator E(s) at the rest frame, and the spin parallel to 
2;-axis (s^ = ±1), 

E(.) = '^^ = "-22? . (A.26) 



\ 2/ \ 2/^ ^ 

This has the following properties: 

E(s) s) = u{p, s) , S(s) s) = v{p, s) , 
S(-s) s) = S(-s) v{p, s) = , 

E(s)E(s) = E(s) , E(s)E(-s) = 0, (A.27) 
E(s)+E(-s) = 1 . 

At the rest frame, = (m, 0) and s'* = (0, s). Thus p-s = 0. From this, projection operators 
for energy and spin commute: 

A±(p)E(s) = E(s)A±(p) . (A.28) 
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Just the same manipulations are possible also for x'^ and x~'^- We show only the results. 
These can be verified also by direct calculations. For simplicity we mean by s, not the 
space component of in arbitrary frame. They are orthogonal x^'^X"^ = = 0. 

Similarly for u and f , 

s) u{p, s) 
v{p,s)v{p, s) 

As we saw above, at high energy, the spin 4-vector of a particle becomes parallel to the 
4-momentum of the particle: s'^{h = ±1) ±p^/m as /? ^ 1. Here h means helicity. Since 
(±y + m) = ±m (±/ + m) for = m^, we have for /3 — 1, 

u{p, h) u{p, h) = ^ ^^^^V ' ^(P' ^) ^(P' ^) = 2^^V 5 (A.31) 

where /i = +1 (—1) for right-handed (left-handed) helicity state. One can see that the relation 
between the helicity and the chirality (= eigenvalue of 75) is opposite between particle and 
antiparticle, or to be precise, u and v. 



2 ) A A 2 



1 - 



2 

1-/75 



+ m) = 2m A_|_(p) , 
— m) = — 2m S(s) A_(p) 



(A.30) 



A. 4 Gordon identities 

In this section, we show that the only CP-odd t-t-gauge boson interaction of dim-5 is that of 



EDM interaction; others are not independent. We closely follow the Appendix. C of ||123|| . Let 
us see ttZ interaction, for example. By integrating by parts, all derivatives on, say, can 
be removed. After this manipulation, dim-5 terms of effective Lagrangian can be written as 
follows: 



C ={gitPd>'^ - g2h^d^^ + ^g3i'cJ^'''^,^ + g^i,ia^'i^d,^)Z^ + h.c. 

Couplings gi^i can be complex in general; but should be real in order for those terms are odd 
under CP. In interaction picture, where the field operators satisfy equations of motion for free 
fields 

{i^ - mt)ij = , (□ + ml)Z^ = , d^Z^ = , 
the following relations hold: 

-{^id^^) + 2mV^7^?/' , 



A. 4 Gordon identities 
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In these relations, arguments of ip and ip are all x, and dy = d/dx'^. The first one is the 
celebrated Gordon decomposition: 

By using these "Gordon identities" , the Lagrangian above can be organized as follows: 

Formulas for u{p) etc. are obtained from 

(p, s\ ^{x) = (0| u{p, s) e'P-^ , ^{x) \p, s) = u{p, s) e^^"- |0) , 
{p, s\ jjj(x) = (0| v(p, s) e^P-^ , \p, s) = v(p, s) e-^^-^ |0) . 

For example, with the convention {p,p\ — p))^ — (ci^a| |0))^ — (0| apUp, 

{p,p\ i^Tid^^ij |0) = +{p - p)^' u{p)rv{p) e'^^P+P^^ , 
(p'l ijUd^'i/j \p) = +{p' + p)^ u{p')ru{p) e^^p'-p)-^ , 
(p'l i^Vid'^i^ \p) = + pYv{p)Tv{p') e^(p'-j^)-^ . 

Here p (p) means the momentum for a particle (anti-particle). Note that the last equation in- 
cludes the minus sign from the anti-commutativity of fermion. The exponential factor becomes 
4-momentum conservation for each vertex: 

J d^x e^(f e-^«-^ = (27r)^5(^) (p + p-q) ■ 

Thus for the process <— \q), 

iqi,u{p)a'^''v{p) = —{p — pY u{p)v{p) + 2mu{p)j^v{p) , 

{p — p)uu{p)a^^v{j)) = iq^ u{p)v{p) , 

iqy u{p)ia'"''j5v{p) = -{p - pY u{p)i-jr^v{p) , 

{p - p)u u{p)ia^'''y5v{p) = -iq^ u{p)i'^r^v{p) + 2mu{pY^'^^v{p) , 

where q — p + p. At the CM-frame of q^^ — {y/s, 0), 

_ _ _ — p) _ 

uip)Y^5v{p) ^ u{p) ia'^ '^^-^v{p) , 



zm Im m 



for example. 
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A. 5 Merging the effects of 7 and Z exchange 

In e+e^ ti, both 7 and Z are exchanged in s-channeL But their effect can be combined to 
single effective couplings. With the SM vertices 



where 



9j 
9 



g sin 6w 



,/7 



1/7 



cos 9 



w 



TzL-Qf sin^ 6 



w 



or 



V 



67 



V 

we have 



i7 



-1 

+2 
3 ' 



,tz 



i + sin^ 9 



ii-lsin^^ 



= -0.01876 
w = 0.09584 



,67 



3L 



(A.32) 








-1 
4 
+1 
4 



^(^^(Pe)Ax^M(Pe))(Mfe)r^f(Pi)) 

'—^ 7. X 



X=-i,Z 

s 



- [f ^a*] (vipehMPe)) {u{Pth^l5v{pt)) 

- [a%^] {v{Peh^l5u{Pe)){u{Pth^'v{Pt)) 

+ [a"a*] {v{Pehf,l5u{Pe)) {u{pt)Yl5v{Pt)) 



where 



[v%'] = v'^v'^ + d{s) v'^v'^ 
[v^a'] = v^^a'^ + d{s) if^a'^ 
[a^w*] = a^^v'^ + d{s) a^^v'^ 



^^a'] = a^^a'^ + d{s) a'^a 



eZ tZ 



d{s)v'^a'^ 
d{s) a^V^ 
d{s) a^^a'^ 



(A.33) 



are energy-dependent "couplings", and d{s) is a ratio of Z-propagator to 7-propagator: 



d{s) = 



9l 



s — m\ + iiTLzT z 



(A.34) 



where 



9l 



cos 9 



w 



{g sin 9wy cos^ 9w sin^ 9w (1 - 0.23124) ■ 0.23124 



5.625 = (2.372)^ 



Extensions for anomalous vertices should be obvious. The width Tz of Z introduces absorptive 
part: Tz/mz = 2.490 GeV/91. 187 GeV = 0.0273. For = 2x175 GeV, its relative magnitude 
is typically s/(s — m| + irrizTz) = 1.073 — 0.002 i. Thus for the threshold region, the effect of 
Coulomb rescattering overwhelms the effect of Tz- Also in the open top region, it is known |90 



that the QCD vertex correction is larger than Tz, as far as the normal component Pn of the 
polarization of top quark is concerned. 



A. 6 Propagators 
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A. 6 Propagators 



A. 6.1 Fermion propagator 

Let us define the momentum configuration at CM-frame of tt as follows: 



Pt='^+P , -Pt = ~+P ■ (A.35) 



It is sometimes convenient to use —pt instead of +pt, since it is the former that enters to the 
propagators. In other words, 



(A.36) 



spin 



spin 



The negative sign for v reflects anticommuting property of fermion. 



± — \- p + rrit = 2mt 
2 



E 



P-1 ^ 1 



2mtC ' 2mfC^ 

,2 1 



(A.37) 



2 ^ 2mt 2mtc^ \ 4 



Thus the Feynman propagators Spipt) = ^i^t + ^t)/{Pt ~ + ''-f^t^t) for t and t are 

■l±70 



5i.(±g/2+p) 



P-7 

2m+c 



2 ^ 2 



2mi 



|( + mt^j-^_+m£ _ 1 + 7° 1-7° 



2mt 
-7 75 



2mi 2 
rut 1 + 7° 



2mt 



-cr 75 



2mt 2 

-'^ + mt 1 + 7° 



icr-^75 

rritc 

-P 



O 



r 



+ o{- 



2mt 



2m, 



^75- 



rritc 



7 



+ 



(A.38) 
(A.39) 
(A.40) 



Note that the last two (and Eq. (|A.36D) are consistent with "Gordon identities" in Section |X74. 



A. 6. 2 Gluon propagator 

Gluon propagator {G'^{q) G'^—q)) = D^j.S"'^ is given by 

for covariant gauge, and 



/^.. = ^^-^.. + (l-0^) (A.41) 



for Coulomb gauge. The later is convenient for non-relativistic calculations. 
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A. 7 Parity, Charge Conjugation and Time Reversal 

Within QFT, Parity operation V etc. are defined on creation-annihilation operators. However 
sometimes it's convenient to consider those operations for c-numbers such as u{p)'y^v{p). 

For explicit calculations, it is convenient to introduce a shorthand notation ~ for flipping 
the sign of space-components of 4-vectors and tensors: 

= (-1) V , = (7^)^ , 

Relations to another notations are 



Oi 



(A.43) 



{an 



Note that the index on (—1)'^ is never summed over. Note also that 



A. 7.1 V, C and T for field operators 

Parity operation etc. for a Dirac field operator iIj{x) [Eq. ( |A.(j| )] are 

C^{x)C-' = vhVc*W'{x) , C^{x)C-' = VcVci^'ix) , (A.44) 

TiIj{x)T-^ = 7]^r]^*ilj{-x)B^ , TiIj{x)T-^ = t]tVtBiP{-S:) , 

where ip'^ is an antiparticle field defined in Eq. (|A.12|) . The symbol r/p etc. are flavor- dependent 
phases l^^pp = 1, and rjp etc. are flavor-independent "over-all" phases \i]p\'^ = 1. In terms of 
creation-annihilation operators, a^, and bp, these transformations are 

ValV-' = +VPVpa-p , = -VpVp*b% , 

Ca;C-' = vcVcK ' ^K^'' = Vhvc*^; , (A.45) 

Note that c-number spinors u{p, s) and v{p, s) are flavor- independent, and 

u{p, s)7° = +u{p, s) 7°m(p, s) = +u{p, s) 

v{p, s)7° = -v{p, s) 7°^(P, s) = -v{p, s) , 

u{p, s)"'"C^ = —v{p, s) Cu{p, s)""" = v{p, s) 

v{p, s)"'"C^ = —u{p, s) Cv{p, s)""" = u{p, s) , (A. 46) 

{u{p, s)YB = u{p, s) B'^ {u{p, s))* = u{p, s) 

{v{p, s)yB = v{p, s) B'^ {v{p, s))* = v{p, s) . 



A. 7 Parity, Charge Conjugation and Time Reversal 
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In order for CVT (and the permutations) to be a symmetry of an action, these phases cannot 
to be arbitrary: 

r]pi = ±1 , rjTi = V*PiV*cH ■ 

If one (or more) of fields is Majorana, the reality condition ip = ip'^ restricts the "over-all" 
phases: 

Note that C and V commutes if and only if rj*p^rpp* = —rjpirpp, orQ rpp = ±i. Note also that T 
is anti-unitary: 

(/I T-'OT \i) = {T{f)\ O \r{i)Y T{cO)T-' = c*Tor-' . 

= (r(z)|ot|r(/)) , (A.47) 

Only spinor-bilinears are relevant for scattering matrices Ai, because they are Lorentz scalars: 

v[M^)T{x)Mx)]v-' = v*PiVP2 P(r) Mx)T{x)M^) 
c[Mx)r{x)Mx)]c-' = vhivc2 C(r) Mx)r{x)Mx) , (a.48) 

T[i)i{x)T{x)tp2{x)]T''^ = r7Jir7T2T(r) tpi{x)T{x)ij2{x) 
Here the parity eigenvalue P(r) for F etc. are defined through the following relations: 

7°r(x)7° = P(r) T{x) 

CT^ix)C^ = C(r) r(x) (A.49) 
B^T*{x)B = T(r) r(-5) 



These eigenvalues are summarized in the Table [A .21 . For example, = —7* [i = 1, 2, 3) 

means P(7*) = —1, and 7°^7° = —gf^z^ means P{di) = —1. Note that, since charge 

conjugation interchanges fields, even though C{d^) = +1, be sure that C changes to (9^, 

since (9^)^ =d^. Thus C(a^) = -1, where AdB = A{d -d)B = A{dB) - {dA)B. Note also 
that T(zr) = -T(r), or T(e*'^r) = e-2^'5T(r), in general. 

In the momentum space, eigenvalues for Parity etc. are defined as follows: 

7°r(p, s,p, s)7° = P(r) T{p, s,p, s) 
CT"^(p,s,p,s)Ct = C(r) Tip,rs,p,s) (A.50) 
B^r{p,s,p, rs)B = J{T)Tip, ~s,~P, ~s) 

For example, 7°(«75(p — p)*)7o = il5{—{p — PY) means P{i'j5{p~py) = +1. This is in accord 

with P(i75i9*) = +1. From the definitions above, we have the following for Dirac conjugates 
F = ^^T^^: 

P(f ) = P(F)* , C(f ) = C(F)* , T(f ) = T(F)* . 



This follows since one can choose riPi = +1 for a certain species i. 
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Table A. 2: Eigenvalues of 7-matrices and derivative. They are not 
transformed under parity V etc., which are defined to operate on 
creation-annihilation operators. It should be understood that they 
are sandwiched by ^ip and if). However, one can define similar trans- 
formations for c- numbers. See the text for the definitions. Note that 
= ie^'^^Vp, with eoi23 = +1. Although ij(d^ + 8^)^^ = 
d^{'ip^), this is listed because each of d and d is not the eigenstate of 
C. The symbol (—1)^ is defined to be -|-1 for = and to be —1 for 
H = 1,2,3. The index on (—1)'^ is not summed always; besides there 
is no distinction whether covariant or contravariant. 
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Table A. 3: A gauge connection should transform the same way to 
id^. For charge conjugation, extra minus sign comes when changing 

d^i to by integration by parts. As for time reversal, (seemingly) 
extra minus sign comes from the anti-unitary nature of T. 



P(r) and C(r) are (usually) real, but T(r) may be complex depending on the choice of the 
phase of T: T(e^^r) = e-'^'^J{T). 

All entries in Table |A.2| are chosen so that F = F. Note (9^)^ =<9^. For a product of 
7-matrices, FiF2---F„ = f „ ■ ■ ■ f sf 1. Likewise for F^ = CT^C\ (F1F2 ■ ■ ■ F„)^ = T^ ■ ■ ■ F^F^. 
Thus, be careful to the eigenvalue C of a product; it is not just the product of each eigenvalue, 
in general; however 7-matrices and derivatives commutes, of course. On the other hand, V and 
T do not reverse the order of matrices. 

We also summarize the transformation properties of gauge fields in Table [A.3|. For non- 



Abelian gauge, = A'^^T"- and F^^ = F'^JT"'. For V and T, the arguments x of fields change 
to X and — x, respectively. For Charge Conjugation, CA^C~^ = —A*^ = —Aj^. In general, field 
dependent phases should also be considered, as in the case of ip. But we do not go into detail 
here. 
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A. 7.2 CVT symmetry 

From the Tables |A.2| and |A.3| , one can see that a Lorentz invariant interaction such as ip'f'^ijj 
is CPT-even: CPT = +1. It's easy to remember, since one can assign —1 for each 4- vector 
index /x; for example, CPT-eigenvalue of a tensor such as a^'^ is (—1)^ = +1. 



A. 7. 3 V, C and T for currents 

When the topology is the same for all the diagrams contributing to we can apply V, C, or 
CV to, say, an initial current and a final current, independently. A (tree-level) amplitude can 
be written as 



M = J2 (^^1 itTAt)iZZ){eABe) \ee) 



A,B 



J2 (iTAt) |0) (0| (ZZ) |0) (0| (eA^e) |ee) 



A,B 



We can rewrite, say, the final tt current ju using V etc.: 

{tt\ (tTp=±t) |0) = {tt\ (ptp)(tTp=±t)(ptp) |o) 
= ±(P(tt)|(tTp=±t)|0) , 

{ti\ (tTcp=±t) |0) = {ti\ {{CVycV){iTcp=±mCVycV) |0) 
= ±{CVmiiTcP=±t)\0) , 



where 



{ti\ (tTT=±t) |0) = {ti\ (TtT)(tTT=±t)(TtT) |0) 
= ±(T(tt)|(tTT=±t)|0)* , 

\V{p,s,p,s)) = - \{p,s,p,s)) , 
\C{p,s,p,s)) = \{p,s,p,s)) , 
|T(p,s,p,s)) = \{p,s,p,^)) , 

\V{p,s,p,s)) = {Va^p,s)V^){Vb\p,s)r^){r\0)) 

= - \{p,s,p,s)) , etc. 

Due to the complex conjugation for T, its not a good idea to consider only the final or initial 
current but the whole amplitude. 



because 



A. 7.4 P, C and T for amplitudes or currents 

Although the parity transformation V etc. in the QFT are defined for creation-annihilation 
operators, the similar manipulations can be applied directly to c-numbers: Combined with the 
definitions Eq. ( |A.49|) [or Eq. ( [A.50| )] for the transformations of 7-matrices, we can analyze V 
properties etc. of matrix elements Ai solely in terms of c-numbers. 
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Parity and Charge Conjugation 

We consider a subprocess Z*{q) t{p, s)i{p,s) for definiteness. Denoted in the parentheses 
are their momenta and their spins. The relevant vertex can be decomposed into pieces, so 
that each of which is an "eigenstate" of parity operation: 

r'^ = + + ■ ■ • 

where = A{q'^)'y'^, for example. Here Alq^) is a form factor. Hereafter, we consider a part 
of the tt current j^^ = Ja + is + ' ' ' • 

J^b' P> ^) = s)T'^v{p, s) 

= (n(p,.)7°)(7°r:;7°)(7°^(p,^)) 
= -P{Ta) ■ j^(-p,s,-p,s) , 

where the minus sign can be understood by the fact that the product of the intrinsic parity of 
a particle and that of its anti-particle is —1. 
From the result above, or 

j'a{p, p)jb\p, p) = P(rA)P(Fi,)* • -p)jri-p, -p) , 

we know that the terms in that are odd (even) under {Pt,Pt) {—Pt, —Pt) should come 

from the interference of the ti vertices with the opposite (same) parity, P(Fa)P(Fb)* = — 1 
(+1), and vice versaQ. 
We can also see that 

jdnfj^fs* = for P{Ta) P(Fb)* = -1 . (A.51) 

This can be understood as follows. For a fermion-antifermion system, its parity is P = (— l)^"*"^. 
Then different parity eigenvalue means different orbital angular momentum L. And the states 
with different L are orthogonal each other: 

-Fl{cosO) Fl' [COS 



J An ^ ' ' ' 2L + 1 ' 

which means partial waves of different L are orthogonal each other. Note that an eigenstate of 
partial wave is a spherical wave, not a plane wave that is a momentum eigenstate \p). In other 
words, \p) is a superposition of the states of definite L. Thus before phase-space integration, 
they do not form eigenstates of definite L, and then need not be diagonal with respect to L. 
We can further operate Charge Conjugation: 

f^ip, s, p, -s) = -P(F^) ■ {vip, sVC^) {CiT^^yC^) {Cu{p, sY) 
= +P(F^)C(F^)-j^(-p,8,-p,s) 

= CP(Fa) -7^(^,5,^,5) , 

or 

j^(5, 5)j- *(5, s) = CP(F^)CP(Fb)* ■ fAs, 5)j- *(5, s) . (A.52) 



We may sometimes write P{rA) as P{j)tu and P{Ta) P{Tb)* as P{jj^)tt- 
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Table A. 4: Transformation properties of momentum and spin. De- 
noted in the parentheses are for when p + p = 0. 



Thus the terms in |A1p that are odd (even) under {st,St) (st, St) should come from the 
interference of the ti vertices with the opposite (same) CP, CP(r^)CP(rB)* = —1 (+1), and 
vice versaPI. We can also see that 



E 



fAfe' = for CP(r^) CP(rB)- = -1 . (A.53) 



This can be understood as follows. For a fermion-antifermion pair, CP = (—1)"^^^, which is a 
compact expression for the fact that spin-0 (spin-1) state is antisymmetric (symmetric) under 
the interchange of s and s. Explicit calculation may be useful: 

$^( A-')^ (x^Vx"^") =^tr [(X-V')(XVV'" 

s,s s,s 

= tr [a'] = , 

where we used the completeness relation for spin states: X^X^'' = 1- They need not be 
zero before spin-sum over both s and s. With the definition x"'^ = —ia'^lx'^)* and the explicit 
calculations, one can easily see that x^^X~^ is a spin-0 state, and x''^cr*X~'^ is a spin-1 state. 

It is easy to see that these analysis can equally be applied to initial current, say, e"*"e~ — > Z*. 
We may use a term "tt-current" Parity transformation Pu for [p^^p^) — > {—Pti—pt)i "e^e~- 
current" for (Pe)Pe) ~^ {~Pe: ~Pe)y cind "over-all" for to flip all the momenta; and likewise for 
CP and C. 

The results of this subsection can be summarized as follows. First we express Ai* M.ji which 
is a part of |A1p, in terms of momenta p and spins s. Then let us define Parity operation in 
terms of p and s as summarized in Table |A.4|. Now the results in this subsection reads 



Mp^^Mp=± ■■■ P-even , 

Mp^^Mp=^ ■■■ P-odd , 

A^c*p=±-^cp=± ••• CP-even, 

Alcp=±-^cp=T ••• CP-odd. 

These results are very plausible. One might think that the same situation holds also for Time- 
Reversal. However in the next subsection, we shall see that this is not the case. 



We may sometimes write CP{Ta) as CP(j)tt, and CP(rA) CP(rB)* as CP(jjT) 
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Time reversal and T 

As we saw, an initial current and a final current can separately be rewritten using either Parity 
and/or Charge Conjugation operations. This is because Charge Conjugation operation invokes 
transpose in the spinor space, while each current is 1 x 1 in that space, which is invariant under 
transpose. As for Time Reversal, it is convenient to consider whole Ai, because each current 
is not real. Of course Ai itself is not real, in general. We mention to this point later. 

For definiteness, we consider the process e~{pe,r) e~^{pe,f) Z*{q) — > t{pt, s)i{pt,s) at 
the tree level, and decompose the final vertex and the initial vertex A^ into pieces that are 
"eigenstates" of Time Reversal. Let A{q'^) and B{q^) be the form factors for the final vertex, 
and C{q^) and D{q^) for the initial vertex: 

V'X^A-V'X, Ac,^C-Ac,, 

and 

Mac = AC -Mac (A.54) 
= AC ■ u{pt, s)T'^v{p, s) v{pe, f)Ac^u{pe, r) , 

up to factors etc. Then 

M^ciP, s) = {u{p„ sTB) (5tf^*5) [B^vip, sY) 

X {v{p,,frB){B^Ac*,B){B^u{pe,r)) 
= J{tA)J{Ac)-MAc{-p,-s) , 

or 

MAciP,s)MBDiP,s) 

=T(f^)T(fB)*T(Ac)T(A^)* ■ MAci-p, -s)Mi^{-p, -s) . 
Here (— p, — s) means to flip all the momenta and the spins. Now 

\M\' = \Mac + Mbd + ---\' 

= ■■■ + MacMbd + MXcMbd + ■■■ 

= ■ ■ ■ + Re {AB*CD*) [MacM^u + MIcMbd] 

+ I Im {AB*CD*) [MacM^b - MXcMbd] +■■■ 

means that the terms in \}A\'^ that are even (odd) under (p, s) {—p, —s) should be accom- 
panied by Re{AB*CD*) when J {V a)^ iV b)*^ {Ac)'^ {AdY = +1 (-1), or by Iyii{AB*CD*) 
when T(ryi)T(rB)*T(Ac)T(A£))* = —1 (+1). Note that one can always adjust the phases of 
f's and A's so that T(f^)T(fB)*T(Ac)T(Az5)* = ±1. Correspondingly, the phase of AB*CD* 
changes counterwise. 

The statement above can be simplified as follows. First let us write T(r^)T(r5)*T(Ac)T(A£)) 
as J{MM*). Since MM* is Lorentz scalar, CPJ{MM*) = +1 [Section |A.7.2|| , which means 



T{MM*) = CP{MM*). As we saw in the previous subsection, CP{MM*) is equivalent to 
the eigenvalue under CP-transformation (s, s) {s, s) [Table |A.4|| . Thus the result in this 
subsection is that 

(MM*) 

cpf=+ °^ Re(product of relevant form factors) , 

{MM*) 

cpt=- °^ Im(product of relevant form factors) , 
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where T-transformation is also defined in Table |A.4|. Note that T do NOT interchange initial- 



and final-states, while the genuine T-transformation do. 

In fact, the argument above holds only for tree- level diagrams. This is because we neglected 
ie in propagators at Eq. ( |A.54| ). The sign of this imaginary part becomes relevant once we across 
a physical threshold in a loop: log(— x ± ie) = log(x) ± in, for example. Thus, to be precise, 
the argument above holds when there is no absorptive part. As we shall see below, the sources 
for absorptive part are rescattering and physical threshold, which are absent to tree level^ If 
there is an absorptive part. Re {AB*CD*) and Im above are replaced by Re (^AB*CD* e**^) etc., 
where e*'' is an effect of absorptive parts. 

Symmetry considerations for |A^p continue to Section [4.3.1. 



A. 7. 5 Absorptive parts and CPT 

Transformation properties of S'-matrix and T-matrix, S" = 1 + iT, is as follows: 

VSP-^ = S pTV'^ = T if VHV'^ = H 

{CV)S{CV)-^ = S {CV)T{CV)-^ = T if {CV)H{CV-^) = H 

TST~^ = TTT-^ = if THT-^ = H 

{CVT)S{CVTY^ = {CVT)T{CVT)-^ = if {CVT)H{CVT)-^ = H 

The last two rows can be understood by S* ^ e~*^* and T is anti-unitary. 
In terms of matrix element, CVT invariance can be written as 

Tf, = {f\T\€) 

= (/I {CVT)-\CrT)T{CVT)-\CVT) \i) 
= {CVT{f)\T^ \CVT{{))* 
= {CVT{i)\T\CVT{f)) 

where i (/) denotes CPT-conjugate state of a state i (/). CPT is defined in Eq. (f4.5|). Thus 

T% = tA = tJ. , from CVT invariance, 

fi if n ' 

= Tfi , if there is no absorptive part, 

where the absorptive part of a T-matrix is defined by its anti-hermitian part[^ T — T'^ (or 
times —i or something). The relation above Tl = Tjj, or |T?-p = |Tjjp, for the case of no 

absorptive part has the following consequences: an expectation value of a CPT-odd observable 
is zero, unless there is an absorptive part, while that of CPT-even observable can be non-zero 
without an absorptive part. This is the basis of the classification of observables based on CPT 
transformation |l89|,p3|. In fact, as we saw in the previous subsection, CPT-even (-odd) term 
is proportional to the real- (imaginary-) part of the relevant product of couplings. 



^Decay width is a part of physical threshold effect, but can be treated to the tree level. Likewise, the 
imaginary part of form factors originate from physical thresholds, but those are treated by effective vertices 
here. 

Hermitian may be called "dispersive part" . 
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The dagger ^" on T (and S) for T complicates symmetry arguments. However T is hermitian 
to tree level; that is, no absorptive part. This can be shown as follows. Unitarity S'^'S = 1 in 
terms of T-matrix is -i{T - T^) = T^T, or 

-t{Tf,-Tj^) = J2T}nTm. (A.55) 

n 

This is a celebrated formula, which relates the absorptive part of an amplitude and the "cut 
diagrams" of the amplitude. If an intermediate state n is a two- (or more) particle state, the 
RHS is loop effect. If n is a one-particle state, it hits the pole of the propagator for n: 

1 



while r„ ~ Im [self-energy of n] is again loop effect. Thus anyway, absorptive part (T — T^) 
comes from loop effect. 

A source of absorptive part can be traced back to the boundary condition for Green func- 
tions: a positive energy particle propagates forward in time, while a negative, backward. This 
boundary condition is embodied by the Feynman prescription for the position of the pole of a 
propagator: 



p2 _ ^2 _|_ 



If one choose principal- value prescription 



1 1 1 



X 2 \x + ie X — ie^ 

instead, there would be no absorptive part. We can see this by analyticity of an amplitude 
A4{s) as a function of CM energy ^/s. Below the threshold, there is no absorptive part, and 
thus M.{s) is real: 

M(s) = [M(s*)]* . 

Analyticity requires this to hold also above thresholds: 

Im[M{s + ie)]^ -lm[M{s-ie)] . 

Thus, no absorptive part for principal- value prescription. We know by experience that the sign 
of ie in propagators are irrelevant to tree-level. This is in accord with the observation above 
that no absorptive part to tree-level. 

A word may be in order for "tree- level" . Higher order effects can be treated by effective 
interactions, and so are absorptive parts. They can effectively introduced as anti-hermitian 
parts of Lagrangian. For example, a decay width F can be introduced to "tree" level: 

£ = (m-i^)^V'H , (A.56) 



2^ 

which is consistent with s — + imV ~ s — (m — iF/2)^. Thus in the statement "no ab- 
sorptive part to tree-level", anti-hermitian parts of Lagrangian are counted as higher order. 
It's interesting that it seems QFT with the Feynman prescription specifies the direction of 
time-flow: 

k/,|2 \ -lEt \2 I -i(TO-ir/2) |2 -Ft 

\yj\ ~e ~e^ '-^ ^e 

That is, simple time reversal t —t may not be a symmetry, if absorptive part is finite. This 
may be a "conjugate statement" to that CPT and CVT differs when absorptive part is finite. 
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A. 7.6 CV violation 

Roughly speaking, CV transformation changes an operator to its hermitian conjugate, while 
the coupling intact [Eq. ( |A.57| )]. T changes the coupling to its complex conjugate. Thus CVT 



changes a term in Lagrangian to its hermitian conjugate. 

There two kinds of complex phase in M. : one is weak phase e and the other is strong phase 
5. Non-zero sin 5 means non-zero absorptive part. In this sense, 5 can also be called "absorptive 
phase". On the other hand, e is related to the hermitian contribution; thus it might also be 
called "dispersive phase"; it is also sometimes called "CP phase". Both of these phases are 
important for CP- violating phenomena. For example, let us consider the amplitude Af for a 
certain process with a final state /, and its CP-conjugated amplitude Aj: 

Af = Ai e'^' +A2 e'^^ 

= \Ai\e'''e'^' +\A2\e''^e'^^ , 

A J = Al e'^' +A*2 e'^' 

= \Ai\e-''' e'^' + \A2\e-''^ e'^^ . 

Here Aj is determined up to overall phase, which is irrelevant to Now one can calculate 

the following CP-odd quantity: 

_ 2lm{AlA2)sm{5i- 62) 

\Af\^ + \Af-\^ ~ |Ai|2 + + 2 Re{AlA2) cos((5i - S2) 







-2\Ai\ 


\A2 


sin(ei 


- 62) sin((5i - 82) 


1^1 


2 + 


A2 


|2 + 2|Ai| 


\A2 


cos(ei - €2) cos((5i - ^2) 



One can clearly see that one needs both of two independent CP phases and two independent 
absorptive phases. 

The fact above that one needs an absorptive part, is well known from K physics. However 
before sum-over final state momenta and/or spins, one do not need an absorptive part for CV 
violating observables. 

In general, one can not say that a certain interaction C^^t violates CV. This is because of the 
degrees of freedom to redefine CV by flavor-dependent complex phases. For flavor-off- diagonal 
interaction, vector-like phase transformation changes CP-property; while for flavor-diagonal 
interaction, chiral phase transformation does. A simple example is a Dirac mass term written 
in terms of two chiral fermions: 

-Crn = m'ipip = m ipn^L + m ^/jl^r . 

Here non-zero imaginary part of m means (effective) absorptive part, not CP- violation. In 
terms of Dirac fermion, this is flavor-diagonal; while in terms of chiral fermion, this is flavor- 
off-diagonal. One can easily see that this term is CP-even. Now we make chiral rotation, which 
is a symmetry of the Lagrangian if it's anomaly-free: 



Accordingly the mass term above changes its form: 

—Cm mijj e^"'"'^ ip = m cos a tptp + m sin a ij^ip 
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Now the term i'j^ip is CP-odd, while ipip is even. Thus the mass term above violates CV 
symmetry, although the original one do not. In terms of chiral fermions, CV interchanges ipR'4'L 
and ipL'ipR with the transformation [Eq. ( |A.44| )] 



{CV)^l/r{CV)-' = Vcpl'{^L/Rr , 

with[^ r]cp = rfprfc = 1. With ijjiTip2 = 4'2^^'^i and {■iIjiTiIj2V = '4'2'^4'i, one can see that 
the non-zero complex phase a seems to mean CV violation. In fact, as can be expected from 
the context, this can be cured by taken into account of the degrees of freedom to choose the 
field-dependent phases rjcp in the definition of CV transformation: 

{CV)MCV)-' = e+-70(^z.)'= , {CV)MCV)-' = e-^'^j^^nY . 

One can see that Cm can indeed be a CV eigenstate. Usually one redefines the field to include 
the phase {r]cp)~^^'^- For the case above, 

{CV) e-*"/2 ^p^iCV)-^ = 7°(e-^°/2 ^ (cp) e+'"/2 ^^(CP)-^ = ^%e+'''/^ ^j^f . 

Note that if)'^ contains i})"^ . The redefinition (e~*"/^ t/^^,, e"*"*"/^ il)^) (ip^, ipji) takes the Lagrang- 
ian back to the original form, where the CV property is transparent. To summarize, non-zero 
imaginary part in the coupling for a flavor-off- diagonal interaction may or may not be a source 
of CV violation; it depends on how many explicitly broken U(l) symmetries are thereQ. 
While for flavor- diagonal such as 

their CV property is definite, since the flavor- dependent phases rjcp cancel in each term. One 
can unambiguously say that a certain flavor- diagonal interaction is whether CP-even or odd. 
Some more examples may be useful. Let us consider the following S — P interaction: 

{tjjtp - ^ i'-f^'4))(l) . 

From Table [A.2|, one can see that tptp is CP-even, while ij^ip is odd. Thus one can not define 



the CP property of so as to preserve CP. The situation is similar for the V — A interaction, 
which violates both C and P maximally. However there is a difference: S — P interaction is 
off-diagonal in terms of chiral fermion, while V — A is diagonal (and conserves CP). Thus 
in principle, there is a possibility to "rotate" ip so that the interaction above preserves CV. 
However such a degree of freedom may already be used to rewrite the mass term, as we did 
just above. Thus if is the Higgs that gives mass to the above S — P interaction may be 
"rotated" to the CP-even interaction iljilj(f) at the same time. On the other hand, if is not 
the Higgs that gives mass to ip, the interaction of (p may violates CP, in general. 

Likewise, MDM and EDM interaction are also flavor-off-diagonal in terms of chiral fermion: 



"'^^Here we dropped field-independent phases r/pi]^, which do not contribute to any fermion bihnears. 

^'^ Unbroken U(l) symmetries, such as fermion-number or electric charge, can not be used to redefine the 
phase of couplings, since those symmetry do not alter the form of Lagrangians. On the other hand, axial U(l) 
symmetry is explicitly broken by a fermion mass, thus can be used to absorb the phase of the fermion mass. 



A. 7 Parity, Charge Conjugation and Time Reversal 
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Since chiral rotation of t is already fixed by the mass term of t, non-zero d means CV violation. 
Non-zero imaginary part of d means (effective) absorptive part, in this case. 

Hermitian conjugation of a fiavor-off-diagonal interaction is closely related to the CV oper- 
ation. For example, 



{cv)c{cvy 
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t 



;-i)w^;* + h.c. + - 



(A.57) 



7! " 



- 1-7^ 



t 



+ h.c. + 



where tilde~means to flip space components [Eq. ( |A.43| )]. Since the gauge coupling g is real, the 
term above is invariant under CV with the phase convention that Vtb is real. CV transformation 
of gauge bosons are summarized in Table |A.3| . 



Effect of rephasing 

We saw that CV transformation of £int looks different for the different choice of the phase 
of fields. Let us consider how M. changes under the phase transformation of the fields. The 
key observation is that C\nt remains intact under field-redefinitions, since couplings are also 



redefined to absorb the difference. Since 



only (/I and \i) 



f 



,t...^t 



(exp(i / £i„t) - l) 



rephasing affects 



a' |0). This means each term in A4 changes with the same phase. 



Thus each term in AiAi* is rephasing invariant; thus it is sensible to consider the real or 
imaginary part of the product of the couplings in CPT-argument. 



Appendix B 



B.l Top quark decay width Ft 



Decay width of top quark is calculated in 

T{t ^ bW+) = TBorn 

where 



3 TT 



Born 



tb\ 



w 



Sttv^ mt 



nit 



-\ 2 



+ 



1 + 



mt 



mt 



mw 
mt 



Sttv^ 



1 - 



m 



2 \ 2 
w 



mt 



1 + 2- 



m 



w 



mt 



where 



Pw = W momentum in the t rest frame 



^J ml - (mw + nibY^J m^ - (mw - / {2mt) 



and 



/ = [tt^ + 2 hhiy) - 2 Li2(l - y)] + [%(1 -y- 2y^) Iny + 2(1 - yf{h + Ay) ln(l - y) 
- (1 - y)(5 + 9y - Qy^)] 1 2(1 - y)\\ + 2y) , 

where y = m^/ml. 

B.2 Coulomb plus potential: explicit calculation 

Confluent hypergeometric function F(a,7;z) 

Confluent hypergeometric function F(a,7; z), or Kummer's function, is defined by 

^A;!(7)fc 7 217(7 +1) 

where (a;)^ = a{a + 1) ■ ■ ■ (a + A; — 1), (a)o = 1. This is one of the solutions of 

d?w / ^ 7\ dw a ^ 
dz^ V z) dz z 
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Its asymptotic behavior is 

F(a,7;C)^l, forC^O, 

Note also that 7; () = e~^ F(7 — a, 7; — C). 
Whittaker's differential equation is 

d^W f I k fi^- l/4\ 

+ -T + ^ + ' ]W = 



One of the solution can be written by F: 

W = g-c/2 F{fi-k + i, 2/i + 1; C) ■ 

Coulomb plus potential: explicit calculation 

This section is continued from Section |2.4.2| . 

Let us consider the case E < for the moment. The case when E > can be obtained 
from analytic continuation. The parameter u is pure imaginary for < 0, thus we introduce 
another variable C: 

C = !£ 

V 

By multiplying Eq. ( p.36| ) by (— z^^), one obtains 

d^ \ iv n 



By comparing this and Whittaker's differential equation, one can see the correspondence k = 
—iu, /i^ — 1/4 = —K, or /i = ±a/1/4 — K = d± — ^. Thus the solution[| is 

9(0 = C'^ p^^^ ^ 2rf±; C) , (B.3) 



or 



g±{z) = z''^ e^/(2p) F(d± + p, 2ci±; -z/p) , p = , ;^ = -pC • 



Any solution of Eq. (|B.2|) is expressed as a linear combination of g±, and so are g>{z) and 



g<{z). Since F(a,7; C) — ^ 1 for C — ^ 0, we have 

g^{z)=g+{z) . (B.4) 



^ For Coulomb (k — 0), these two degenerate, and independent solutions are 

C F(l + zz., 2; C) , C e"'^/' + zz., 2; C) In C + (1 + ii., 2; C)] 

where F* is some other function. Thus it may not a good idea to let k = at this stage. However, wave 
function etc. can be obtained with this method. 



B.2 Coulomb plus potential: explicit calculation 
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While the limit C oo is more complicated. For < 0, Re^ < 0, since 



CpOts / rn II- 

■ ■ — "jIz^I + e , or p ~ + «e . 



2 \ E + ie 

Thus one can see the the first term in Eq. ( JB.ID dominates: 

r(d±-p)Vp; ^ r(d±-p) 

With this expression, we have 

g^i^z) = g4z) - p^-^+ r(^^-)^(^+-/^) ^^(^) . (B.5) 

The Wronskian W can be calculated from g±{z) ~ z'^'^ for 2; — 0: 



W = l-2d_ = Vl -4k , (B.6) 

Now all ingredients for the Green function in Section |2.4.2| are obtained. 

Next we want to obtain the location of pole and its residue F(a,7; () diverges when 7 is 
zero or negative integer. However since d± > 0, g± do not diverge. Thus the divergence of the 
Green function comes solely from T{d^—p). The poles sit at d^ — p = — (n — 1) (n = 1, 2, 3, ■ ■ ■ ), 
or 

p^ = {n — d-Y , while 
~ 4 E + ie ■ 

Thus 

4(n - d^y 

Especially with k ^ or (i_ — 0, we have energy levels for Coulomb potential 



y^F^s III, . inn \ 

^n = -— T^-«e (n= 1,2,3, ■■■) • 



{CFasfm (Ci7a,m/2)2 1 

= — = 9 = J-,2,3, ■ ■ ■ j . 

One can also read off energy eigenfunctions from the residue of the pole, since 

5 



{E + iT) 



where the wave function for bound states are properly normalized: (ipnl'^m) = ^nm- With 
n' = n — 1 + n' = n — dJ) 



(-1)"' ^, (-1)"' {CFas) 



m 



^^"^^-P^^ d^-p + W n'\ n'\ 2{n-d.fE-E^ + ^e ^ near the pole. 



since with p{E) = 

dp 2{n — d. 



d\E\ {Cpasym 



near the pole. 



132 



Thus 



-If (C^«,)2m 



r(2rf_) 



or 



2 (n - r(2 - 2rf_)r(2rf_ -n)E-En + ie 

{-if Cpasm 1 r(2d_) 



1 



47r(l - 2rf_) n'\ 2 - rf^)4-2<i- r(2 - 2(i_)r(2ci_ - n) 
X ^ i , while, 



= rr'G{r, r') ^ 
which means 

CpdsTn (— 1) 



E - En + ie^ 



E - E„ + ie 



r(2(i_ 



1/2 



[ Stt (n-l)!(l-2d_)(n-rf_)4-2'i-r(2-2(iV(2rf--ri)J ^^^""^ ' 
^+(z) = 2'='+e-^/(2p)F(l-n,2d+;z/p) , p = n - c/_ , c/+ = 1 - . 

Here we used 7; Q = F{j — a, 7; — ^). For Coulomb, using r(e — n) = - '■"V , 



we have 



.2 / 1 



r^n(O) 



-2 -2 , -, 
1 - 7^ + 1 + 



where 



IV^n(0)| 



2 12 V 2^2 



2 {CFasm/2f 1 



These are for < 0. For £" > 0, one obtains 



ImG(r,r') = ^iV'uWP , where 

= 1 - iCi.a,m(r + r') + 0{r^) , 



4^ in 



ur 



where 



/(x) = Re [^(^a;) + Ofe] = X] C(2n + 1)( 



n=l 



and Tq'^'^^'' is given in Eq. (|2.45| ). Note ImGo(0,0) = m^M/(47r). The factor {ipu 



1 — e~^^ 



u 



B.3 Results for Coulomb potential 
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is called Sommerfeld-Sakharov factor. One can see that for both E < and E > case, 



Im Gc{r, r') ocl - ^Cpasmir + r') + 0{r^ etc.) , 

irrespective of E. 

The Sommerfeld-Sakharov factor is derived as follows. Consider Im G(0, 0) for Coulomb. 
By expanding to Taylor series, one have 



g>{z) = 1 — zlnz + 



+ 0{z',k) , 



where ip{z) = r'{z)/T{z) is digamma function. Thus 



G{r, r') 



Cpagrn^ ( 1 



An z 



-In 2; + 



2 z zp 



+ 0{z) , 



or 



ImG(0,0) = ^q^Im 



This can be reduced more. Let us concentrate for E > 0. Using 
2 (Cfcks)^ m 



4 E + ie' 



I TrCpag E 

or p = h e , where u — \ — 

ZTT u \ m 



and 



~2~z ~ "^^^^ " ^ + i + 5Z C(2n + 1) 



2n 



n=l 



one obtains 



Thus we obtain 



TT 



1 TT 1 + e 

Imlnp = — , ih(—p)—i h real 

^ 2 ' 2p ^' 2 1 - e-^- 



Im (7(0,0) = 



for E>0 . 



4 1 - e-^« 47r 1 - e"^- 
The second factor in the last expression is the Sommerfeld-Sakharov factor. 



B.3 Results for Coulomb potential 

Here we summarize some results for Coulomb potential —CpOLs/r. Noting that the reduced 
mass p is m/2, we define 

Bohr momentum • • • pb = Cpasmj^ , 

Bohr radius • • • re = 1/pb , (B.7) 

Bohr energy • • • E^ = p^/m . 
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Bohr energy may also be called Rydberg energy. Typically pb — 20 GeV. See Figure p.4|. For 



^<0 (n= 1,2,3, ■■■), 



and 



where 



£^ - En + «e 



V-nlr) = Vn(0) e-^/(^Bn) _ 2; 2r/(rBn)) 

2 



While for E > 0, 



im? E 



ImGc(r,r') = ^A/-|^ij(0)|2 
47r V 



1_ I Jl + ZL ) +C'(r2 etc.) 
re tb ' 



l^i^(0)P 



1 -e-^B ' y E 



Thus irrespective of energy 

ImG'cf'", r') oc 1 



— + — ) + Cfr^ etc.) 

Tb Tb/ 



With this relation, we have 



1 d 



rGc(r 



masT dr 

where Gcij-) may either be Gc{r,r') or Gc{r,p). 



ma^r I V c 



(B.8) 



(B.9) 



B.4 Expressions needed for matching 

As was explained in Section |2]^, perturbative expansions with respect to and with respect 
to q;s//3 are both valid when -C /? -C 1. Here we collect both of them. These are used 
to determine the matching coefficients for NRQCD calculations. In the following, we use the 
quantity u that is equal to /5 to the leading order. Relation between them are as follows: 



u _ I E / ^/s — 2mt 
c V rritc^ y rrit 




B.4 Expressions needed for matching 
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4m^ 1 1 u'^ 



,2 



1 — 



7^ - 1 = «M 1 + X ) • 



B.4.1 NRQCD calculation of R ratio 



When cts <C /3 <^ 1, the Green function G defined in Eq. ( |2.26| ) can be expanded perturbatively: 



10 



lmG{r,r') = Im Gi{r, r') 



where with 



1 1 
E — le 



m 
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lmGo{ro,ro) =Im(ro|G'o |ro) 



ImGi(ro,ro) = Im(ro| Gq 

m^u CpasTi 
47r u 2 

ImG2(ro,ro) =Im(ro|Go 

-M\2 5 



1 + 



r L 4nc 
2 a. 1 



{2/3oln(/.V) + ai} Gq |ro) 
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rri^u Cpds — H 



-llTrGiraj 



47r 



c 6mroC 



+ C(ro) , 



Go |ro) 



ImG4(ro, ro) = Im (ro| Gq 



Goiro) 



rn^u 
An 



c ymroc c 2 



ImG5(ro,ro) = Im (ro| Gq^^^— ^Gq |ro) 



hi 



An & 
ImG6(ro,ro) =Im(ro|Go 
nn?u {CpasY n^ 



2mr^c^ 
( 2'murQ 

CEO'S ^ Cpd, 

-(jQ- 



An V? 12 
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■Goko) 
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Thus we have 
47r 



lmG{ro,ro) / 



u 



y u 2 \ u / 



TC 

12 



+ 



a., \ as „ I 



-Cf 
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ai — 2/?o In 



2'mu 
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/ 2murr 
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asTJiro 
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asTJir 



U\2 51 a, r CTT 
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u mrQC 



e-- J+^-i^^' + i 



(B.IO) 



The two expressions above are exactly the same each other. The R ratio and ImG is related 
in Eq. ( |2.8| ). By comparing the result in this section and that in the next section, one can 
determine the matching coefficients Cj'^"'^^ and C2^^^\ For this purpose, it is convenient to 
rewrite u with To each order in a^, their relations are 



u 



One can just substitute j3 for u in other cases, since there difference is and higher. 



B.4. 2 Perturbative QCD calculation of R ratio 

•> 7* — s> gg was done in |H0|. Here we summarize 



Perturbative calculation of R ratio for e~^e 
the result for reference: 



R=^NcQ'i/3{l-^ 



1 + Cf 



TT 



TT 



A(2) + 0(«3) 
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where m is mass of g, and 



where 



A^,^^ = — — - 2— + — + ttM ln/5 + -ln2 + — 

^(2) 3l7r2 2 /^179 , ^ 8 , \ 151 13, 

A^;^^4 = +7rM In/?- - In 2 U , 

72/3 V 72 ^ 3 y 36 2 ' 



Stt^ 11 

is^^y ' 

44 47r2 



(2) 



9 9 ' 



aS/9=-— /5oln2/?, 



4/3' 



B.5 Top quark polarization 

Neglecting C(/5t) and higher, the production cross section for ti pair can be written as follows: 



da ^ Nca^V t 1 - Pe+Pe- ,^,2 



,4 



X |G|'(ai + xa2) (B.ll) 



Here we sum over the spins of t and t. 

D^G-F (B.12) 



|G|^(ai + xa2) - ( [no Sj, St dependence] + [st dependent] 



4 



+ [st dependent] + [s(, St dependent] | (B.13) 



[no Sti St dependence] 

= 1^12(01 + xa2) + 2 Re [G*F(as + xa4)] 



PePt 

ml 



|G'|'(ai + xa2)<ll + 2Re(CFB^ ) Acos^te}> (B.14) 



B.5 Top queirk poleirization 
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[st dependent] 



rut 



- Re [G*F(a4 + xas)] 
+ Im [G*F{a^ + xai)] 2 



St-Pt ^ St P^ p^ Pt 
rrit rrit mf 

Sf(PeXPt) 



+ [ - luv{dtgG*D){a^ + xa2) + Im [G'*F(a5 + x^e)]] 
+ [ - Re(cit,G'*L>)(a2 + x«i) + Re [G*F(a6 + x«5)]] 
= |G'|2(ai + xas) |c;?s// + 2 Re (^C10 cos Ot, 
+ Re l^d^ j sin 6**6 + Im /^Cn^ j snA sin 6**6 



g^-(PeXPt) 



+ 



+ 



Bnuiidtg^ +Im r^^^ 



sj_/3t sin Ote 
s^(5t sin 



(B.15) 



[st dependent] = — > Sf, djg — > — c?tg, dt-y — > — c?t-y, d^z — — c^tz at [sj dependent]^ (B.16) 
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Note that d 



dt-y and dtz 



-d^ mean B^/^ 



-BJ/^ and B^J^ 



-B. 



7/Z 



In general, 



[sti St dependent] 



St Pe St Pe 



^\G\^{ai + xa2) 

rrit rat 

+ Re [G*F{a^ + xcia)] 
-Im[G'*F(a4 + x«3)] 
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mt mt rrit "n^t 
+ [ - Re{dtgG*D)(ai + xas) + Re [G*F(a5 + xag)]] 

'SfPeStiPeXPt) SfP^Sfip^XpJ 
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rrit 



rrit 



\G\'^{ai + X«2) + 2 Re {^^^ "^//"^//A ^ 

+ Re ^^N^^ (•S//SL + s//s±)Asin6'te 
+ Im ^Q-^^ (^//^N + s//Sn)A sin dte 
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S^Im ( dta^ 1 +Im 
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Rc I dtg- I +Rc I Bx- 



(s//sj_ - sj_s//)/3tSin6'te 
(s//Sn - snS//) A sin 9te 



{St-Pe) {St-Pt) - {St-P^) {St-Pt) = {StXSt)iPeXPt) 

By replacing Pt with {PgXpt) we have 

(St-Pe) [StiPeXPt)] - (St-Pe) [Sf{PeXPt)] 
+ {StXSt)-[pt\pf - Pe (PePt)] 

= {Pe)i{Pe)j{Pt)k{St)l{St)m ' -^Siljeklm] = 



(B.17) 



(B.18) 



(B.19) 



Thus one can see that the anomalous contributions to the [st, St dependent] -part is proportional 

to StXSt- 
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